{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Tolerancing - Compensators"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "\n",
    "from optiland import optic\n",
    "from optiland.tolerancing.core import Tolerancing\n",
    "from optiland.tolerancing.monte_carlo import MonteCarlo\n",
    "from optiland.tolerancing.perturbation import DistributionSampler"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "class Singlet(optic.Optic):\n",
    "    def __init__(self):\n",
    "        super().__init__()\n",
    "\n",
    "        # add surfaces\n",
    "        self.add_surface(index=0, radius=np.inf, thickness=np.inf)\n",
    "        self.add_surface(\n",
    "            index=1,\n",
    "            thickness=7,\n",
    "            radius=19.93,\n",
    "            is_stop=True,\n",
    "            material=\"N-SF11\",\n",
    "        )\n",
    "        self.add_surface(index=2, thickness=17.0695)\n",
    "        self.add_surface(index=3)\n",
    "\n",
    "        # add aperture\n",
    "        self.set_aperture(aperture_type=\"EPD\", value=25.4)\n",
    "\n",
    "        # add field\n",
    "        self.set_field_type(field_type=\"angle\")\n",
    "        self.add_field(y=0)\n",
    "        self.add_field(y=10)\n",
    "        self.add_field(y=14)\n",
    "\n",
    "        # add wavelength\n",
    "        self.add_wavelength(value=0.48613270)\n",
    "        self.add_wavelength(value=0.58756180, is_primary=True)\n",
    "        self.add_wavelength(value=0.65627250)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgUAAAFfCAYAAADeYb8BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/TGe4hAAAACXBIWXMAAA9hAAAPYQGoP6dpAACSjElEQVR4nOzdd1hT1x/H8XcmewgkIAgCshUFgrNurdqhdVStnXZX7bKtVru3be1utbu1vy47rN3b0daqVYIoywUqLgh7j4z7+yNMRUUZATyv5+FR5Cb3BJX7yfee8z0ySZIkBEEQBEE478ltPQBBEARBEDoHEQoEQRAEQQBEKBAEQRAEoZYIBYIgCIIgACIUCIIgCIJQS4QCQRAEQRAAEQoEQRAEQailtPUAWstisXDs2DFcXFyQyWS2Ho4gCIIgdDqSJFFaWoqvry9y+anrAV0+FBw7dgx/f39bD0MQBEEQOr3Dhw/Tq1evU369y4cCFxcXADIyD9T/XhAEQRCEBqWlpfQJDjrjdbLLh4K6WwYuLi64urraeDSCIAiC0Hmd6Ta7mGgoCIIgCALQzqHg77//ZvLkyfj6+iKTyfj222+bfH3u3LnIZLImH5MmTWrPIQmCIAiCcArtGgrKy8sZMGAAK1asOOUxkyZN4vjx4/Ufn3/+eXsOSRAEQRCEU2jXOQUXXXQRF1100WmPsbOzw8fHp8XPWV1dTXV1df3nJSUl5zw+QRAEQRAa2HxOwcaNG9FqtYSHhzNv3jzy8/NPe/yyZctwc3Or/xDLEQVBEAShbcgkSZI65EQyGWvXrmXq1Kn1f7Z69WocHR0JCgoiIyODBx54AGdnZ7Zs2YJCoWj2eZqrFPj7+2PIzROrDwRBEAShGSUlJWg1XhQXF5/2WmnTJYlXXHFF/e+jo6Pp378/ffr0YePGjYwbN67Zx9jZ2WFnZ9dRQxQEQRCE84bNbx80FhwcjJeXF/v377f1UARBEAThvNOpQsGRI0fIz8+nZ8+eth6KIAiCIJx32vX2QVlZWZN3/QcOHCApKQkPDw88PDx4/PHHmTFjBj4+PmRkZLB48WJCQkKYOHFiew5LEARBEIRmtGsoSEhIYMyYMfWf33PPPQBcd911vPnmm+zatYuPPvqIoqIifH19mTBhAk8++aSYMyAIQqeUX5nPH+m/YyyrIdQhDJVcdcpjIyIicHJy7sDRCULrtWsoGD16NKdb3PDbb7+15+kFQRBaJa8yj0SDHr1Bj96QwMGSg/VfU1qUeFd507OyJz0re+Jgcaj/msViISsri2nTpttg1IJw7rr8hkiCIAhtJa8yF70hsTYINISAQNdA+qhD8MnzYYDnACLjItEX6tlesB19iR4JiVCXUOI9BjLQYyCHth/EZDTZ9sUIwjnosD4F7aWkpAQ3NzfRp0AQhLNWFwL0hgT0OXoOlR4EIMg1iDitDp1WR5xWR/HRItas+Rp/f3+GjxiBXN4wR7vEWIK+wBoQdhQkUmGuwElyIloezXMzluOkcrLRqxOEBi3tUyBCgSAI543citz6WwF6g56s0kOANQTotPHovHXEauLwcvCqf8zBAwf46qsv8fX1ZeSoUU0CwYlMFhNpJWmsSfqaFHkK7g7u3DngLi4KvPiMW9YKQnsSoUAQhPOeNQQk1H4k1oeAYNdgdN46dNp4YjVxeDp4giRB8WHkhlRkuenIDGmYj+1CXnyIfJUfiqhLqAoYRU2PcDjDBX7D+vVUq6vZ57efPw//wQCvGBbpFhPhEdERL1sQTiJCgSAI5x1DhaF+PoC1EpAFQLBbH3TauIYQgByZIQ1Zbjry2l9luenIasoAkOzdqHEPISUXqhx86O9RjdPxLciN5ZgcvanoNYJy/5FU+QxCUtqfNI4N69djZ2fPzFmz2J69jeWJyzlQnMn0kBnM6z8fdzv3jvy2CIIIBYIgdH85FTkNqwNyEjhcdhioCwE6dF4D0Clc8Sw6iiw3DbkhHZkhFVlZNgCSQo3kGYakjUTSRCJpo7Booii2OPLxx/9DpVIxYeJE1Go1mI3YGxJxOvw3jkf+RlV6GIvCnsqeg6joNZKKXiMwO3kDTUMBgMli5Kt9X/F28lvIZXLm9Z/P9D4zUMib3+NFENqaCAWCIHQ7ORU56HOsVYBEg74+BPRx64POPZyBCjfiqo145mcgM6QjK8hAJpkBkNwCsGgjkTRR9QFA8ugDiqa9BiorK/n4fx9RU1PDxEmTcHR0PHkgkoSq5BCOR/7G8fBf2BuSkElmqj3Cqeg1kn+yHSlxCWXm7CuaPKygqoAVO9/g+8zvCHUPZZFuMbHauPb5ZglCIyIUCILQ5WWXZzfpE3Ck7AgAIY6+6FSeDDSa0RUZ8Mzdg6ymHADJ3r3Ru/5IJG1fJE0E2Lmc8Xwmk4nVn39Gbm4uky66qMU/U+TVJTgc22wNCUc2oagpoUjdE8cZK5CCRp50fGp+Ks/rnyM1P4WJvSdxV8zdaB21Z/GdEYSzI0KBIAhdTnZ5dv18AH3Odo6WHwMgVOFCvFnOwOI84oty6GGxNC3915b9JW0UOPuccSJgcyRJ4rtvv2Xfvr1cOGECGo3m3F6ExUT6rx8wsPQ3tFUZmCOmYBr3BLj1anqYZOHHAz/wRtLrVJorubHvTVwZfhVqhfrczisIpyFCgSAInV52+XEScrajP/I3O3J3cKSmEIAws4z48lIGVVYSV1WNu0svLNooa+m/7haAR/BJpf/W+GvjRrZs2cyoUaMI6N27Vc+1Yf167NR2zO4rR7n+MagqwTzsbsxDFsAJExPLakp5O+Vtvtz7BX7OftwTdx/DfYe36vyCcCIRCgRB6HSO56WjP/AziTl69GUHOSpVARBeXUN8VTXxFiVxrsG4aaIb3f9vWem/NZKTk/npxx+Ii4ujb79+rX6+JhMNq0tR/Psiim1vgWsvTOOfxBI66aRqRkZxBi/on2d7znZG+I7gnrh78XcJaPVYBAFEKGiVSa9t4mhRVZs8lyCcj9TU0IejBKj2oHLcS6m9gYMOVWSr5MgkidAaE32qVLhXaDBV9uGQJZi9UgA59AA6tsmPJElgNhHqamZafCDDAl1xULVuV/kTVx8AyPL2ofzjAeQHNmDuMx7z+KeQPENOGsv6I+t5OfFF8qvyuTriam7oexMOSocTTyEIZ0WEglb4NukYZdWib7kgnJFkwaXqGD3K9uNRto+a8nQOGA+QqihDb2fHUZUSmSQRZLEjVO6Dv30UHq4XYHSKRJLbfuuVqsoqtm7dAioHDGof9uRWYaeUMTjAmTF93Bjs74z9OQSE5kIBAJKEfO8vKP98CEqPYx48H/MFC0HddDfFKlMlH6Wv4qO0j3C3c+eu2IVMCJgguiIK50yEAkEQ2lZFPjJDGvLcdGS5aWBI41jBPhKUFhLs7UhwdOSYQo4MCLfTEOfVnzj/scT6DsPNzs3Woz+JyWTik08+pqS4hIsvuRgHBweOldSwMaOEjZkl7M+rwl4pY2hvF8b0cWWQvzNqZcsCwilDQR1jJYqtr6PY8ho49MA09nEsUdNOuqVwtOwoL+94iY1HNhCn1bEobhGhPcJa+9KF85AIBYIgnBtTFbK8PdZ1/nUNf3LToCyHo0oF2x2d2e7mRYJKxnGpBhkywlyD0fUcYu0YqI3FVd35/y/+/NNPpKamMHHSJDw9PU/6+uGiav7KLGFjRgmZBdU4quQM6+3C6D6uxPs7oVacOiCcMRTUKcpC+efDKPb+hCVgGKYJy5C0fU86bOvxLSzXL+dwWRaXh8zktv7zusT3WOg8RCgQBOH0JAsUHaq/6NeFAGvDHwsScMSjN9s8erJdrUJvKiTbWIIMGeE9wus3EIrRdI0Q0NjOpCR++eVnhg0bRp+QkDMef6jQGhA2ZJRwqLAaJ7Wc4YHWgBDn54xK0fQdfotDQS3ZgY0of1+KrCADs+5GzCPuBwf3JscYzUa+2Luad1PeQaVQsaD/7UwJvkx0RRRaRIQCQRAanFD6t/b934PMWNvwx8EDiyaCLM8gEuzt0FvKSCjNJKfSgFwmrw8Bcdo4YjVxuKjbdzVAe8rOzubj/31EcJ8+DBky5Kwff6Cgio0Z1oBwpLgGFzs5I4JcGRXsSqyvE0qF7KxDAQDmGhTb30GxaTko7DGNeQjLgKtA1rQikVeZy+tJr/PTwR+J9IhkkW4x/b0GnPXrEM4vIhQIwvnIWIksf6/1Xb8h1RoCDGnIyg0ASAo7JK8wa7Mfrwiy3HxIoAJ90T4Sc/XkVOQ0CQE6rbUS0JVDQGNVVVWs+vAD5HIFky6ahEJx7u+yJUkis6C6NiAUc6zEiKu9gpFBLriXZBLmLmf27LMIBXXKslGufwJFypdYesZimvgcku/JrZB35u5kuf45dhfu5pLAS7kj5s4mWz4LQmMiFAhCd1Zf+k+r3+1PZkhDVpiJTLJYD3EPbFjrr43C4hXJYbUKfV5S/f4BhtpKQESPCOsGQtp4YjQxOHeTENBYXcfCjIz9XHLppbi4tN1rlCSJfXlV9ZMUs0uNeNuZeefG4UT4nNt5ZIf/Q/n7EuQ5yZgHXIVp9EPg1LTLotli5vvM73hj5xuYLEZu7ncLs8OuQNWGTZ2E7kGEAkHoLsrzmkz4s27z27T0L2kisWj7NnT704QjqZw4XHYYfU5C/f4BDSEgsjYE6LptCDhRUtIOfv3lF0aOHEnvwMB2O48kSXz88yZ+zXMn36jkrrEh3HBBIAr5OSwntJiRJ/0P5V/PgMWMecT9mONvhBOWc5bUlPDWrjf5ev9XBLgEcF/cIob0HNpGr0joDkQoEISuxlhpnfVf+67/pNK/0t5a+q9t9Wup/RUnb5DJkCSpPgTUbSCUW5nbNAR411YCVM5nGEz3kpeXx6oPPyAwKIihQ9v/Yrlh/XoUKnuOeMTw/r8HiQ/owbPT+9Grxzk2IaooQPnXM8h3fISkibCuUuh9civkfYV7Wa5/nsTcREb3GsPC2Hvwc/Zr5asRugMRCgShs5IsUHiw4aJvSLNWAAoP1M76l0GPQOsOf3Vb/GqikDyCmrxDlCSJrNKshg2EDAnkVeahkCmI8Iisvx0wQDPgvAsBjZnNZj5a9SHV1dVcdPHFqFTtX1pvPNFw+8EC7v8mhaJKIw9fHMHUGN9zbkIky96J8velyI9swxx5mXWjJdemF31Jkvg963de3fEyxTXFXBt5LddFzsVedEU8r4lQIAidwYmlf0OatRpgrABAcvBs8q6/rvR/Yoc7sP6wP1R6yBoCcvQk5urrQ0CkR2Tt6gDdeR8CTrRxwwa2bfuPSRdd1Gw/gvZw4uqDsioTT/28m7VJx5gQpeXxyVF4OJ3jboiShDzlK+tGSzVlmIctxDx4PijtmhxWYazgw7QP+GT3x3jae3J37D2M8x8nuiKep0QoEISOZKxopuFP+qlL/9q+SJpIcNKecptfawg4SEKOvn5OQH5V0xCg0+oYoInBSeXUka+2yzhy5DCffvIJA2JiiI6O7rDznmpJ4m+pOTzyQxpKuYxnpvZlVNg5bs8M1o2WNi1Hsf0dcPPHdOHTWEImnHTY4dIsXkp8kX+O/cMg70Hcq1tEH7c+535eoUsSoUAQ2oPFjKzoYEOjn7pfTyr9N3rnr41E6hEMZ2gy0zgE6A3WyYH5VfkoZAqiPKLqmwX19xogQkAL1NTU8MH776FSqZgwcSJyees2OTobp+tTYCit5sFvU/l7Xx5XDOzF/RPDcFSf+z4Qsry91lsKB//CHHIhpvFPg0fwScdtOraJF/XLOVZ+jFlhs7m1363nxQRTwUqEAkForfJcZIZ05LmNlvzl7kZmqgTqSv9RWLRRDff+vcJB3bILtiRJHCw5UD8fINGQWBsClPT1jEKn1VlvB3jF4KhybM9X2i398fvvJCXt4NLJkzv8Z8OZmhdJksTqhCM89+sevF3teX56Pwb4u5/7CSUJ+Z6fUK57GMpyMA9egHnY3Sf9W6wx1/DZnk95P/U9HBQO3B5zB5cGTUYu67jAJNiGCAWC0FKNS//1DX/SkVXkAnWl//DaCX+R9SHgdKX/5kiSxIGSA/VVAH2OnoLqgiYhQKeNp7/XABECWunw4cN8+snH6OLjiYqK6vDzt7Sj4YG8chavSSb1eCm3jQxi3qhgVKfZU+GMjBUotryGYsvr4OiJadzjWCKnnvTv1FBh4LWkV/n10C/09ezHIt1i+nn2O/fzCp2eCAWCcKLGpX9Dau16/zRkBQeQISEhQ+oR2FDyP4vSf3MahwB97byAhhDQt8nqAAcxM7zNmEwmPnj/PRQKRYffNqhzNm2OTWYLb/19gJV/ZRLV04XnZ0QT7NXK20NFh2o3WvoZS8AFmCY8a12+eoIdhkSW659nb9FepgRfxu0D7sDD3qN15xY6JREKhPNbeW7tWv9G9/1z9zSU/h29Gt7115X/z6L03xxJksgsyazvE5Bo0FNYXYhSrqSvRz903jriNDoRAtrZX39t5L+tW7nk0ktxd3e3yRjOZe+DXUeKWbQmmeySKhZPDOPKgf6tXikgy1xv3Wip8CDm+Jswj1gM9k23sTZZTKzN+IY3d63EIlm4Nfo2ZobORCkXXRG7ExEKhPODscJ6n79Jw5/mSv99G5b+aSLBWdvqU1skC5nFmfXzAZoLAdbbAf1FCOggeXm5fPD++/Tr148BMTE2G8c5bYgEVNSYWP77Pj7bdpgRIZ48M60fWhe7Mz/wdMw1KLa9hWLTi6BywDTmYSz955y00VJRdSErd61k7f5vCHILZrFuMfHeA1t3bqHTEKFA6F4sZusM//oJf6co/dcu9at/998j6JxK/80OoVEIqKsEFFUXoZQr6efZr36JYH+v/qJRjA1IksSnn35CSXExl06e3KrNjlrrXENBnb/25vLgt6kYzRJPTIliYl/v1g+q9DjK9Y+hSF2DxTfOekuhmY2WdhfsZrn+eXbmJTHe/0Lujr0bH6eerT+/YFOdIhT8/fffLF++HL1ez/Hjx1m7di1Tp06t/7okSTz66KO8++67FBUVccEFF/Dmm28SGhra4nOIUNANlRmQ5dbO+q/r+Je3t1HpX1P7rj+yzUr/zbGGgAwSDPra/QMSKa6xhoBoz+j6EBDtFS1CQCeQnJzMTz/+wPgLL6RnT9texFobCgAKymt47Ic0fkszMDWmJw9dHIGLfetL+rKsLSh/X4LMkIZlwJWYRj8MTk13V5QkiV8O/sxrO1+ltKaM66Ou55rIa7FTtLJqIdhMS0PBuS+ObYHy8nIGDBjADTfcwPTp00/6+vPPP89rr73GRx99RFBQEA8//DATJ04kLS0Ne3v79hya0BmcWPqvqwBU5AEgKR2sG/to+2Lpd3mblv6bY5EsZBTtr10iqK8PASq5in6e/ZgZOlOEgE6qqqqKDevXERgYaPNA0FY8nNS8OnsA3+08zpM/7WbbgUKem96PQUGtmwgoBQzFeMM65Ikfofz7GdR7fsQ0cimWuLn1bbRlMhkXB13CqF6jeS/lXd5LfZcfMr9nYdy9jPIbJboidmMddvtAJpM1qRRIkoSvry/33nsv9913HwDFxcV4e3uzatUqrrjiimafp7q6murq6vrPS0pK8Pf3F5WCzqyu9F/X47++4c/BhtK/R1Bti9/I+vv/kntgm5X+mx1WkxBQVwkoRiVXWSsBtXMC+nlGY68UIbUz+/OPP0hK2sFlU6fi6Gj75ZxtUSlo7GhRJfevSSEhq5Drh/Vm4bhQ1Mo2WFVRnofyr6eRJ32CpInENPFZpIBhJx12sOQgLyW+wObjmxnqM5R7dYsIdA1s/fmFDtMpKgWnc+DAAbKzsxk/fnz9n7m5uTF48GC2bNlyylCwbNkyHn/88Y4apnA2JAnKa0v/hkYNf5or/YdObFjy5xUOHbAu3yJZ2F+0v35OwI7GIcCrP7PCZteGgH4iBHQheXl5JCbqGRAT0ykCQXvwc3fgo+vjWbX5EC+v28em/fksnxFNhE8rOxI6eWG6+GVksdei/H0p6k+mYI6ajmnc4+DSUHEJdA3k1VGv8/fRv3kp8QVm/zyLK8PncGO/m8U+G92MzUJBdnY2AN7eTSfQeHt713+tOUuXLuWee+6p/7yuUiB0sJry2oY/aQ0hwJCGrDIfqCv9RyB598PSb2ajbX5b0ev9LFkkC/uK9tUvEdyRm0hJTQlquZp+XtEiBHQT69b9iZOTE5GRJ6/D704Uchk3Dg/kghBPFq1J5vK3t7JwXChzh/VGIW9dOV/qGYvx2p+RJ3+Bcv0TqN8agvmCezAPuq1+oyWZTMaoXqMY0nMIH6d/zIdpH/DzwV+4M+ZOLg68RNxS6CZsFgrOlZ2dHXZ2YrJLh7GYkRVmNur1b132R13pXyZH6mEt/Vvib2x499/Opf9mh3qaEBDt1Z8rwuag87aGADFhqnvIzMzkQGYmo0aNsulqg44U4ePC17cM5pX1+1n+x1427s3l2en98HNv5TwXmRxL/znUhF2CYtPzKP56BvmuTzFf+AyWPg0VXTuFHTf1u4lLgy7hlR2v8OjWR1izfw2LdYuJ8Ojewex8YLNQ4OPjA0BOTk6TiUE5OTnE2HB98XmrSek/tSEE5O1FZqqyHuKkRdJEYrZB6b85ZovZGgIa3Q4oNZailqvp79WfOeFXotPq6CtCQLdksVhYv34dWq0W/4AAWw+nQ9mpFNw/MZzRYRqWfJPClBVbePiSCC4b0LP179jtXTGPfwrLgKtR/r4U1RdXYA6diGn8U9AjqP4wH6eePDv8ObZnz2B54nKu+e1qpvaZxoIBC3C369HKVyjYis1CQVBQED4+Pqxbt64+BJSUlPDff/8xb948Ww3r/FBTXjvrP61hwx9DekPpX+VobfjjHY2l36zaXv9RJy1b6mjWELDXukTQkECSYQelxlLsFHb09+rPlRFXiRBwHklJSSEvN5dJF1103pauBwd58P2CoTz1827u/yaF9bsNPDY5Cg8ndaufW9JEYLzyG+S7f0C57mHU7wzHPOR2zMPuavJGYKDPID6b9Dlf7/uKt5LfZN3hP7kteh7TQ2aglHe5YvR5r13/xsrKyti/f3/95wcOHCApKQkPDw8CAgK4++67eeqppwgNDa1fkujr69ukl4HQChaTtblP3Xr/3PTTlP5vsu74p4mEHoEndTuzBbPFzN6iPehz9OhzrZWAMmNZfQi4KuLq+hCgVrT+h6DQdZhMJv75+y969+6NRtNx81Q6Ixd7Fc9Nj2ZMuIZHv09nyorNPDOtHyND2yDEy2RYIqdQ02ecdaOlrW+gSF6NadyTWCIm12+0pJQruSJ8DhN6T2TFzjdYrn+etRnfsEi3mDitrvXjEDpMuy5J3LhxI2PGjDnpz6+77jpWrVpV37zonXfeoaioiOHDh7Ny5UrCwsJafA7RvIja0n9O03f9db3+zdblm5KTtv6i36Thj6rzrLdvEgIMCezI3VEbAuwZ4NWfOG088d46ojz6ihBwntu2bRsb1q9jymWXdcr/9229JLGlDKXVPLA2hX/253PlIH8WTQjFUd2G7/0KD6D88yEU+37D0nuEtSuiJvykw1LzU3le/xyp+SlM7D2JO2PuwtuxDboyCuesU3Q07AjnXSioKbNe7Osm/dUt/assAGpL/5oI60VfE4Wldrc/W5f+m2OymNhbuLd2TkACO3KTKK8PAQOI0+pECBBOUlNTw5srV+Dn58fQYSevqe8MbBUKwNoD5vPtR3jutz34uNqzfEY0/Xu5nfmBZ0Ge8SeKPx5AVpRl3Whp+GKwb/rz1yJZ+OnAj7ye9BqV5kpuiLqRqyKuFv+XbUSEgq7upNJ/GnJDOrKigwC1pf/ghnf9naz03xyTxcSewj21jYL0J4WAumZBfT36olKIHdqE5m3ZvJl//vmbqVOn4uTcOdfI2zIU1MnMK2fxmmTSjpcyf1Qwt40MQqlow58NpmrrRkv/vgRqJ0xjHsESPeuknz9lNaW8k/IOX+xdja+TL/fG3cdwvxFtNw6hRUQo6CrqS/9pJzf8abb0X9vtzzOsU5X+m9M4BOgNepIMOyg3lWOvsGeAJgad1hoCojyiRAgQWqS6upo3V64goHdvBg8ebOvhnFJnCAUARrOFt/7K5M2/D9C3pwvPz4gmyKtt9wih5BjK9Y+iSFuLxS/eekuhZ8xJh2UWZ/KC/nm25WxjuO9w7o27D3+X82vViC2JUNAZ1ZRZZ/03bvhzqtK/tm/9/X8cPW088JYxWUzsLtxd3ydgZ25SfQiI0cSg08YTp9WJECCcs/oqwbRpODm18cWtDXWWUFBn5+EiFn+TQk5JFfdPCueK+F5tvmJDduhflL8vRZabjiXmGkyjHzzpZ5ckSaw/sp6XE18kvyqfqyOu5vqoG3G00bLm84kIBbZkMSEryGyyxW/zpf/ad/219/5x791pS//NMVlM7C5Ir99AKCl3BxWmiiYhQOdtDQFKuQgBQuvUzSXw9/dn8JAhth7OaXW2UABQUWPiud/2snr7EUaGevH01L5oXdp46a7FhFz/Icp/ngVktRstXVe/0VKdKlMlH6Wv4qO0j3C3c+fOmLuZ2Hviebu0tCOIUNARJAnKspu2+c1Nb1r6d/ZuMuGvq5T+m2OyGNldsJsEg55EQwJJuUlUmCpwUDoQ4xWDztu6lXCkR6QIAUKb2759G+vXrWPqtGk4d9K5BHU6Yyio89feXB74NhWTReLJKVFMiGqHVQHluSg3PoV852dI2r6YJixDChh60mFHy47yyo6X2HBkA3GaOBbpFhPao+Wrz4SWE6GgrTUq/cvrdvvLTUdWWQiApHKylv7r3vl3sdJ/c0wWI+l1lYCcBJLykqg0VeKgdCBWE0ucNh6dNk6EAKHdmc1m3npzJRqNhguGD7f1cM6oM4cCgILyGh75Po0/0g1Mi/HlwYvDcbFv+//DsmOJKH+7H/nxHZj7zsA09rEmGy3V2Xp8Ky8kPk9WaRYzQi7ntuh5uNm17YqJ850IBa0gy9uLrFGrX3luOrKiQ0Bt6d+jT8O7fm0UFk0UuAd0qdJ/c0wWI2kFaehz9CQa9PUhwFHpSIwmxtonQKsjwiNChAChQyUnJ/PTjz8wecoU3N3dbT2cM+rsoQCs9/e/TTrGkz/vxt1BxXPT+zEw0KMdTmRBvutzlBueBFMl5gvutW60dMLSRKPZyBd7V/Nuyjso5SoWDFjAZcFTUXTwHirdlQgFraB+IwZZyZGmpX9t39qGP2HQTXbUaxwC9IYEdubtbBQCYq2rA7zjiegRIdqVCjYjSRLvv/8edmo7xo4ba+vhtEhXCAV1jhRWcv83yeizirhhWCB3jwtBrWyHNzhVxSj+fg6F/n2kHoHWWwrBJ/995lXm8cbO1/nxwA9E9Ihgke5+BmgGtP14zjMiFLSCzJCG5OwDju2Qmm3IaK4NAYaG1QFV5iqclE7EaK0hIE6rEyFA6FQOHDjAF6s/58IJE+o3UuvsulIoADBbJD749yCvrt9PsJcTyy+PJtzbpV3OJTOkofx9KfKsfzGHXWTdaMm990nH7crbyXL986QXpHNx4CXcGXMnXg7nd0vr1hChQMBoNpJakGoNATl6duXtPCkE6LTxhPcIFyFA6LS+/PILCgsKuOTSS7vM7PSuFgrq7M4uZdHXyRzIL+ee8aHMHdobubwdvueShDz9W5TrHoWKfMxD78A89M6Tdlw1W8x8n/kdK3a9gdFs5KZ+N3NF2ByxpPkciFBwHjoxBOzM20m1uQonlTOxdUsEtfGE9QgTIUDoEgoKCnjn7bcYNmwYfUJCbD2cFuuqoQCg2mjmlXX7+XDLIQb27sFz0/vh695Oq6VqylFsfgXFfyvASYtp/FNYwi+p32ipTklNCW8nv8VX+77E3zmA+3T3MbRn52xx3VmJUHAeqDHXkFaQSkKOdYngzrxd9SEgThNLXG0lQIQAoav644/fSU1JYfqMGSgUXWfCWVcOBXW2ZhawZG0KpVUmHrk0gin9e7ZfpaYgE+WfD6LY/weWoNGYLlyG5BV60mH7i/bxvP55Eg16RvmNZmHcPfRy7tU+Y7IRSZKoqq5CIVegUqna7HsuQkE3VGOuITU/pbZZUAK78pKpNlfhrHImVhNb3ywozD1czNgVuryamhpWvPE6ISEhxOm61va73SEUAJRUGnnq5918t/M4k/p689jkSHo4tt+GRvJ9v6H88yEoPox54C2Yhy8Cu6ZzGyRJ4o+s33k16RUKqwq5NvI65kbNxV7Z9Xq/NKeyqpJXX34ZgHvuvQ+1um2+3y0NBeLtYydWY64hJT+ZRENibQjYRbW52hoCtHHM6z+feK2OUPcwEQKEbic9PY3q6mrCwk/emlfoGK4OKp6fEc2YcA2P/ZDO5BVbWDa1LyNC22fXVUvoRGqCRqH4700Um19GkfI1prGPYek3s/6WgkwmY0LviYzwG8kHqe/zUfoqfjzwA3fH3sM4/3FdZt5JcyyShef1z2JHG3eaPAsiFHQidSGgrllQcn4y1eZqXFQutSFggQgBwnkjaccO/Pz8On33wvPBRf18iAtw54G1qdz0cSJXDfJn0YQwHNTt8HNIaY/5goWYo2ehXPcIqh/mY9mxyrrRkk//+sMclA4sGHA7U4Kn8FLiiyz5dzEDvQdyn24xfdz6tP242pkkSSzXP88vB39hKlNtNg4RCmyo2lxNSn5K/QZCKSeEgPn9F6DTxhPqHipCgHBeycnJ4fjx44waPdrWQxFqebva8961cXy27TDP/baXzZn5LJ8RTbRfO3UedPXDNO19zLFzUf6+FNUH47DEXodp1ANNlov7uwTw8qhX2XRsEy8lvsCVv1zBrNBZ3BJ9Gy7q9llW2dYkSeL1na/x1b4veSD+QfZm7bHZWEQo6EDV5mqS85JJrN1AKDlvFzWWGlzVrsRq4ljQ/3Z03vGEuIWIECCc13bt3ImDgwO9enWvSWRdnUwm46rBAQwJ9mTxmmSueHcb80cHc+uIIJSK9unoKgWOwHjjBhT6D1D88yzq3d9hGvUAlphrodHPyeG+wxnkPYjP9nzG+6nv8uuhX7l9wB1MDp6CvJN3m30/9X3+l/4R98Tey6VBk3kJEQq6pboQUNcsKCUvuT4ExGl13BFzJzptPCHuIZ3+H60gdBSTyURqagrBwcHI5eL/RWfUR+PE6psHsXJjJm9syOCvvXk8P6MfgZ7ttJ21QoV50K2Yo6ah3PgUql8XYdnxP+stBf/B9YepFWrmRs3l4sCLeS3pVZ7c9gRr9q9hcfz99PPs1z5ja6VPd3/CW8krmdd/PldGXEVNTY1NxyNCQRuqMlXVzgmw9glIzk/GaDHipnYjVhsnQoAgtEBGxn6qqqq6VF+C85FKIeeucSGMCvNi8ZoUpq7cwpJJ4cyO79V+k/2ctZgufQ1z7HUof78f9ceXYO43C9PYR8C5odul1lHLU8OeZkbIDJbrn2fu79cyOWgKtw+4A0+HzrNJ3dr93/DyjpeYG3k9N0TdaOvhACIUtEqVqYrkvF3oDdYNhBqHgDhtHHfF3I1Oq6OPCAGC0GIpySl4enp2iY2PBIjxd2ftvCE899teHv0hnfW7c3l6al80Lu03g17y02Gc+zvynZ+i3PAU6r0/Yx5+H+aBtzTZaClWG8fHEz9lbcY3vLlrJeuPrOfWfrcyK2yWzTd1+/XgLzyz/Wlmhc5mwYDbO82qCREKzkKVqbL2doC1T0BKfkptCHC3hoDYhcRrdQS79REhQBDOQVVlJZmZGcTGxtp6KMJZcLJT8sSUKMaGa3jwu1QuXbGZp6ZEcWGUd/udVCbHEnMNNeGTUfzzHIoNTyDf+al1o6Wg0fWHKeQKLg+dyfiAC3lz10pe3vES32Z+y6K4RQz0GdR+4zuNjUc28OjWR7g0aDL36RZ1mkAAIhScVpWpkl21lYC6EGCymEQIEIR2smfPHiwWC4FBQbYeinAORodr+GHBMB7+Po3bV+9keqwvD14UgbN9O15qHNwxT1iGJeZqlL8tQf355ZjDL8U07gnrlva13O3cWTrwAab1mc5y/fPM23Ab4/zHszB2IT5OPdtvfCfYcnwzS/9dwlj/cTw06OFOd+0QoeAEkiTxTsrbbM/ZVh8C3O3cidPqWBh7DzptPMFuwZ3uL1IQuoO0tFS8fXxwcOge3enORx5Oat64YgDf7DjGUz/v5r8DBTw/PZr4wB7tel5J2xfj1d8jT/sG5bpHUb8zDPPQOzEPuQNUDf+eIjwieG/8+/x66BdeTXqFGT/NYG7UXK6JuBZ7pX27jjHRoOe+f+5jiM8Qnhz6ZKdcZSZCwQlkMhmp+Sl42ntxT+y96LzjCXINEiFAENpZeXk5WVlZDB48+MwHC52aTCZjRpwfg4J6cP+aFK7+cDs3Dw/ijjF9UCvb8WepTIal7wxqQiei+Pdl68eu1daNlsIuatIV8aLAixnpN4r3U9/j/dT3+CHze+6Ju5dRfqPbpZyfkp/Cwr/upr9XNM8Of97mcxpORVzpmvHa6Dd4bvjzzAqbTR9xa0AQOsS+vXsB8A8IOMORQlfh38ORj28YyD3jQ/lw80FmvvMfe3NK2//EamfMYx7GePM/SJ6hqNZci2r1LGT5+5oc5qRy4s6Yu1h90ZcEugZy3z/3csfGBRwsOdCmw9lXuJc7N95OH/cQXhzxMnYK27UxPhNxtRMEoVPYs2c33t7e2Nu3bwlX6FgKuYxbRgTx1S2DMVsszHj7P1ZtPoTF0v578UmeIRhnr8Z4+cfICjNRvTsSxfrHoLppMAl0DeTVUa/z4oiXOVx2hNk/z+aVHS9TZixr9RgOlhxkwYb5+Dr58uqo13BUObb6OduTCAWCINhcVVUVhw4dElWCbiyypytrbh3CVYP8WfbrHuZ+lMCxosr2P7FMhiXsImpu+Rfz8HtRJLyP+u0hyFO+gkabBMtkMkb1GsWXF3/FLdG38NW+r5jx43R+PPAjFslyTqc+VnaMBRvm4W7nzhtjVnSJtssiFAiCYHOZmRlYLBb8/f1tPRShHdmpFCyZFM6quToO5VcwZeUWvt95HElq/6oBSnvMw++j5tbNWHoNRvX9PFQfT0aWk9x0jAo7bux7E2su+YY4bSyPbX2Em/68gd0F6Wd1utyKXOZtuA2VXMXKsW/ibte+Ey3biggFgiDY3L59+/Dw9MTJqZ3a5AqdytBgT35YMIzRYRoWrUlm4Ve7KKowdszJ3fwxTf+AmjnfQGUBqg/Gofx1MVQWNjnMx8mHZRc8x1tj36bcWME1v13N09ueoqi68BRP3KCwqpD5G27DZDHx5ti38HLQtNeraXMiFAiCYFMWi4XMjAx6+fnZeihCB3J1UPHC5dG8PLM//+7PZ/KKzfy7P7/Dzi8FjcR401+Yxz2OPOVL1G8NRp64CizmJsfFew/k00mfcW/cffyR9TvTfpjKF3tXY7KYmn3e0ppSFmyYT0lNCW+OfYueTr4d8GrajggFgiDY1NGjR6mursZXhILz0sXRPvx4+zBCtc7c8D89T/20m8oa85kf2BYUKsyD5lFz239YQiag+vU+VKsuRHZke5PDlHIlV4TP4ZtLv2VcwHhe0C/nmt+uQm/QNzmuwljBnRtvJ6cimxVjVhLg0vXmyNg8FDz22GPIZLImHxEREbYeliAIHeRAZgb29vZ4eXnZeiiCjXi72vPeNXE8dHEEX+qPMP2traQcK+m4ATh7Y5r8BjXX/gKA+n8XofzhdijLaXKYh70HDw16mFUT/oedwo5b193MA/8uJacihypTFff8fTeZxZm8PnoFIe6hHTf+NtQpmhf17duXP//8s/5zpbJTDEsQhA6QmZmJj49Pp+r/LnQ8uVzGNUMCGNbHg0VfJzP7nf9YMLoPt4wIRKnomPevUq+BGOf+gXznJyg3Po1670+YRyzGrLsJFA3Nhvp69uWDC1fx04EfeT3pNWb8OA1vR2+Ol2ezcuxKojyjzvrclaZKdubuRH88oS1f0lnrFFdfpVKJj4/PmQ8UBKFbqayoIDs7m6HDhtl6KEIn0UfjzOqbB7Pyr0xe37Cfv/fl8tz0aHp7dtD6frkCS+x11ERMQfnXMhTrHkWe9AmmC5chBY1sOEwmZ3LwFEb4juT6P+ZyqPQQXvYaSmta1pypPgQYEtAb9KTmp2KWTHiqPBnDmPZ6dWdk89sHYJ157OvrS3BwMFdddRVZWVmnPLa6upqSkpImH4IgdE2Hav+v9xRvCoRG1Eo5d48L4bMbB5FfXsPUN7fwZcKRjlm6WMehB6ZJz2O8YR3Y90D9+XSU39wAxUfqD7FIFl5Jeplj5UdZFHc/wW7BLPz7bu76606ySptexyqMFWw9voUVO1/nhj/mMvrrUdy+cT7fZXyL1kHLfbr7+PLir/l+8o8d9xqbYfNKweDBg1m1ahXh4eEcP36cxx9/nBEjRpCSkoKLy8mNHpYtW8bjjz9ug5EKgtDWsg4dxMXFBSdnZ1sPReiEYgPc+XbeUJ79dQ8Pf5/G+t0GnpraFy/njmsTLHlHY7zmB+SpX6Nc/xjqt4diHnY3psHzWb7zFX468CNPDXuGib0nMitsFhuOrOflHS8x++eZjO41Gq2DN8n5u0jNT8MsmfCw80DnrePiwEvQaXUEugY1uXVmPJbaYa+tOTKpQ6PXmRUVFdG7d29eeuklbrzxxpO+Xl1dTXV1df3nJSUl+Pv7Y8jNw9XVtSOHKghCK7337ju4u7t3u9sHG9avx87OnpmzZtl6KN3G+t0GHvouDYsk8dRlfRkfqe34QVSXovj3ReTb3uZl756ssoeHBj7M1JBpVBgr2JmXhD5Hz/acbaQVpCEhIUdOlGdfLgm8lIE+8fR2CTzl/BnZ0QSkL67h2eprALjn3vtQq9VtMvSSkhK0Gi+Ki4tPe620eaXgRO7u7oSFhbF///5mv25nZ4edXefdTEIQhJapqKggLy+PcLHaSGiBsRFafujlxkPfp7Hg8yQuj/Nj6UXhONt14GXMzgXz2Md429meVRlfMrOklKPb32Du3tWkl2Rilsx42nuh0+qYHHwZAS7+fLXvKzYcWY9KrmSAZsCpA0HmelRr5lKtiYHjHfeSTtTpQkFZWRkZGRlcc801th6KIAjt6MgR671ZrdYG7/iELsnT2Y6Vc2JYs+MYT/+8m60HCnh+ej90vdu/hXC5sZyk3CT+l/4RekMCMmR85eqCl6WQ+GPHuKznEGKHLKK3R2STC/8gn8FsPb6VFxKf5+rfrmRGyAxui56Pm51b/THy1G9Q/rAAS/AYjJe8Ca+taPfXcyo2DwX33XcfkydPpnfv3hw7doxHH30UhULBnDlzbD00QRDa0dGjR3B0dBStjYWzIpPJuDzOj8GBPVj8TQpXf7Cdm4cHcfuYPqiVbTd3vsxY1mR1wO6CdMyStalSsGswV4TPQaeNJ8Beg/K/FSi2vAqZCZjGPoYlajo0CgZDeg5h9UVfsHrvF7yb/Da/H/qD+QPmMzV4GqrEVSh/X4Ilehami18Bs23v6Ns8FBw5coQ5c+aQn5+PRqNh+PDhbN26FY2m6/SKFgTh7B09cgSNRiP6EwjnxN/DkU9uGMh7mw7w2voM/t6Xx/LLownVntuk1TJjGUm5SehzakNAYToWyYKXgxfx2oH0cQ3m+wPfMyt0Not0i5v8uzWPWIQ5ejbKdY+g+u5WLImrME18Fknbt/4YpVzF1RFXM6n3JN7Y+TrLtj/D2qSVPHBoN9GD52Ee+xjI5GCuaeV3pnVsHgpWr15t6yEIgtDBzGYz2dnZDIiJsfVQhC5MIZdx68hgRoR6cd/XyUx/ayv3jg/l2iEByOWnD5tlNaXWEFBXCSjcjUWyoHHQoNPGM63PNHTe8fg7+7PxyAaW/Hs/U4Iv4z7douaDrHsAphmrMB/YiPL3pajeH4Ml7gZMI5eAg3v9YV4OXjw26BFm5h3nubx/udbXh4scLNxVld8pNk6yeSgQBOH8k5eXh8lkwsvT09ZDEbqBqJ6ufHPbEF78Yx/Lft3Dhj25PDu9Hz3d7OuPKaspZUfuDhIN+iYhQOugRafVMT1kBjqtjl7O/k0u+puP/cvSzUsY6z+OBwc+hFx2+lsUUtBojDf9hWL7Oyg2LUedthbTmIewDLiqvhKg/H4+sbu/53+TXuRbV2fe2Pk6f/24kZv63czlgZe317epRUQoEAShwx0/fgyZTEYPDw9bD0XoJuxVCh68OIIx4RqWrE1h8pvrmDPChMIxE70hgT2FexpCgHf8KUNAY3qDnkWb7mNoz2E8OfRJFHJFywajUGMecjvmfpejXP8Eqp8XYtnxEaaxj6Lc/CqyrC2Ypn8I4ZcwFRjrP463k99ixc43+GHfD+iIa7Pvy9kSoUAQhA6Xk52Nm5sbKpXqzAcLQguU1pSyIzcRfbGeXn0T2Fu0h8+zJNR4MLLXIGaEzCReq8PPuVeL5rGk5Kew8K+7GOA1gGcveA6l/Bz+rTr7YJqyEnPsdSh/vQ/1p9OQ5EqMU99DCr+k/jBXtSuLdIuZ1mcay7ctP/vztCERCgRB6HDZ2dmiSiC0SklNCTsMdbcDrJUACQlvRx90Wh2zw2ZRXBjAa78XsvmokinT+tHLpWW3q/YW7uWODQsIdQ/jxZEvY6doXW8cya0XmI1IKieQy1H9dAfmkqOY466HqgJkhnTkuWlEGNJ505DB8wS26nytIUKBIAgdymKxkJuby4ABA2w9FKELqQsBdRMD9zYKAfHaeGaFzUanjcfXybdJJWBCaBVL1qZw/Ud6rhkSwH0XhmKvOvVtgIMlB1iwYR5+zn68OupVHJQOrRq3LG8fqtWXAzJMl74GpcdR7PwUxZ8Polj3CLLaZY6S0gFJEw7afmBo1SlbRYQCQRA6VFFhISaTiR492r/hjNB1FVcXsyO3NgTk6NlXtBcJCR9HH3TaeGafIgScyMfNng+u1fHJtixe+H0fmzPyWT4jmr6+J7f6PVp2lPnr59HD3oM3xqzAWX3y/jtnZDEhKziALDcNecY65ClfASCzGFGtvRFJJkfqEYSl9wjkhZlQchRzn/GYJj5vXcFQUwMpL5z9eduICAWCIHSo3LxcANxFKBAasYaARPQGfZMQ0NOpJzptPHNqmwX5Ovue9XPL5TKuHdKbYcGeLFqTzKx3/uOOMX24aXggSoV1NYGhwsD89bdhp7Bj5ZiVuNud4d+nJEF5DjJDGnJDGrLcdGSGNGR5e5GZrfvzSMisrZH7zkDqGYukjUTyCgeVY/1zyFO+sm609O5wzMPuhtibz/r1tSURCgRB6FB5uXnY2dlhb29/5oOFbqs+BNQ2C9pXtA8JCV8nX3RaHVeGX0mcVndOIeBUQrTOfHHzYFZszODV9fvZuDeP52f0w8WxmvkbbsMsmXl73Dsn9wuoKUOWu9t60c9NbwgBlQUASCpHJK9wJJ/+WKJnQ2Uhis2vYgkejWn6Bw0h4EQyGZboWdSEXYRi03IU/zyPeueXwJQ2e81nS4QCQRA6VF5+Hm5ubqKT4XmmqLqIHYbaSoAhgX1F+wDwc/IjThvHleFXtXkIaI5aKWfh+FBGhnpx/zcpXPbWOnpGrMIsK+XdsW/Ts7IU2aH/rBf93DTkhnRkRQcB6kv/krYvlvibkbRRWDSR0CPQ2oMAkOs/QPnvS1j6XY7pktdA0YJVC3YumMc9gWXA1Ui/PgjF7ff6z0SEAkEQOlR+Xr7Y5vw8UFRdSGJtCEisrQRAXQjQcVXENei0cfR0at8Q0CxJQudRxZcTi7kx8WXyKkpYnicn5K0RDaV/Jy2SNgpz2EVI2igkTV3p/xQTDyUJxaYXUP7zHKaBt2Ie/2R9UGjxsLzCMM78BF56sbWv8JyJUCAIQoeRJInCwgJ6+fey9VCENtY4BOhzEthfvB8AP+de6LQ6ro64Bp1Wh49Tz44d2Eml/1RkuelUVxXxoLeGQjs1r5e7cqzShxcYycjhI4nRDQMnr5afQ7Kg/P0BFPr3MI160Do34FwrYTauoIlQIAhChykvL8doNOLicg6zuoVOpbCqkMRGcwIyakNAL+de6LTxXBN5bceGAIsJWUGm9eJvSLOW/nPTkRUdAupK/8FI2iiqdTeysHQHyeVZvDF6BTHaOHqVVfP9d2ms+DOXmYU5LJnkjrNdCy6R5hqUP9yOPP1bjBe9iCX2unZ+oe1LhAJBEDpMUVEhgAgFXVBhVWH9fIBEQ2J9CPB39idOq+O6yOuI0+rwcfJp34FIEpRlN7zrN6Rb7/83nvVfV/oPvwRJE2Wd9e8ZBioHTBYTD2xeyrayg7wy6jVitNaWwl7Odrx5ZQxfJx7lmV/2sDWzgOdnRBMX4H7qsdSUo/rmemSHNmGa9j6WiMnt+9o7gAgFgiB0mKIi6wwqZ+dz295W6DgFVQW1twOslYDM4gzAGgJ03vHMjZpLnFaHt6N3+w2iutRa+q+b8JdbN+vfGi4llROSJgKp5wAs/edg0UQiaaPAsfnOhRbJwhP/Pc5fR/5i+YjlDPYZ3OTrMpmMmbpeDA7yYPGaZK56fxu3jAhiweg+qJUnzA+oLET15RxkuekYZ69GChzZLt+CjiZCgSAIHaa4uAh7e3ux50EnZA0B+vpmQZklmQAEuASg0+q4Pur69gsBFhOy/AzrRb/24n9S6d+jD5ImCsvAkdZftZHg3rvFk/kkSeL5hOf45eDPPD3sGUb6jTrlsQEejnxyw0De3XSQNzZk8Pe+PJbPiCZEWxtmS46hWj0TWXkexqu+ReoZ2+pvQWchQoEgCB2mpLgEJycnWw9DAPIr80nMtTYKSjQ0DgG90WnjuKHvjcRpdWgdtW130rrSvyENeV0AMKQhy9+LzFxjPcTZG0nTuPQfheQVBspz72shSRKvJb3C1/u/4uFBjzCh98QzPkapkDNvVDAjQr1Y9HUy09/ayn0TQrmmTw12X8wEwHjtT0ieIec8rs5IhAJBEDpMSUkxjo6naOQitKu8yjwSDfr6WwIHSg4AdSFAxw19b0Kn1aFx1JzhmVroxNJ/7ax/WVUR0Kj07xuDecBV1vv+mihwbPuNst5LfZePd3/MfbrFXNZn6lk9tp+vK2vnDeGFP/bx/S+/MNthOXI3b6Sr14BLB6+k6AAiFAiC0GFKS8vw9BS7I3aEuhBQNznwYMlBAHq7BKLz1nFTv1vQaeNO7t53tpqU/tMaQkBxFtCo9K+NwhI0GkkbiUUTBe4BZ72O/1x8kv4xbye/xYL+t3NF2BXn9Bz2KgUPRxqQpzxDurkX8wqXcO9BGRdHt/FgOwERCgRB6DBlZaWiR0E7yavMRV9bBUg06OtDQKBrIHFaHTf3u7V1IaC29G+d8V/X6z/9hNK/j3XWf8TkRg1/Wlf6b42v933NK0kvc33UDVzf94Zzfh55+vcov78NqfcItBe9Q//fDrHwq12s22PgkUsicXPoPnNkRCgQBKFDmM1mqqqqcBB7HrSJ3Ipca7fA2nkBh0oPAtYQoNPGc0u/W4nT6vByOIsmPHWqS+s3+JHXbfTTuPSvdrJe8H1j2730f65+OvAjzyUs44qwOczvv+Ccn0e+4yOUv9yHJWoapslv4K5Q8/Ks/oyL0PDYj+lMObiZZ6f1Y2if5lc8dDUiFAiC0CEqKioAsHdo3f7056u6EFC3RDCr1DozP8g1iIHeA7mt/23EauLOLgSYjcgKMhou+ieV/hW1pf/I2tJ/FBZtFLj5d0jp/1ytO7yOx/97jCnBl3FP3L3nts+GJKHY/DLKv57BFH8z5gufbvKaL+3fE13vHiz5JoW5H+m5bmgA94wPxV6laMNX0vFEKBAEoUPUhwJRKWgRQ4WhYYmgQU9WqfVCHewazGCfQczrP484jQ5Phxa8Q5UkKD1eO+P/dKX/vpgjpzRq+BNqs9L/ufr32L88uHkp4/0v5IGBDyI/l/AiWVD8+TDK7W9jGrkE8wX3Ntt+uKebPR9ep+N/W7N48c99/Ls/n+WXRxPVs+vu7SFCgSAIHaKy0hoK7OzsbDySzslQYaifD3ByCBjM/P4LiNXEnTkENC791078k+XuPqH0H4XkG4c55mrrbQBtFDj0aOdX2P70Bj2LN93HsJ4X8MTQJ1DIz+Fdu9mI8qc7kad8jXHi81h0p5+LIJfLmDusNxeEeLLo62RmvfMfd4zpw03Dg1DIu95OoCIUCILQIaoqqwARCurkVOQ0rA7ISeBw2WEAgt36MNhnCPP7306cNg4P+1Pcp29S+k9ruP9fbH0eSaZA8gxB0kRiCR5r/bULlP7PVUpeMgv/uosBXjEsu+BZlPJzmPxnrED5zQ3ID/yFaeq7WKKmtvihoVpnvrxlMK9vyODldfvZuDeP56f3w9+jay3BFaFAEIQOUVVlDQXnazfD7PLsJksEj5QdAaCPWwhDew5jgVbXfAioK/3XrfOvK//n72so/bv0tDb8ibysaa9/5fkRwPYW7uWOjbcT6h7GiyNfwk5xDq+7shDVV1chy0m1ti0OOnXHw1NRK+Xce2Eoo8K8WLwmmSkrt/DgxRHMiPU9t3kNNiBCgSAIHaK6phq1Wt1lfji2VolUwk8HfrSuEDDo60NAiFsIw3pegM47njhNHD3sG5Xtq0qQHf7v5F7/VdY9IyS1c+2sfx3mmGtrZ/1HdovS/7k6WHKABRvm0culF6+OehUH5TlMZC09jmr1LGRlORivWovkG9eqMcX37sH384fxzK97ePDbVNbvNvDklCg8nTt/SBOhQBCEDlFdXd2tqgSSJFFsLCa7KpvsyuNkV2VzvDKb7KrjZCmyqLBUwFZrCLjAdzi62kqAu10Pa+k/fz+yjPVNl/6dWPrXRmEJHtfQ8MfNv9kJb+erI2VHmLf+NjzsPXl99Aqc1We/+6asIAPV55eDZMF4zU9IXqFtMjZneyXPTO3L2HAND32XyuQVW3jqsijGRrRh2+h2IEKBIAgdwlhjRKnsWqHALJnJq8rleFU22bUX/Oza3x+vyqbKXFl/rJvKDR+HnvS098GtyI2eCl9uv2wB7jWV1gl/R1KRJX5V2+t/HzKLEQDJxReLNgpz5NT69f7WWf+d/12lLeVU5LBg/TzsFfasGLMSdzv3s34OWfZOVKtnIzl4YJzzFbj6tfk4x0dqifF348FvU5n3WRKzdH4smRSOk13nvPx2zlEJgtDtGE1GlMrOt4a72lzd7Lv97MpsDNUGzJIZADlyNPYafOx9CHMNZ6R2FD3tvPCT7PCzgGtNKYqKPBQVueTm7cbdnIhm5aqmpX9tFFKvgZhjr2to+OPgbsNX3zUVVBUwf/08zJKZt8e9c04NmmQH/0H19TVIXmEYZ61u18ZLXs52vHVVLF/qj7Lsl91sPVDA89OjiQ1wb7dznisRCgRB6BAmowmFouNDgSRJlJpKG97p177LrwsBBTUF9ceq5Wp62vfEx07DEJcI/Fxj8LMo8DeZ8asux66qAGVhDorKVBSVuShqSpueS67C5OCF0aiiQq3FrLuioeGPay9R+m8DJTUl3L5hPmXGMt4b/x4+Tme/KZF8z08ov70ZKeACjDM+BLVzO4y0KZlMxuz4XgwJ8mDxmmSufH8bt44MZsHoYFSKzrMaRIQCQRA6hNlsQn4u68ZbwCJZyK/O53jtRT+7/ldrqb/cXF5/rKvCCV+lG35yB+LlGvxVnvjXVBNQWYZ3RT7KygTkpsqmz690wOzghclBg9nRi5oeIfW/r/9zBy8sdm4gk7Fh/Xrs7OyZecGsdnm956tyYzl3brydnAoD74x/F3+XgLN+DnnSxyh/uRdLxBRMU1aCQt0OIz213p6OfHrjQN755yArNmbw9748ls/oRx9N+weTlugUoWDFihUsX76c7OxsBgwYwOuvv86gQYNsPSxBENqQ2WxG3opmLjWWGnKqcuov+sfrL/rHyanKwSSZAJABWpkDfqjoa5aYZDQRUGWid0UhAdVVuEhSw5jUrtaLuqMGs5MfFV4DGj538MLsoMHk6IWkcmrtyxdaqcpUycK/7+ZAyUHeGvsWfdz6nN0TSBKKLa+h3Pgk5rgbME1YBu0UUs9EqZAzf3QwI0I9WbQmmWlvbmXRhDCuGuRvk/E0GZutB/DFF19wzz338NZbbzF48GBeeeUVJk6cyJ49e9BqO/csTUEQWs5isSCXn75MWmYsa3i3X3GEnPJDZFceI7s6jzxzGXWXc7UEfhYZ/kYTo2sq8Dca8Tea8DeZ8DWaUdj3sF7UHa3v5k2e1gt8haMXJbUXe7ODJ1IXa+F7vqox17B40yLS8lN5Y8xKIj2izu4JJAuKdY+i3PYmphGLMQ9f1Clu5UT7ubH2tqG8+Mc+nvp5Nxv35PL4pWE2HZPNQ8FLL73EzTffzPXXXw/AW2+9xU8//cQHH3zAkiVLbDw6QRDaikWSQAYFZUcxlOwhp+wA2RVHOF5l4LipiKOWCkox1x/vajbjbzLhbzQxxGSil8mCr9yJnqoeeDhokRw0mHp41ZbwNZgcrKX8ow6eILf5jzahjZgsJh7c/AAJOQm8MupVYjQxZ/cEZiPKn+9GnvwlxgnPYYm/sV3Gea4c1AoeuiSC0eFeLF2bysx3tnKZDcdj0/85NTU16PV6li5dWv9ncrmc8ePHs2XLlmYfU11dTXV1df3nJSUlbT6uyhozmXnlZz5QEITmSRKqiuNY8naSkbeZtLI9JNXkcFAJ1YkN79C8ay/2oZKSETiiVXjhpfTCQ+2LyqUn1XaeVNl5UW3niVHtRpFMTtEpzwlUABVGwNjuL/F0jlcpUVsUpB5r+59P5xOLZOHd3cv4L/cv7oh6EidL5Fl9T2WmSoI33oHLsX84OPIVCn0nQyf9O+nhqGb5jGheX7cH8qx/VmU0o+7YKQ+2DQV5eXmYzWa8vb2b/Lm3tze7d+9u9jHLli3j8ccfb9dxZeaVM/2tre16DkHoLlyoIEx2mAj5YXorDoJjFrkORexyULBbrcYik9HDBL5VrkRVeWGq0VBp9KG4xo8ci4b9OGCdCXA6RbUfXYV1idzKDPFz5NxJ2Pl8i8p9G1VHr+DJNDnQ8u+nK+W8p34Blewg1xnv5Z/fNWf1eFtRYia6tinjgbwKBjh17FbjXa7GtnTpUu655576z0tKSvD3b9vJGcFeTnxz25A2fU5B6PIsRuyLM3Eo3IND4W4cCvdgLNpDirmQ7Q52bLe35we1GkkGXjJP+jqFMNxzCCHaUWgd/Fj35x/UKIzEj4q39Stpdwnbt6NW2zH+wgttPZQuSZIkVme+ya9H/uPG8PsZOfqis3q8ssJA6B/XoarIYf/4z1moiWFhO421LZRWmXjr70y2ZBYwvI8nHLP+eZBXx2+mZNNQ4OXlhUKhICcnp8mf5+Tk4OPj0+xj7Ozs2n2XNQe1gr6+XXc/bEFoFUmC0mPWDXjq+u8brBvwlGAi0d6e7a6eJDg4sFujREKDt50HOp9BXO49CJ13PH5OfiftcbDbUaKysoYwTce+87GFo/Ym7OyU4ufIOXon+W1+PfIli3SLmR02++weXJCJ+tvZYDZivPYngjTh7TPINvLPvjyWfpuK0WThlVn9GRfmwUsv/gmAvarjV0fYNBSo1Wp0Oh3r1q1j6tSpgHWG8rp167j99tttOTRBOD9UlTTddrcuBFRb77sW27uSoAkkwbMHCZ792VOdh4SEt6MPOq2Oy7W6U4aAEykUciwW6bTHCMLH6f/jnZS3uX3AHcwOu+KsHivLSba2LbZzxXjVd+DWq51G2XqVNWZe+GMvn/x3mOF9PHlmWl+8Xe2pqamx6bhsfvvgnnvu4brrriM+Pp5BgwbxyiuvUF5eXr8aQRCENmCusW7AU/euPzfN2o+/5CgAklyJ5BlKkVcoet8wEmRGEqqOsa/kIBJF+Cjs0WkHM8tbh04bj6/T2W8Fq1AosEiW9nh1Qjfx9b6veDXpFW7oeyNzo87uGiA79C+qr69G8uhjbVvsdPatjztK8tFiFq1J5lhRFY9cEsGVg/w7ze6hNg8Fs2fPJjc3l0ceeYTs7GxiYmL49ddfT5p8KAhCC0gSlBxF3ujib92AZ3/DBjyuflg0kZj7zqCoRyB6lYS+Mht9bhL7inYiVUn4OPoQ7z2QOZHXWUOAs2+rh6ZQKrGYzWc+UDgv/XjgR55NWMacsDnMi55/Vo+V7/0Z5dqbkfwHY5zxEdid/W6JHcFktvD2PwdYuTGTCB8X1s4bSh9N52qMZfNQAHD77beL2wWCcLaqihu23TWkWX+fm15f+pfsXKw77vUajDnueiRtFIWuvuwo2U+iQY/eoGdf+ndISPg6+RKn1TEnfE6bhYATqZRKTCIUCM1Yd3gdT/z3GFP7TOOeuPvO6l2zfOdnKH++G0v4pZimvNlpd5c8lF/B4jXJ7DpazG0jg5nfyfY8qNMpQoEgCKdhrrFutVt34a+7/39C6V/SRGIJmYCkjcSiiQJXP4pqitlhSERv0KNPeZl9RfsA8HXyRafVcWX4VcRpde0SAk6kUqsxmUztfh6ha9l0bBMPbl7KhQETWBr/wFkFAsXW11Gufxxz7FxME5+zWdvi05EkqX53RK2LPZ/fNIgYf3dbD+uURCgQhM5CkqDkSMO7/rpf8/chs1gvppKrHxZtFOa+l9dvvSt5htRv6lJUXUiiIRH9/s/Q5ySwv3g/AH5OfsRpdVwVcQ06bRw9ndo/BJxIrVJjMtq2qZDQuSTkJHD/pkVc0HM4jw95HEVLL+qShGL9Yyj/W4Hpgnsxj1zSKdoWnyi3tJqHvktl4948Zsf34v6JYTjZde7LbucenSB0V1XF9Rf9htJ/GrJq61a8kp2r9YLvP6S+9C9pIsHercnTFFYVknjsH+vtgMYhwLkXOq2OqyOvtVkIOJHaTm3zmdVC55Gct4uFf99FrCaWZRc8i1KuatkDLSaUPy9EsetzTBc+g3ngLe070HP0Z7qBh75LRSGX8fZVsYwO19h6SC0iQoEgtKfGpf/GIaDU2p1EkquspX/tyaX/5t75FFYVkpibiD4nAb1BT8YJIeCayGvRaXXntMd8e7Ozs8dsNmM2m1EoOl+ZV+g4ewr3cOfGOwjvEc4LI15E3dLti42VKL+9GXnGnxgvewtL38vbd6DnoKzKxNO/7OabHccYH6nlySlReDh1cK/iVhChQBDaQuPSvyGtYelfwf5Gpf9eWLSRmPvNrH3nH4Xk2ee0+7kXVhVa5wMYrCEgszgDgF7OvdBp47ku8jritDp8nJpv9tWZODhYdySsqanBwaH7NzASmnegOJPbN8zH38WfV0a+ir2yhf8WqkpQfXUVsuNJmGZ+gqXP+PYd6DlIOFjI4m+SKaww8szUvkyPPfulu7YmQoEgnK3KovqLvjy38az/E0r/AUMxx99o/b0mEuzP3N2uoKqgfmVA4xDg7+yPzjueuZFzu0wIOFFdEKiprhah4Dx1pOwI8zfMw8Pek9dGv4GzuoVLB8tyUK2ejazkCMYrv0HqNbB9B3qWakwWXt+QwbubDhDn785H18fj36PjWxS3BREKBOFUTNXW0n/9mv9Tlf6jsIROss7+10aBi2+LJz01hIAE9Dl6MksyAQhwCSBOq+P6qOuJ0+rwduz6fTscHKw/JKuqq3E7w7FC95NTkcP89bfhoHRgxZiVuNu5t+yBhQdRr54JxkqM1/yIpIlo13GerX2GMhZ9ncz+3DLuGR/KjRcEopB3repAYyIUCIIkQfHhkxv+FGQ0lP7d/K0Nf6JnIWkaz/pv4eSoWvmV+STm6tHnWIPAgZIDgDUE6LQ6buh7I3FaHVpHbZu/TFtzdKwNBVVVNh6J0NHyK/OZt/42JCRWjnkLL4eWdRuUGVJRfT4Tyc4Z47U/g3tAO4+05SwWif9tzeLFP/fR28OBr24ZTGTPrr/XhQgFwvmlcenfkNpQ+q8pA0Cyd7Ne9AMuwBx/01mV/puTV5lHokFvXSbYJAT0RqfVcWPfm7ptCDiRvb09crmcqspKWw9F6EDF1cXcvnE+FaZy3h33XotvfcmytqD66iqkHoEYZ38BTp1n9v7x4iqWfJPC1gMFXD+sNwvHhWBng82L2oMIBUL3VFf6r530J6+rAJQeB2pL/15h1pJ/2EXWWwCaKHDp2ar1znUhoG5y4MGSgwD0dgkkThvHjX1vRqfVoXHsPD/gOopMJsPJyYlKEQrOG+XGcu786w4MFbm8O/49/F1a9k5fvu9XlGtvQvIbiPHy/3WatsWSJPFjcjaP/5iOk1rBqrk6hgZ72npYbUqEAqFrkyzW0v9JDX/2I5OsLXXrSv+W6CusVQBtFJJHn7Mu/TcnrzIXfW0VQJ+j51DpQcAaAnTeOm7udys6bRxeDudfCGiOs7MzFRUVth6G0AGqTJUs/OsuDpYc5O2xbxPsFtyix8l3rUb5011Ywi7CdNlboLRv55G2TFGFkcd/TOfnlGwm9+/JI5dE4OrQ+p8hnY0IBULXUVnYfMOfmnKgrvQfhdR7OOb4mxsa/rThu4zcitz6KkCiIbE+BAS6BhLvHc+t0bcSJ0LAKbm4ulJaUmLrYQjtrMZcw+JNi0grSGPFmDeJ8Ihs0eMU/61Aue5RzDHXYJr0QqdpW/zv/nyWrE2hymjm5Zn9uTi6663+aSkRCoTOx1SNLH9voy1+05EZUpGVZQMgKdRInmHWRj9tWPpvTuMQoDfoySo9BECQa1CjEKBr8cSp852rqys52dm2HobQjkwWEw9uXkpCTgKvjHqVAZoBZ36QJKHY+CTKLa9hGrYQ86gHOkXb4iqjmRd+38fH/2VxQR9Plk3ri7dr56hctBcRCgTbOan0n1rb8CejUek/AIs2Ekv/OW1e+m+OocLQsETQoCerNAuwhoBB3oOY138esZo4EQLOkZurG2VlZUiS1OWaughnZpEsPL71Uf4++g8vjHiBQT6DW/AgE8pf7kWx81NM45/CPOi29h9oC6QcK2HR18kcLarkoYsjuGqQP/IuvNSwpUQoEDpGRUHzDX/qS//u1ot+4AjMA29B0va1rkdu5wlGhgpDfQBIbBQCgl2DGewzmPn9FxCricPToXtNJrIVN3c3zGYzVVVVooFRNyNJEsu2P8NvWb/x9LBnGOE38swPMlWh/PYW5Pt+wzh5JZboWe0/0DMNyWzh3U0HeWNDBuHezqydN4Q+GmdbD6vDiFAgtC1TFbK8vU23+DWkNVP6j8ISfgkWTRSSNgqcfTqkXJhTkdOwOiAngcNlhwEIdusjQkAHcHfvAUBpaakIBd2IJEm8vOMl1mZ8w6ODH+PCgAlnflB1KaqvrkZ2LBHTzI+xhLTgMe0sq6CCxWuS2XmkmFtHBjN/VDBqpdzWw+pQIhQI50ayQFHWCQ1/miv9R2Hpf2XDNr8ewe1W+m9Odnl2kyWCR8qOANYQMKTnUBZo7yBOG4eHvUeHjel85u7uDlhDgVbb/XsznC/eSXmLz/Z8ymLd/UwOnnLmB5QZUH0xG1lxFsY5XyP5t+A2QzuSJImv9EdZ9usevJzVfHbjIGID3G06JlsRoUA4sxNL/4Y0ZHm7G0r/Dj1qS/8jMQ+6tbbhT/uX/puTXZ7d5HZAXQjo4xbCsJ4XoPOOJ04TRw/7Hh0+NgHUajVOzs5iBUI38r/0j3g35V1uH3Ans8Jmn/kBRYdQfX45MmMFxqt/sFYKbSivrJqHvk1lw948Zun8WDIpHCe78/fSeP6+cuFkTUr/qQ0Nf8pygNrSv1e4teFPxKUdXvpvTnb58frNg/QGPUdrQ0CIWwgX+A4nTqsTIaCT8fTwoESEgm7hq31f8lrSq9zY9ybmRs094/EyQxqq1TORVI7UXPszuPdu/0Gexp/pBh7+PhUZMt68MoaxEaJ6JULB+aiu9N94i9/cNGQFmQ2lf/fe1oY//a+ylv61fa2lf7lt/8lklx8noXY+QKJBz9Hyo4A1BAz3HY5OqyNOG4e7nQgBnZWnpxcHDmTaehhCK/2Y+QPPJTzLnPAruS163hmPlx3+D9VXVyK5BVjbFjvb7gJcVm1i2S97+DrxKGMjNDw1JQpPZzubjaczEaGgu6vIb5jwVx8AdiMzNi39W4JGIw2aZ13y5xXeadqKHi8/Vr95UKIhsT4EhLqHMsJvJDqtjlhtrAgBXYiXlxc7dyZhsViQy8+vSVzdxZ9Zf/LEtseZ2mca98Tee8blpfL9v6P85kYk31iMl39yznuJtAX9oULu/yaF/PIanp7alxmxvmJ5bCMiFHQXpipkeXuaNvxpUvq3s/b610ZhiZiMpa7bnw1L/805VnasSZ+AY+XWbYrD3MMY6TeSOK2OWG1cy7ddFTodL40XFouF0pIS3GonHgpdx6aj//Dg5geYEDCRpfEPnDkQpHyJ8oc7sIROxDT1HZu1La4xWVixMYN3/jlAjL87H16nw9/D0SZj6cxEKOhqJAsUHWq46BtqPwozkUkW6yHugdaGPwOubtTwx/al/+YcKztWHwD0hgSOlx9HhoxQ9zBG+Y1C5x1PjCZWhIBuRKOxlo0Li4pEKOhitmdvY/GmRYzwG8FjQx5DcYY2xIptb6H88yHMA67CdNGLNvsZtN9QxqI1yezNKePucSHcNDwIxXnQiOhcdL6rhNCgrvRfd+8/Nx1Z7p5GpX8Pa+k/eAySZn6nK/2fSJIkjpUfq988KDFX3yQEjO41Bp02nlhNLG52brYertBOHB0dcXJyorCwkMDAQFsPR2ihXXk7ueefhcRp43hm2DKU8tMsLZYkFH89g3Lzy5iG3ol59MM2qUhaLBIf/5fFC3/sw7+HA1/dOpionra7ddEViFDQGRgrraX/Exv+lBuAE0v/U7Bo+yJpI8HJu1OV/k90YgjQGxLIrshGhoywHuGM6TXWejtAhIDzjrePD4UFBbYehtBCuwt2c+fGO4joEcELI15ErVCf+mCLGeWvi1Ak/Q/TuMcxD17QcQNtJLu4iiVrU9iSWcC1QwK498JQ7FWdY4OlzkyEgo4kWaDwYMNFv+7X5kr/MdfU7vIXheQR1ClL/yeSJImj5UfR5yTULxHMqQ0B4T3CGes/zloJ0MbiqhZp/Xzm4+3Njh07xB4IXcCB4kxu3zifAJfevDzqVeyVp+lEaapG+d1tyPf+jPHSN7D0v6LjBtrIj7uO8/iP6TioFXx4nY5hfUSH0pbq/Fearqo8r+mEv9y6Wf/WveQlBw/rO//gsUjaBbUNf8JB3XV6bEuSxNGyI9YlgrWrAxqHgPH+49F564jRiBAgNOXTsyeVmzdTUVGBk5OTrYcjnMKR0sPM23AbXvZevD76DZxVp/n5VF2K6utrkR3djunyj7CETuq4gdYqrjTyxI/p/JiczaXRPjxyaSRuDh3XQbU7EKGgtU4s/ddVAOpK/0p7a+lfE4Ul8rLahj+dv/TfHEmSOFJ2uL4KkGjQk1ORg1wmJ7xHOBcGXEicNo5YTRwu6s45r0HoHHr29AUgPz9fhIJOKrs8m3kbbsNR6cSKMW+e/hZfeS6qL65AVngA4xVfIQUM7biB1tqckc+StSlU1ph58fJoLu3fs8PH0B2IUNBSdaX/Exv+FB5AJlmQkEGPQGvDny5Y+m9OkxBQe0vAUGlALpMT0SOCCwMmoNNaKwEiBAhnw8XFBWdnZ/JycwkICLD1cIQT5FfmM3+DtSHRyjFvnn6DsOLD1rbF1aUYr/4eybtfB43Sqspo5sU/9vG/rVkMDfbg2Wn98HGzzbLH7qBrXq3aW3lu0y1+DWnWakB96d8TSRuJJXhcQ7c/r3BQd+13PJIkcbjscH0ASDwhBEzsPRGdNp4YTQzOIgQIreTr50deXp6thyGcoLi6mAUb5lFhKue9ce/j4+RzymNlubtRfX45KO2tbYt7BHbcQIHUYyUsWpPMkcJKHrw4nKsHBSAXSw1bRYSCZqjfH4OsLLuh9K/tiyVqWkPDHydtlyv9N0eSJLJKs+rnA+gNCeRW5taGgEgm9p5k7RPgNUCEAKHN+fn1ImP/ftHZsBMpM5Zxx8bbyavK451x79HLxf+Ux8qObEf15Rwk114Yr/gCnL07bJxmi8S7/xzg9Q0ZhHk7881tQwjRdp35WJ2ZTUNBYGAghw4davJny5YtY8mSJTYakZVx6rvg5IXUIxjO0JyjK2kcAuqaBeVV5tWHgEm9L6ptFhRz+glFgtAGevXqhdlsJj8/H41GY+vhnPeqTJUs/OsuskoP8dbYdwh2Cz7lsbKMdai+uR7Jpz/GmZ+CfcctKT5cUMGiNSnsPFLELSOCWDC6D2qlCJVtxeaVgieeeIKbb765/nMXF9u/I7XFJJn2IEkSh0oPNWkWlFeZh0KmIMIjkosDL0GnjWeAZoAIAUKH8/HxQalUYsjJEaHAxmrMNSz65z52F+5mxZiVRHhEnPJYeerXKH+4HUuf8Zimvguq0yxRbEOSJPF14lGe+WUPnk5qPrlhILreYs+TtmbzUODi4oKPz6nvWQktZw0BB0nI0dfuH6Anv8oaAiI9Irkk8FJ0Wh0DNDE4qbr2/Aeh61MoFPj5+ZGTk0Pffh07OU1oYLIYeWDzUvQGPa+Oeo3+XgNOeax8+7uo/liKuf8cTBe/3GGTqPPLqnn4+zTW7c5lps6PJZPCcbaz+eWrW7L5d/XZZ5/lySefJCAggCuvvJKFCxeiVJ56WNXV1VRXV9d/fj7vy944BFjnBejJr8pHIVMQ5RHFpUEiBAidW+/egWzZslnMK7ARs8XMY1sf45+j//DiiBcZ6DOo+QMlCcXfz6L890VMgxdgHvtYh82rWr/bwEPfpSEhsfLKGMZF2G7L5fOBTUPBnXfeSVxcHB4eHmzevJmlS5dy/PhxXnrppVM+ZtmyZTz++OMdOMrOQ5IkDpYcqJ8PkGhIbBICJgdNQeeto7/XABEChC6hd+/e/P33X2JegQ1IksSzCc/we9ZvPDNsGcP9RjR/oMWM8rf7UexYhWnso5iH3NEh4yuvNvHsr3v4Un+UMeEanrosCi9nuw459/mszUPBkiVLeO655057THp6OhEREdxzzz31f9a/f3/UajW33nory5Ytw86u+b/8pUuXNnlcSUkJ/v6nniHblUmSxIGSA/VVAH2OnoLqAhQyJX09o5gSPAWdNp7+XgNwVIktQIWup6evL2q1muPHj4tQ0IEkSeKlHS+yNmMtjw1+nPEBFzZ/oKka5ffzke/5AeMlr2EZcGWHjC8xq4jFa5LJL6/hqcuiuDzOT7TD7iBtHgruvfde5s6de9pjgoObn9U6ePBgTCYTBw8eJDw8vNlj7OzsThkYurrGIUBfOy+gIQT05bI+U9FpdSIECN2GXC6nd+/eHD92jP79+9t6OOeNt5Pf4vM9n3F//BIuDZ7c/EHVpajWzEV2eCumGauwhF3c7uOqMVlYsTGDd/45wIBebnxwnY4AD/GzriO1eSjQaDTnnPiTkpKQy+VotefHPSNJksgsyWzSLKiwuvCEEGBdHeBwuk1IBKELCw7uw++//0ZNTQ1q9Wl23xPaxKq0VbyX+i53xtzFzNBZzR9UkW9tW1ywH+MVXyL1vqDdx5WRW8air5PZk1PGnWNDuHl4IEqFmGfS0Ww2p2DLli38999/jBkzBhcXF7Zs2cLChQu5+uqr6dGjey4zsUgWMosz6+cD1IUApVxJX49+TAuZXns7oL8IAcJ5I7hPHyRJ4vjx4/Tu3dvWw+nWvtz7BW/sfI2b+t7MtZHXNX9Q8RFr2+KqYoxXfYfk074VHItF4pNtWbzw+z783B344pbB9PMVG6jZis1CgZ2dHatXr+axxx6jurqaoKAgFi5c2GS+QFfXOATUVQKKqotQypX08+zH9JAZxGl1IgQI5zU3Nze8vLw4euSICAXt6MfMH3he/xxXhl/FrdG3NXuMLG+vtW2xQkXNtT+Bx6kbGLWFnJIqlq5N5d+MfK4ZEsB9F4Zir+o+DeO6IpuFgri4OLZu3Wqr07cLawjIsG4lnGOtBhTXNISAGSGX184J6H/6PckF4TwTEhJKUtIOsTSxnfyZ9QdPbHucaX2mszD2nmYn7cmO6lF9eQWSc0+Mc74E5/btH/NzcjaP/ZiGnVLBB9fquCDkNJsuCR3G5n0KujKLZCGjaH+jrYQbQkC0ZzQzQ2ei0+qI9ooWIUAQTiMkNJStW7eQl5d33swp6iibjv7Dg5sfZGLARJbEL20+EGRuQLVmLpJ3P2vbYgf3dhtPcaWRJ35K58dd2Vzcz4dHL43E3VHVbucTzo4IBWehaQioqwQUo5KriPaMZlbYTOI0IgQIwtny9fXF0dGRw1lZIhS0oe3Z21i8aREj/Ubw6JDHUTSzl4s8bS3K7+djCR6Nadr70I4rm7Zk5LNkbQrlNWZeuDyayf17ttu5hHMjQsFpWCQL+4v2188J2HFSCJiFThtPP89o7JVi/25BOFdyuZzQsDAy9u8nTqcTa9LbwM7cndzzz0J0Wh1PD1uGspmWxHL9Byh/ux9L9CxMF78CivZ5x15tNPPSn/tZteUQQ4I8eHZ6P3q6iZ+ZnZEIBSeQJIkv9q4mobZhUElNCWq5mn5e0cwKm10bAvqJECAIbSwiIpKdSUkUFhTg4SnuL7fG7oJ07vrrDiJ7RLJ8xAuoFScs9ZQkFP88j3LTckyD5mEe9zjI2mcuR9rxEhZ9nUxWYSVLJ4Vz7ZAA5HIR+jorEQpOIJPJ+OXQL9gr7LkibA46b2sIsFN0z4ZJgtBZBAQEYG9vz8GDB0UoaIXM4kwWbFhAgEtvXhr1ysm3MiULyt+XotC/j2n0w5iH3tku+xiYLRLvbzrIaxv2E6Jx5pvbhhCqFbuxdnYiFDTjwwtXIW+n1CwIQvMUCgXh4RFkZOwnNi5O3EI4B0dKDzN/w21oHLx4ffQbJ2+Jbq6xti3e/T3Gi1/GEnNNu4zjcGEF969JIfFwETcPD+KOMX1QK8XP1K5AhIJmiEAgCLYR1TeKnTuTyMvLE3shnKXs8mzmbbgNR6UTK8a8iZudW9MDaspQrbkeWda/mKZ9gCXi0jYfgyRJfLPjGE/9vJsejio+uWEg8b27ZzO67kqEAkEQOg1//wCcnZ05kJkpQsFZyKvMY/6G25Ah482xb+HpcMLtl4oCVF/OQZa3B+PsL5ACT7EjYisUlNfw8Pdp/JluYEacHw9MCsfZXlxiuhrxNyYIQqchl8uJ6tuXnUlJ6OLjUShEd7szKaouYsGG+VSaKnl33Ht4O3o3PaDkKKrPZyKrLMB49XdIPgPafAwb9uTy4LepWCSJFXNiGB8plpV2VaJOLghCp9KvXzTV1dUcPXrU1kPp9MqMZdy58Xbyq/JYOeZNerk03UZelrcP9f8uRmaqxHjNj20eCMqrTTzyfRq3fbqDaD9XflgwTASCLk5UCgRB6FS0Wi3e3t5kZmQQEBBg6+F0WpWmSu7+604Olx7mrXHvEOTWdJ8C2bFEVF9cgeSkxTjnK3Bp20ZBO7KKWPxNMrml1TwxJYpZOj8xObQbEJUCQRA6nf4DBnDkyBEqKyttPZROqcZcw33/3Muewj28Nvp1wnuEN/m67MBfqD6dhuQZgvGaH9o0EBjNFl5Zt58r39+Gh6Oa7+YPY3Z8LxEIugkRCgRB6HT6RvVFLpeTkZFh66F0OiaLkaX/LiEpdwcvj3yFaK+mWxvL079D9cUVWAKGYpzzNTi03ez/jNxyZr+7jXf/OcAdY0L49MaB9PZsv7bIQscToUAQhE7H3sGByMhI9u/bhyRJth5Op2G2mHls62P8e3wTzw1fTrz3wCZflyd+iHLtTVgiL8N0+cdtto+BJEl88l8W097cQkWNidU3D2L+6GCUCnEJ6W7EnAJBEDqlmNg4UlJSOH7sGL5+frYejs1JksSyhGf4Pes3nhn2LMN9hzf+Iop/X0T597OYBt6KefyTbda2OKekigfWprIpI5+rB/tz34VhOKjFqpDuSoQCQRA6JT8/PzQaDXv27D3vQ4EkSbyU+ALfZqzlsSFPMD5gfKMvWlD88SDKhHcxjXoA87CFbda2+JeUbB79IQ21UsF718QxItSrTZ5X6LxEKBAEoVOSyWTodPH89tuvlJWV4ex8/vbNfyv5TT7f+zlL4pdyaVCjToTmGpQ/3oE89RuMk17AEje3Tc5XUmnkyZ938/3O40zq681jkyPp4ag+8wOFVlOpVNxz7331v+9o4oaQIAidVlTfvqjt7Nize7eth2Izq9I+5P3U97gr5m4uD53Z8IWaclRfXYN89w+Ypr3fZoFga2YBU1ZuYcOeXJbPiOaVWf1FIOhAMpkMtVqNWq22yYoOEQoEQei01Go1MTEx7N+/H6PRaOvhdLgv9q7mjZ2vc3O/W7gm8tqGL1QWovr8cmRHtmKctRpL5JRWn6vaaObZX/dw3aoEAjwc+H7+UKYM6CmWGp5nRCgQBKFT0+niMRqN7N+/39ZD6VDfZ37Hcv3zXBV+Nbf0u7XhC6XHUX08GVlBJsarvkUKGtnqc+3OLuXyt//jk/+yWDIxjFXXxePr7nDmBwrdjphTIAhCp+bq6kpkZCS709MJDw9HLu/+72X+yPqdp7Y9yfSQGdwdu7D+3bosfz+qzy8HwHjtj0ieoa06j9ki8cG/B3l1/X76aJz55rYhhHm7tHr8QtclQoEgCJ3e4CFDSUtL49ChQwQFBdl6OO3qn6N/89Dmh5jYexJL4pc2BILjSahWz0Zy8sJ4xVfg6tuq8xwprOT+b5LRZxVx4wWB3DU2BLWy+wcu4fREKBAEodPz9vYmMDCI1JQUAgMDu+197m3Z/3H/psWM9BvBo4MfQ17ba0B28G9UX1+DpInEOPMzcPQ453NIksTapGM89fNu3B1UfHx9PAMDz/35hO5FxEJBELqEocOGUVhY2G13T0zKTeKevxcS7x3P08OWoZRb37PJd/9g3dio12CMc9a0KhAUlNdwx+qdLF2bysQob76fP0wEAqEJUSkQBKFLCAgIwM+vF8m7duHn17125NtdkM5df91JlGdfnh++HLXCugRQvuN/KH+9D0vkVEyT3wDFuS8N3Lgnlwe/S8VkkXj9igFMiPJuq+EL3YioFAiC0CXIZDIuGD6cvLw8jh87ZuvhtJmM4gwWbFhAkGsgL498BXulQ23b4pdR/XIPlrgbMF321jkHgooaE4/+kMatn+6gb09XflwwTAQC4ZREpUAQhC4jKCgIX19fdu7cSU9f3y5fLThcmsX89fPQOmp4bfQbOKmcrG2L/3wY5fa3MY24H/Pw+865bfHOw0Us/iaFnJIqHp8cKbY4Fs5IVAoEQegyZDIZI0aMJC8vj2NdfG5Bdvlx5q2/DWeVM2+MXomr2hXMRpQ/LECx/R2ME5/HPGLROQUCo9nCa+v3M+f97bg5qPh2/lCuGOgvAoFwRqJSIAhClxIYFESvXr1ISkrCt4vOLcirzGP+hnnIZXJWjn0TTwdPMFagXHsj8syNmKa+iyVq6jk9d2ZeOYvXJJN2vJT5o4K5bWSQ2OJYaDHxL0UQhC5FJpMxavRoCgoKOHTwoK2Hc9aKqotYsGE+laZKVo59C29Hb6gsQvX55cgPbcY4+/NzCgSSJPHptiymvbmF0ioTq28exO1j+ohAIJwVUSkQBKHL8fcPoE+fPiQlJeEfEIBCobD1kFqkrKaUOzYuoKAqn3fHv0cv517WtsWrZyEry8F41Vok37izft6ckioe+DaVTfvzuWqQP4smhOGg7hrfE6FzabcI+fTTTzNs2DAcHR1xd3dv9pisrCwuueQSHB0d0Wq1LFq0CJPJ1F5DEgShGxk1egxlZWXs27vX1kNpkUpTJXf/fRdHSo/wxpiVBLoGISvIQP2/S5BVFWO85sdzCgS/pmYzZcUW9mSX8u41cTxyaaQIBMI5a7dKQU1NDTNnzmTo0KG8//77J33dbDZzySWX4OPjw+bNmzl+/DjXXnstKpWKZ555pr2GJQhCN6HVaukXHc2uXbsI7tMHtbrzbu9bba7mvn/uZW/hXlaMWUl4j3Bk2busbYsd3DHO+QFc/c7qOUurjDz5026+23mciVFaHp8SJbY4Flqt3SoFjz/+OAsXLiQ6OrrZr//++++kpaXxySefEBMTw0UXXcSTTz7JihUrqKmpaa9hCYLQjYwcOQqz2UxycrKth3JKJouRpf8uISl3By+NfIVor/7IDm1C9ckUJHd/jNf8eNaBYNuBAqas2MK63bk8N70fr84eIAKB0CZsNgNly5YtREdH4+3d0ERj4sSJlJSUkJqaesrHVVdXU1JS0uRDEITzk4uLC0OGDmV3enqn/Flgtph5ZMsjbD7+L88NX068dzzyPT9bKwR+8Riv/AYcPVv8fDUmC8/9todrVyXQq4cD3y8YytSYrt+vQeg8bBYKsrOzmwQCoP7z7OzsUz5u2bJluLm51X/4+/u36zgFQejcBg8egpOzM/qEBFsPpQlJknhm+9P8efgPnh62jOG+w5EnfYLym7lYwiZhnPkpqJ1b/Hy7s0uZ8fZWPt6axaILw1g1Nx4/d4d2fAXC+eisQsGSJUuQyWSn/di9e3d7jRWApUuXUlxcXP9x+PDhdj2fIAidm0qlYty4cRw5coSjR47YejgASEisk/7ku8xveXTwY4zrNRbFltdQ/Xw3ltjrMF32DijtWvRcZovEe5sOcPnbWwH4+tYh3Dg8EIVcVAeEtndWEw3vvfde5s6de9pjgoODW/RcPj4+bNu2rcmf5eTk1H/tVOzs7LCza9l/JkEQzg/h4RH07t2b7du349Ozp02XKB6vPMbv8t9Jl9JZEv8AlwRejGL9oyj/W4lp+CLMIxa3uEvhkcJKlnyTQkJWITcMC+TucSGolaLvgNB+zioUaDQaNBpNm5x46NChPP300xgMBrRaLQB//PEHrq6uREVFtck5BEE4P8hkMiZMnMj7771HSnIyA2JiOvT8kiSRVpzKd0e/Y1v+Nuxl9kyUTeLy4MtQ/ngn8uQvME54Fkv8TS1+vu92HueJn9Jxs1fxv7nxDAoSWxwL7a/dliRmZWVRUFBAVlYWZrOZpKQkAEJCQnB2dmbChAlERUVxzTXX8Pzzz5Odnc1DDz3EggULRCVAEISz5unpxeAhQ/hv61YCg4Jwc3Nr93OaLCY25W7i+6PfkVGWQS9Hf+aHLkCWDq5yNco11yHPXI/psrew9J3RoucsKK/hsR/S+C3NwLQYXx68OBwXe1U7vxJBsJJJkiS1xxPPnTuXjz766KQ/37BhA6NHjwbg0KFDzJs3j40bN+Lk5MR1113Hs88+i1LZ8qxSUlKCm5sbhtw8XF1d22r4giB0QUajkfffexc7e3suvPDCdpuVX2os5bfjv/HzsZ/Ir8knpkcsl/lNIbZHHDKZjH///ImJ+R+gMR/HOGMVUvDYFj3v3/vyeGBtCkaLxBOTo5jYV2xxLLSNkpIStBoviouLT3utbLdQ0FFEKBAEobGDBw6wevXnDBk6lNDQ0DZ97mOVx/jh6Pesy16HRbIwyns0U/ym0NupNwDK0qO47v0au7SvkCMhu2YNkl/8GZ+3osbE8t/38dm2w4wI8eSZaf3QuoiKqdB2WhoKxN4HgiB0K4FBQfSLjiZRr8fPzw9HR8dWPZ8kSaQWp/Ld0W/Znr8dV5Ur0/2nM6nnRbir3cFixjFrI657v8Lh6GYsamfSHWPZ53khF7cgEOw6UsyiNclkl1Tx6KWRzBnYS/QdEGxGhAJBELqd8ePGcyAzk/+2bmX0mDHndJE1Woz8m7uJ745+R2ZZJv6O/iwIXcAo79Go5WoUFbm4pL+L6741KMuzqfLsS+6wRykPmsh/f2/BTmV/+uc3W3j77wOs/CuTqJ4uvHnVUIK9nM71JQtCmxChQBCEbsfewYGJky7imzVfk5mZSZ8+fVr8WOt8gV/56dhPFNQUENsjlsf6PU5MjxhkgH32dlz3fIlT1kYkuZKy4IsoCZ9JjWfLV00dyCtn8ZpkUo+XMm9kELeNCkYltjgWOgERCgRB6JbCwsLo27cf27dtw8fHByen078LP1pxlB+O/sD6HOt8gdHeY5jiN4UApwDk1SU4p3+K656vUZccpMYtmPyB91LW5xIs6pbPZZIkic+3H+G53/bg42rP5zcOZIC/eytfqSC0HREKBEHoti6cMIGsrENs/vdfxjezGkGSJJKLk/n+yPckFNTNF5jBpJ6TcFe5YZeXguuOR3A6+BsyyUx5wDjyhj5ElbeuxQ2I6hhKq3nw21T+3pfHnIG9WDwxDEe1+BEsdC7iX6QgCN2Wvb09l1x6Kas//5z09PT6xmhGi5FNuZv47si3HCg/QG/H3twedgcjtSOxM5txPvAzrnu+xq4gHaNTT4r630xp6DTMDi3fvKix39NyePj7NFRyGe9cHcuosLZpAicIbU2EAkEQurXAwCAGDhyEXp+As5cz/1Vt5adjP1NYU4Cuh465wXMZ4B6DuigD1+0v4ZLxIzJjORW9hlMQ+xqVvheA/NzaJpdWGXn65z2sTTrGhCgtj0+OwsNJbHEsdF4iFAiC0G2ZLWb2Fe3lmO8xth3fxpr0NSjkCkZ7j2ay3xQC7H1wOrQO1/9uwiFHj8neg+KI2ZSGTcfk7Neqcx8oV3DZyi0UVRp5dlpfscWx0CWIUCAIQrdRFwISDHr0hgR2GBIpM5Zhp7AjUhuJ2wE3RvYYyaSeUbjsXoPr/rUoqgqp9NaRM/I5ygPGguLcWwqbLRJJx8r5PtuVXaUO6Hrb87/rB9Krh9jiWOgaRCgQBKHLMlvM7C3agz6nNgTk7qgPAQO8BnB1xLXoPPrSzyLHLn8/BQV/UZnxCv4Zx7ConCgLmUxJ2OUY3Vu+ZPHkMUjsOl7BxswS/sksoajKTA+Vmkk+Vbx4/UCxxbHQpYhQIAhCl2GymNhbuBe9IYFEg57E3B2UG8uwU9gzwCuaa3pfQjyO9KsowS53L7I9byMrPIgMCQkZGo8gjrv68FNZP/wuvAs3je85jcMiSaRkV7Aho4S/M0sorDTj7axiYrg7o/u4cnTXFuzt7UUgELocEQoEQei0GocAayUgyRoC5HbEuPRmrnM48VU1RBccRn3wZ2SmSgAkRw2SNhJLyAQkbRSSJhLJKxzUTrgZjRz6aBX7/t3GRRdfjFrdsol/FkkiLaeSjRkl/JVZQn6FCY2TknEhbozp40qE1qF+zsAxkQWELkqEAkEQOg2TxcSewj21IUBPkiGRclMF9jIlMUo3bqhWEl9YSXRJFir2ISkdkDThSNp+mPvNxKLpi6SNBKdTL/lTqVRMnz6DVR9+wJbNmxk5atQpJwBKksTu3Co2ZhTzV2YJhjITno5KRgW7MrqPK1HeDsjF5EGhGxGhQBAEmzFZTOwu3I0+ezuJxzaRVJBKuaUGe2TEGiVuKC1iUFUVUdVGlB5BSJoopAGTkDRR1GgjkdwDz2m5oIeHB5OnTGHN11+TkpJCdHR0/dckSWJfXhUbMkrYmFFCTpmRHg4KRga7MjrYleiejiIICN2WCAWCIHQYk9nI7mNbSMz6E31uEkmVxynHjL1FIq6qipuqqonHgcgeYSh8+1pDgNZa+jeqWrfb4YlCQ8O44ILh/PvvJtzd3alx1LBhfzEbM0s4VmLEzV7ByCBrRaB/T0cxP0A4L4hQ0IxbPk4kp7Ta1sMQhC7NTqoiwJyJi3wXVcr95Khz2WdnpFIuw8FiYUCVkak1jvQwBmAyhXFYEcwfit58LXeHIqwfewGqgJ3tNEolxVIsX/2ST5lUiqudguFBLtw9wpVYXycRBITzjggFzejn50qvCqOthyEIXYJcMuFRfRSfqky8qvZRbtpDFsdIVVWzzd6OCrkce4tEmMmRcaZAeij7IrcbRJGrP6UyBaW1z6MAQmwwfrOfK/v27sZPnsNNU0bi7Ch6CgjnLxEKmnHnWFv8aBKETk6SoDwHmSEdeW4aMkM6ptxU0ksOkKCS85e9HYkO9lTKZTjKVMQ4BnOjdzy63hOI0A5AKT/3pkDtraioFx+t+pB///6L8RdeiEJxbm2NBaGrE6FAEIST1ZQhy92DrPbibw0BaZgqC0i1U7PdyYXtzu4kOZipdPDAUW5HjFd/buo5FJ13PBE9IlDKu86PF3d3dy6/fCafffYpm//9l+EjRoiWxMJ5qev8rxUEoe1ZTMgKDtRe/NOQ5aYhN6QjKzoIQI1MTrJXb7a7aUnoHUKSMZ8qixEnpRMDNDHcpNV1yRDQHL9evZg8eQrffrsWJ2dn4uLibD0kQehwXft/sSAILVNf+k9DbkhDlptuDQF5e5GZrZNqJSctNdpIkoMGk6C+gARjATtLMqgyV+EkKyHGM5ZbtLPRaeMJ7xHe5UNAcyIiIxlbWsL6detwcnQkPCLC1kMShA7V/f5XC8L5rqYMWe7u2nf+6Q0hoLIAAEnliKSJQPLpT3W/mSQ7OJJgLkFfmM7OvJ1UFx3ASeVMrCaGW3vfhk4bT1iPsG4ZApozcOAgSopL2LZtG/b29vQODLT1kAShw5wf/8sFoTuymJAVZNZf/E8s/UsyOVKPYCRtFJb4m5G0UVR5hZJmLiPBkEiiIYGdhz6h2lyFk8qZOE0st0WffyHgRDKZjHHjx1NeUc6mTZtQq9X09D23PRIEoas5P//XC0JXIklQlt30XX8zpX9JG4U5/OKGhj+eYdTIFaTmp6A3JJJ49Ad27tzVJATMi56HzltHmHs4inPoDNhdyWQyLr10MtVVVWzcuJHxF16IRnPq1smC0F2IUCAInUmj0r+8duKftfRfCDQq/fccgKX/HCyaSCRtFDh6Wh9urqkNAdvQp61kV94uqs3VOKucidXGMa//fOK1OkLdw0QIOAOFQsG06TP4YvXnrF+3jgsvvBAPT09bD0sQ2pUIBYJgC41L/3Wz/nPTkRUdAmpL/x59kDRRWAaOtP6qjQT33iCT1z9NjbmGlPxk9Jlr0OckkJyfXB8C4rQ65vVfIEJAK6hUKi6fOYvVqz/nzz//ZMLEibi7u9t6WILQbkQoEIT2VFf6N6TVN/yR5aY3Lf07eyNpojCHX9Kk9I/q5M561eZqUvJTSDTom4QAF5ULsdo45vdfgE4bT6h7qAgBbcTe3p7Zs6/gs88+5c8//uDCCy/ETQQDoZsSoUAQ2kp1qbX0Xzfh76TSv1PT0r82CkkTWV/6b/Ypa0OAPse6lXBy3i5qLDW4qFyI0+pY0P92dN7xhLiFiBDQjhwcHJgzZw6fffYZf/zxBwNiYlDI5ac8vqKiAjs7+w4coSC0DREKBOFsWUzI8jPqu/1ZQ0AasuIs4MTS/ygkbSQWTRS4BzQp/Ten2lxNcl4yeoM1BKTkJVNjqcFV7UqcVsftA+4QIcBGHB2duPLKK/ns00/ZumXLGY/v1y/6jMcIQmcjkyRJsvUgWqOkpAQ3NzcMuXm4urraejhCd3JS6b82BOTvRWausR7i7GNd8lc74U/SRCJ5hYGyZe8Sq0xV1jkBhgT0hsSTQoBOq0OnjSfEPQT5GQKF0DEsFgst+bEpl8tFq2Sh0ygpKUGr8aK4uPi010pRKRAEaCj9G1KtE/7qlv5VFQEgqZ2sF3zfWMwDrrLe99dEgaPHWZ2mSQjI0ZOcn4zRYsRN7UacNo47Y+5Cp9XRR4SATkt+mtsGgtDVtVsoePrpp/npp59ISkpCrVZTVFR00jHNpejPP/+cK664or2GJZzvzEZkBRknN/ypL/0rrKV/bSSWoNFnVfpvTpWpiuS8XegNevSGBFLyU2pDgDtx2jjuirlbhABBEDqNdgsFNTU1zJw5k6FDh/L++++f8rgPP/yQSZMm1X8ulvsIbUKSoPR4Q9k/N/2UpX9zxGRr6V8bheQZ2uLSf3OqTJXsahQCUvNTm4aA2IXEa3UEu/URIUAQhE6n3ULB448/DsCqVatOe5y7uzs+Pj7tNQzhfFBdWt/l7/Sl/zjMMVdbf6+JPOvSf3OqTJXszNtlXSJYWwkwWUy427kTp9Vxd+xCdNp4gt2CRQgQBKHTs/mcggULFnDTTTcRHBzMbbfdxvXXX3/ayTnV1dVUV1fXf15SUtIRwxQ6gxNL/3X3/4sPA7Wlf88QJE0kluAx1l+1UeDmf06l/+ZU1lcCrHMCUgsaQoBOG8/C2HtECBAEocuyaSh44oknGDt2LI6Ojvz+++/Mnz+fsrIy7rzzzlM+ZtmyZfVVCKGbOrH0XxcC8vc1lP5delob/kRe1rThj9KuTYdSaapkZ+7O+iWCaQWpmCwmetj1IE6r457Ye9F5xxPsGixmmguC0OWd1ZLEJUuW8Nxzz532mPT0dCIa7UG+atUq7r777mYnGp7okUce4cMPP+Tw4cOnPKa5SoG/v79YkthVVZc2tPlt3PCnqhgASe1cX+6va/YjaaPAoUe7DKfCWMGuvJ21cwL0pOanYpasIUCnjUfnbV0iGOQaJEKAIAhdRrssSbz33nuZO3fuaY8JDg4+m6dsYvDgwTz55JNUV1djZ9f8Oz47O7tTfk3oxMxGZPn7Gy76dff/my39j60NAX1rS//td/FtCAEJtSEgDbNkwsPOA523josDL0an1REoQoAgCOeBswoFGo2mXbcPTUpKokePHuKi35VJEpQea7rFryHNWvq3GK2H1Jf+p9av97fO+m//v/cKYwU785LQ51gnBqYVpGOWTHjaexKnjePiwEtECBAE4bzVbnMKsrKyKCgoICsrC7PZTFJSEgAhISE4Ozvzww8/kJOTw5AhQ7C3t+ePP/7gmWee4b777muvIQltraqk6Vp/QxqyvBNK/9oopF4DMcdeVxsAItut9N+ccmM5O3OTSDQk1oaANMySuTYE6LgkaDLx3jp6uwSKECAIwnmv3ULBI488wkcffVT/eWxsLAAbNmxg9OjRqFQqVqxYwcKFC5EkiZCQEF566SVuvvnm9hqScK4al/4b3/8vOQKAJFcieYRYG/30GdfQ8KedS//NqQsBdX0C0gvSa0OAFzqtjkuDJ6PTxtPbpbcIAYIgCCcQex8IDepL/6n1G/1YG/40Lv37NprwF4mk7YvkEdIhpf/mlBnLmqwO2N0oBMR7xxOnjRMhQBCE857Y+0A4vaqS+nf+8rpuf7lpyKqtfR9OLv3XbvPr4G7TYZcZy0jKTUKfk0BibmJ9CPBy8EKnjeey4MvQaeMJcAkQIUAQBOEsiVDQ3ZlrGpX+G23zW3IUOKH0HzK+oeGPa68OL/03p3EI0Bv07C5MxyJZ0Dho0GnjmRo8FZ13PP7O/iIECIIgtJIIBd2FJEHJ0dqGP+kN9//z9zeU/l39sGgiMfedUb/eX/IMBYXaxoNvUFZTag0BdbcDCnfXh4B4bTzT+kwTIUAQBKGdiFDQFVUVN6z1r9/tL72h9G/nYl3m12sw5rjrG3r927j035yymlJ25O5Ab9CT2CgEaB206LQ6pofMQKfV0UuEAEEQhHYnQkFnZq6xTvIznNDwp3Hp3zPUWvIPubCh4Y+rX6co/TentKaUpNoQoDcksKdwT0MI8I5nesjlxGt1+Dn3EiFAEAShg4lQ0BlIEpQcadrwp67Xv8VkPaRx6V8bVdvwJ6RTlf6bU1pTyo7cxPpmQXsK9yAh4e3ojU4bz4yQmSIECIIgdBIiFHS0quL6i35D6T8NWXUp0Ezpv27Wv72bjQfeMiU1Jeww7EBvSCDRoG8UAnzQaXXMDJ2FzjsePyc/EQIEQRA6GREK2ktzpX9DGrLSY0Bd6T+sdtb/hIaGP5249N+ckpoSEg2JJNZuILS3UQiI18YzK2w2Om08vk6+IgQIgiB0ciIUtFbj0n9ttz+ZIR1Zwf5Gpf9eWLSRmPvNbNTrv/OX/ptTXF1cOzEwAX2Onn1Fe5GQ8HH0Id57ILNFCBAEQeiyRCg4G5VF9Rd9eW7jWf91pX9X6wU/YChm3Q1drvTfHGsISKztE5BYHwJ8nXyJ0+qYEz7HGgKcfW09VEEQBKGVRChojqnaWvqvX/OffkLpX2Wd9a+NxBIyEUkbZW344+LbpUr/zSmqLmJHo9sB+4r21YcAnVbHleFXEqfViRAgCILQDYlQ0Az1mwMbAkBd6T96Vu16/ygkzz5dsvTfnLoQULdEcF/RPgD8nPyI0+q4MuJqdNo4ejqJECAIgtDdiVDQDNPE55AcPGpL/91rk6Wi6sLabYT16HMS2F+8HwA/517EaeL+3979xjR172EAf4oXCij/lP4BBS2ay+bc3Dgqt0zNiARczDI3w6sbp8OL6HDJNrKII5PthYHAMpMRB/jiVnNzMxZf6LJldwlhii+GLgKT4QaZU/4IK+gmLReRInzvC+Fs1crkWnpaeD7JSTinP+nDNw08np62+Ptj21kCiIjmKJYCD8b/+rzWEbzm5u2baLo+8XTAPSVAMSrY/vgrUIwKzPPjNE5KRERaYymYZSZLwOQHCP08UQKWLFgCxbgGrzy+AylGBeb5Zo2TEhGRv2EpCHC/3f5t4umAuyXgiuNnAEDCggSkGBXsYAkgIqKHxFIQYO6WgMaJEtDkVgIU0xq8uvJVpBgVmMJNGiclIqJAw1Lg59xKQF8jrjivAAASIxKhGBXkrMxBilGBMdyocVIiIgp0LAV+5tfhX9F0vVH9AKGrzqsAgMSIpVCMKch5YhdLABERzQiWAo3dGL6hvlFQU3/jPSVAwa4ncqEYFRjCDRonJSKi2Y6lwMduDF9Ho/qOgRfQ4ewAACyNWAbFpOAfq3ZDMaYgNowlgIiIfIulYIZNloDJawI6BzsAAMsil0ExrkHuqjyWACIi8gssBV52/dZ19S2Dm/qb3ErAGtMa5D2ZhxSjgtiwWG2DEhER3YOl4BHdLQEXJopAI7oGOwEAlkgL1prWYs9Te/CMIYUlgIiI/B5LwTT13+r/w/sENKJrsAsAkBSZhFTzOux9ai9SDAoWhS3SOCkREdH0sBT8if5b/b+fCei7gO7/dgMAkqKWI9WciteeysczhhSWACIiCngsBfcQEXzV+R+PJeBvcVbkG19HijEFC0MXapyUiIjIu1gK7qHT6WC79E/odEGwxqUh36iwBBAR0ZzAUuDBvzb/G/p5eq1jEBER+VSQ1gH8EQsBERHNRSwFREREBIClgIiIiCbMWCno6OjArl27YLFYEBYWhuXLl6O4uBgul8ttXUtLCzZs2IDQ0FAkJCSgrKxspiIRERHRFGbsQsO2tjaMj4+juroaK1asQGtrK3JzczE0NIQPPvgAAOB0OpGZmYmMjAxUVVXh+++/R05ODqKjo7F79+6ZikZEREQe6EREfHVn5eXlqKysxJUrVwAAlZWVKCoqgt1uR0hICACgsLAQp06dQltbm8fvMTIygpGREXXf6XQiISEB/ddvIDIycuZ/CCIiogDjdDphNMTC4XBM+bfSp9cUOBwOLFz4++v9GxoasHHjRrUQAEBWVhba29tx8+ZNj9+jpKQEUVFR6paQkDDjuYmIiOYCn5WCy5cvo6KiAnl5eeoxu90Ok8nktm5y3263e/w+Bw4cgMPhULfu7u6ZC01ERDSHTLsUFBYWQqfTTbnde+q/p6cHmzdvRnZ2NnJzcx8psF6vR2RkpNtGREREj27aFxoWFBRg586dU65JSkpSv+7t7UV6ejrS0tJw9OhRt3Vmsxl9fX1uxyb3zWbzdKMRERHRI5h2KTAYDDAYDA+1tqenB+np6VAUBTabDUFB7icmrFYrioqKMDo6iuDgYABAbW0tkpOTERMTM91oRERE9Ahm7NUHPT09eO6557B06VIcP34c8+bNU2+bPAvgcDiQnJyMzMxM7N+/H62trcjJycHhw4cf+iWJDocD0dHR+PnKVURERMzEj0JERBTQBgcHsTzJgoGBAURFRT14ocwQm80mADxuf3Tx4kVZv3696PV6Wbx4sZSWlk7rfrq7ux94P9y4cePGjRu337fu7u4p/6b69H0KZsL4+Dh6e3sREREBnU6ndRzNTL5fQ3d3Ny++9BLO1Ls4T+/iPL1vNs9URDA4OIj4+Pj7nsr/o4D/6OSgoCAsWbJE6xh+g6/I8D7O1Ls4T+/iPL1vts50yqcNJvADkYiIiAgASwERERFNYCmYJfR6PYqLi6HX67WOMmtwpt7FeXoX5+l9nKmPPxCJiIiI/BfPFBAREREAlgIiIiKawFJAREREAFgKiIiIaAJLAREREQFgKQh4HR0d2LVrFywWC8LCwrB8+XIUFxfD5XK5rWtpacGGDRsQGhqKhIQElJWVaZQ4MBw6dAhpaWkIDw9HdHS0xzVdXV3YsmULwsPDYTQa8fbbb+POnTu+DRogjhw5gmXLliE0NBSpqan49ttvtY4UMM6ePYsXXngB8fHx0Ol0OHXqlNvtIoKDBw8iLi4OYWFhyMjIwE8//aRN2ABQUlKCtWvXIiIiAkajEVu3bkV7e7vbmtu3byM/Px+LFi3CggULsG3bNvT19WmU2LdYCgJcW1sbxsfHUV1djUuXLuHw4cOoqqrCO++8o65xOp3IzMzE0qVL0djYiPLycrz33ns4evSohsn9m8vlQnZ2Nvbu3evx9rGxMWzZsgUulwvffPMNjh8/jmPHjuHgwYM+Tur/Pv30U7z11lsoLi5GU1MTVq9ejaysLPT392sdLSAMDQ1h9erVOHLkiMfby8rK8NFHH6Gqqgrnz5/H/PnzkZWVhdu3b/s4aWCor69Hfn4+zp07h9raWoyOjiIzMxNDQ0PqmjfffBOff/45Tpw4gfr6evT29uLll1/WMLUPTesjCSkglJWVicViUfc//vhjiYmJkZGREfXY/v37JTk5WYt4AcVms0lUVNR9x7/88ksJCgoSu92uHqusrJTIyEi3OZPIunXrJD8/X90fGxuT+Ph4KSkp0TBVYAIgJ0+eVPfHx8fFbDZLeXm5emxgYED0er188sknGiQMPP39/QJA6uvrReTu/IKDg+XEiRPqmh9//FEASENDg1YxfYZnCmYhh8OBhQsXqvsNDQ3YuHEjQkJC1GNZWVlob2/HzZs3tYgY8BoaGvDkk0/CZDKpx7KysuB0OnHp0iUNk/kXl8uFxsZGZGRkqMeCgoKQkZGBhoYGDZPNDlevXoXdbnebb1RUFFJTUznfh+RwOABA/Z3Z2NiI0dFRt5k+9thjSExMnBMzZSmYZS5fvoyKigrk5eWpx+x2u9sfLwDqvt1u92m+2YIzfTg3btzA2NiYx1lxTo9ucoac7/9nfHwcb7zxBp599lmsWrUKwN2ZhoSE3Hct0VyZKUuBnyosLIROp5tya2trc/s3PT092Lx5M7Kzs5Gbm6tRcv/1/8yUiGav/Px8tLa2oqamRusofuMvWgcgzwoKCrBz584p1yQlJalf9/b2Ij09HWlpafddQGg2m++7cnZy32w2eydwAJjuTKdiNpvvu4J+Ls70z8TGxmLevHkeH3+c06ObnGFfXx/i4uLU4319fXj66ac1ShUY9u3bhy+++AJnz57FkiVL1ONmsxkulwsDAwNuZwvmymOWpcBPGQwGGAyGh1rb09OD9PR0KIoCm82GoCD3E0BWqxVFRUUYHR1FcHAwAKC2thbJycmIiYnxenZ/NZ2Z/hmr1YpDhw6hv78fRqMRwN2ZRkZGYuXKlV65j9kgJCQEiqKgrq4OW7duBXD3lG1dXR327dunbbhZwGKxwGw2o66uTi0BTqcT58+ff+ArZ+Y6EcHrr7+OkydP4syZM7BYLG63K4qC4OBg1NXVYdu2bQCA9vZ2dHV1wWq1ahHZt7S+0pEezbVr12TFihWyadMmuXbtmvzyyy/qNmlgYEBMJpNs375dWltbpaamRsLDw6W6ulrD5P6ts7NTmpub5f3335cFCxZIc3OzNDc3y+DgoIiI3LlzR1atWiWZmZny3XffyVdffSUGg0EOHDigcXL/U1NTI3q9Xo4dOyY//PCD7N69W6Kjo91euUEPNjg4qD7+AMiHH34ozc3N0tnZKSIipaWlEh0dLZ999pm0tLTIiy++KBaLRYaHhzVO7p/27t0rUVFRcubMGbffl7du3VLX7NmzRxITE+Xrr7+WCxcuiNVqFavVqmFq32EpCHA2m00AeNz+6OLFi7J+/XrR6/WyePFiKS0t1ShxYNixY4fHmZ4+fVpd09HRIc8//7yEhYVJbGysFBQUyOjoqHah/VhFRYUkJiZKSEiIrFu3Ts6dO6d1pIBx+vRpj4/FHTt2iMjdlyW+++67YjKZRK/Xy6ZNm6S9vV3b0H7sQb8vbTabumZ4eFhee+01iYmJkfDwcHnppZfc/qM1m+lERHx4YoKIiIj8FF99QERERABYCoiIiGgCSwEREREBYCkgIiKiCSwFREREBIClgIiIiCawFBAREREAlgIiIiKawFJAREREAFgKiIiIaAJLAREREQEA/gdkHRSMaBjDugAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 1000x400 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "optic = Singlet()\n",
    "optic.draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "tolerancing = Tolerancing(optic)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Define perturbations:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "for k in range(1, 3):\n",
    "    # X-tilt\n",
    "    sampler = DistributionSampler(\"normal\", loc=0, scale=0.01)\n",
    "    tolerancing.add_perturbation(\"tilt\", sampler, surface_number=k, axis=\"x\")\n",
    "\n",
    "    # Y-tilt\n",
    "    sampler = DistributionSampler(\"normal\", loc=0, scale=0.01)\n",
    "    tolerancing.add_perturbation(\"tilt\", sampler, surface_number=k, axis=\"y\")\n",
    "\n",
    "    # X-decenter\n",
    "    sampler = DistributionSampler(\"normal\", loc=0, scale=0.1)\n",
    "    tolerancing.add_perturbation(\"decenter\", sampler, surface_number=k, axis=\"x\")\n",
    "\n",
    "    # Y-decenter\n",
    "    sampler = DistributionSampler(\"normal\", loc=0, scale=0.1)\n",
    "    tolerancing.add_perturbation(\"decenter\", sampler, surface_number=k, axis=\"y\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Define operands:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "input_data = {\n",
    "    \"optic\": optic,\n",
    "    \"surface_number\": -1,\n",
    "    \"Hx\": 0,\n",
    "    \"Hy\": 0.0,\n",
    "    \"wavelength\": 0.55,\n",
    "    \"num_rays\": 5,\n",
    "}\n",
    "tolerancing.add_operand(\"rms_spot_size\", input_data, target=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Run Monte Carlo without compensator:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "monte_carlo = MonteCarlo(tolerancing)\n",
    "monte_carlo.run(num_iterations=100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa8AAAGGCAYAAADMyA7/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/TGe4hAAAACXBIWXMAAA9hAAAPYQGoP6dpAABJPUlEQVR4nO3deXxTZb4/8E/SbG2zdU1a6MZadpC1wuiozDDqHUFw/bmg4zIL4ACDo9x7FXVUHL0qOoN49SoMow6Kg/uGoqKyL1M2oWylLW2TrkmapE3b5Pn9UZppoUBbkpwsn/frldfLnnNy8s0B++F5znOeRyaEECAiIoogcqkLICIi6imGFxERRRyGFxERRRyGFxERRRyGFxERRRyGFxERRRyGFxERRRyGFxERRRyF1AUEm8/nQ0VFBXQ6HWQymdTlEBHRWQgh0NDQgMzMTMjl525bRX14VVRUICsrS+oyiIiom8rKytC3b99zHhP14aXT6QC0XQy9Xi9xNUREdDYOhwNZWVn+39vnInl4lZeX44EHHsBnn30Gt9uNAQMGYOXKlRg3bhyAtmbkkiVL8Oqrr8Jms2Hy5MlYsWIFBg4c2K3zt3cV6vV6hhcRUQTozi0eSQds1NfXY/LkyVAqlfjss8/w448/4tlnn0VSUpL/mKeffhovvvgiXn75ZWzbtg2JiYmYNm0ampqaJKyciIikJJNyVvkHH3wQmzZtwvfff9/lfiEEMjMz8Yc//AGLFi0CANjtdphMJqxatQo33XTTeT/D4XDAYDDAbrez5UVEFMZ68vta0pbXhx9+iHHjxuH6669Heno6xowZg1dffdW/v7i4GBaLBVOnTvVvMxgMmDhxIrZs2dLlOT0eDxwOR6cXERFFF0nD6/jx4/77V1988QV++9vf4r777sPf/vY3AIDFYgEAmEymTu8zmUz+fadbunQpDAaD/8WRhkRE0UfS8PL5fLjooovw5JNPYsyYMbj33ntxzz334OWXX+71ORcvXgy73e5/lZWVBbBiIiIKB5KGV0ZGBoYOHdpp25AhQ1BaWgoAMJvNAACr1drpGKvV6t93OrVa7R9ZyBGGRETRSdLwmjx5MoqKijptO3z4MHJycgAAeXl5MJvN2LBhg3+/w+HAtm3bUFBQENJaiYgofEj6nNeCBQtw8cUX48knn8QNN9yA7du345VXXsErr7wCoG2s//z58/H4449j4MCByMvLw0MPPYTMzEzMmDFDytKJiEhCkobX+PHj8d5772Hx4sV47LHHkJeXh2XLluGWW27xH/PHP/4RLpcL9957L2w2G6ZMmYLPP/8cGo1GwsqJiEhKkj7nFQp8zouIKDJEzHNeREREvcHwIiKiiCP5xLwUvkpsNnx5/Ci8woe7xoyD4jzr6xARhQrDi7q0o+Ik/t8/30GcTI5W4UOhxYIXfnE1VHFxUpdGRMRuQzpTi9eL//r6K+QYjXhu2pX47dgJ2FB8DEt/2Ch1aUREABhe1IXXC3fheH0dbh0xGuo4BUabM/AfAwfjH/v3otrlkro8IiKGF3XW6vPhlV07cUl2LrL0Bv/2y3LzECeT4bXCXRJWR0TUhuFFnWw5WYr6pkZMycrptD1BqcKlOXl4c28h7FwIlIgkxvCiTj49chjpiYnINhjO2HdFXn80trbi82NHJKiMiOjfGF7k1+L14vOjhzE2ow9kMtkZ+40aDQYkp2A9w4uIJMbwIr/NZaWwezwYn9nnrMeMMZnxQ2kJGjyeEFZGRNQZw4v8NpYUIy0hEX11Z59TbIw5Ey0+HzaWFIewMiKizhhe5Le9ohwDkpO77DJsl5KQgByDEZ8fZdchEUmH4UUAAGdzMw7VVGNgcsp5jx1tMuO7khNo9flCUBkR0ZkYXgQAKLRUwicEBiSdP7zyU9PgbGnGj9VVIaiMiOhMDC8CAOysKIdOpYJZqz3vsbnGJKji4rC9/GQIKiMiOhPDiwAAOytOol/Sue93tVPI5eiXlIyt5WUhqIyI6EwML0Krz4d/WSq7db+r3eDkFOwoL4eX972ISAIML8KR2ho0traiX1Jyt98zKCUVDc0eFNXWBLEyIqKuMbwIh2raAuhcz3edLs+YBKVcjm2870VEEmB4EQ7VViMtIRHxSmW336OMi0OeMQm7K8uDWBkRUdcYXoSDNdXoo9P1+H05RiMKLZYgVEREdG4ML8Khmmr01Z85i/z55BmTUN7gQK3bHYSqiIjOjuEV42rcbtS43T2639Uu15AEANhbxdYXEYUWwyvGFdVWAwD69KLllZqQAJ1Khb1WhhcRhRbDK8YdqqmBKi4O6YmJPX6vTCZDjiGJ4UVEIcfwinGHaqrRR6eHvBsza3Ql99SgDSFEgCsjIjo7hleMaw+v3so1GlHf1IjyBkcAqyIiOjeGVwwTQuCErb5bk/GeTc6pQRv7qzjDPBGFDsMrhtW43XC1tMCU2PvwMqjVMKjVOFjD8CKi0GF4xbBiWz0AwHQBLS+ZTIY+OgMOVlcHqiwiovNieMWw4vo6yGUypCX0fKRhR331ei5MSUQhxfCKYcW2eqQmJEAhv7C/Bll6AyqcDXB4mgJUGRHRuTG8Ytjx+nqkJ/S+y7Bd1qkHnNtnpyciCjaGVww7bquDqRcPJ5/OpNVCKZez65CIQobhFaO8Ph9K7fYLGqzRTiGXI1Onx8EaDtogotBgeMWo8gYHWn2+Cxom31FfvR4H2PIiohBheMWo4vpTw+QDFV46A47W1cLr8wXkfERE58LwilHFtnoo5XIkxccH5Hx99Ho0e70osdsCcj4ionNheMWoUocdaQmJvZ6Q93TtKzEfrq0NyPmIiM6F4RWjTtrtSE5ICNj5dCo1dCo1DtdyuDwRBR/DK0aV2m1IjQ9ceLVNE6VjeBFRSEgaXo888ghkMlmnV35+vn9/U1MT5syZg5SUFGi1WsyaNQtWq1XCiqODEAInGxxIDWDLCwAydDocquVweSIKPslbXsOGDUNlZaX/9cMPP/j3LViwAB999BHWrl2LjRs3oqKiAjNnzpSw2uhQ39QId0tLQFteANBHp0eJzQZPa2tAz0tEdDqF5AUoFDCbzWdst9vteO211/DWW2/h8ssvBwCsXLkSQ4YMwdatWzFp0qRQlxo1TjraFo5MvcAJeU+XqdPDKwSO2+oxJDUtoOcmIupI8pbXkSNHkJmZiX79+uGWW25BaWkpAGDXrl1oaWnB1KlT/cfm5+cjOzsbW7ZskarcqHDSYQeAgHcbZmrbRxzyvhcRBZekLa+JEydi1apVGDx4MCorK/Hoo4/iJz/5Cfbv3w+LxQKVSgWj0djpPSaTCRaL5azn9Hg88Hg8/p8dDi5Pf7oyhx3xCiUSlMqAnjdRpUKyJp7hRURBJ2l4XXnllf7/HjlyJCZOnIicnBy88847iO/lw7NLly7Fo48+GqgSo1KZw4G0hATIAvSMV0cZOh2O8FkvIgoyybsNOzIajRg0aBCOHj0Ks9mM5uZm2Gy2TsdYrdYu75G1W7x4Mex2u/9VVlYW5KojT5ndhuQAzaxxugytDkfqGF5EFFxhFV5OpxPHjh1DRkYGxo4dC6VSiQ0bNvj3FxUVobS0FAUFBWc9h1qthl6v7/Sizk46HAEfrNEuQ6tDmcPOEYdEFFSShteiRYuwceNGnDhxAps3b8a1116LuLg43HzzzTAYDLjrrruwcOFCfPPNN9i1axfuvPNOFBQUcKThBfAJgfIgPOPVLkOng08IFNvqg3J+IiJA4nteJ0+exM0334za2lqkpaVhypQp2Lp1K9LS2oZZP//885DL5Zg1axY8Hg+mTZuGl156ScqSI161y4Vmrzd44XVqxOHRulrkc7g8EQWJpOG1Zs2ac+7XaDRYvnw5li9fHqKKol9FQ9voy5QAP6DcTqtSwaDW8L4XEQVVWN3zouCrdDYAAJI0wRmwAQAZOi2O1dUF7fxERAyvGFPR0AB1nCLgz3h1lJHICXqJKLgYXjGmwtmAlIT4oDzj1S5Dp8MJuw0tXm/QPoOIYhvDK8ZUNjQEtcsQaBu00erzoZSrKhNRkDC8Yky5wx6S8AKAo/W870VEwcHwijGVzoagza7RTq9WI1GpxFEO2iCiIGF4xRBPaytqGxuDHl4ymQxmrQ7HOFyeiIKE4RVDLE4ngOAOk29n1mpxlOFFREHC8Ioh7c94JQfpAeWOMrQ6HLfVQwgR9M8iotjD8IohlQ2nHlCO1wT9s8xaHdwtLf7AJCIKJIZXDKlwNkCnUkEdF/xZwdpHHHKmDSIKBoZXDKlocISkyxAAUhMSoJTLcYzD5YkoCBheMaSyoQFGTfC7DAFAfmrEIZ/1IqJgYHjFkIoGB5JDMNKwnSlRi6O1HHFIRIHH8IohVpcTxiA/49VRhlaLo/UMLyIKPIZXjGhqbYHd4wlZtyHQNkFvXWMjbE2NIftMIooNDK8YYXW6AITmAeV25vYRh7zvRUQBxvCKEVZX2+waoWx5mRK1kIHD5Yko8BheMcLiX0E5dOGliotDakIiW15EFHAMrxhhdTkRr1BAowjeCspdMWu1bHkRUcAxvGKExekM6f2udmatFkc4QS8RBRjDK0ZYnU4YQthl2C5Dq0N5gwOe1taQfzYRRS+GV4yodIZudo2OzFodfEKg2FYf8s8moujF8IoRVqcTRgm6Ddsn6D3OQRtEFEAMrxjgEwJVbldIRxq206pU0KnUOMpBG0QUQAyvGFDX2IhWn0+SbkOA00QRUeAxvGKA9dQzXlJ0GwJt9704QS8RBRLDKwZYTs2uIUW3IdDW8iq21cMnhCSfT0TRh+EVA6xOJ+QyGfRqacLLrNXB4/Wi3OGQ5POJKPowvGKA1eWEUaOBXCaT5PPbRxzyvhcRBQrDKwZUuVwwSNTqAoCk+Hio4+I4TRQRBQzDKwZUuZzQq9WSfb5cJmsbtMFpoogoQBheMcDidEo2TL4d5zgkokBieMWAKpdT0m5DoO2+17H6OgiOOCSiAGB4RblWnw91jY2St7wytDo4PB7UNLolrYOIogPDK8rVuF0QQFi0vADwYWUiCgiGV5SrcrkAQPKWV1piIuJkMhzlBL1EFAAMryhXdWp2DalbXgq5HGatDsc4aIOIAoDhFeWqXC7IZTLoJBwq344jDokoUBheUc56aqShVLNrdJSh1TG8iCggGF5Rrsrlkvx+V7tMnQ41bjdsTY1Sl0JEEY7hFeWknl2jI/+IQ04TRUQXiOEV5SxO6R9QbmdK1EIuk3GaKCK6YGETXk899RRkMhnmz5/v39bU1IQ5c+YgJSUFWq0Ws2bNgtVqla7ICFTtDp9uQ2VcHNITE9nyIqILFhbhtWPHDvzv//4vRo4c2Wn7ggUL8NFHH2Ht2rXYuHEjKioqMHPmTImqjDytPh9q3e6waXkBbWt7HamrkboMIopwkoeX0+nELbfcgldffRVJSUn+7Xa7Ha+99hqee+45XH755Rg7dixWrlyJzZs3Y+vWrRJWHDn8s2uEScsLaLvvdZizbBDRBZI8vObMmYOrr74aU6dO7bR9165daGlp6bQ9Pz8f2dnZ2LJly1nP5/F44HA4Or1iVbW7bR7BcOk2BIBMrQ5WlxMOj0fqUogogkkaXmvWrMHu3buxdOnSM/ZZLBaoVCoYjcZO200mEywWy1nPuXTpUhgMBv8rKysr0GVHjH/PrhEeow2BtuHyADhog4guiGThVVZWht///vd48803oQlgy2Dx4sWw2+3+V1lZWcDOHWmqXS7IAOhU4RNeGVod5DIZjtTyvhcR9Z5k4bVr1y5UVVXhoosugkKhgEKhwMaNG/Hiiy9CoVDAZDKhubkZNput0/usVivMZvNZz6tWq6HX6zu9YlW12wW9Wo04ueS9w37tIw4Ps+VFRBdAIdUHX3HFFdi3b1+nbXfeeSfy8/PxwAMPICsrC0qlEhs2bMCsWbMAAEVFRSgtLUVBQYEUJUecKpcrrEYatuOgDSK6UJKFl06nw/DhwzttS0xMREpKin/7XXfdhYULFyI5ORl6vR7z5s1DQUEBJk2aJEXJEafa5YI+DMMrU6fD9vKTUpdBRBFMsvDqjueffx5yuRyzZs2Cx+PBtGnT8NJLL0ldVsSwhtHUUB1lavWodrthb2oKq2H8RBQ5wiq8vv32204/azQaLF++HMuXL5emoAhX7XIh0xx+9/zaRxwerqvB+My+EldDRJEofO7kU0AJIVDtdoXVMPl27XMcHuF9LyLqJYZXlHJ4PGjx+cJywIYyLg5mrRaHOVyeiHqJ4RWlqlwuAOE1NVRHGVodDtUwvIiodxheUarafSq8wrDlBQB9dXoU1VZDCCF1KUQUgRheUcrf8grDe14A0Eevh93j8ddJRNQTDK8oVe12Il6hgFoRVgNK/froDACAIt73IqJeYHhFqWqXK2zvdwFAakICNHEKFNVWS10KEUUghleUqnKH5+wa7eQyGTJ1OhRx0AYR9QLDK0pVOcPzGa+OMnV6HKxhy4uIeo7hFaWq3OE5NVRHfXR6HKuvQ6vPJ3UpRBRhGF5RqsbtDtth8u366PVo9npRYquXuhQiijAMryjkaW2Fw+MJ+/Dqe2qtNT6sTEQ9xfCKQjVuNwCEfbehTqVGkiae972IqMcYXlHIP7tGGA+Vb9dHr8fBmiqpyyCiCMPwikLVrvCeGqqjLJ0eP1az5UVEPcPwikLVbhfkMhm0KpXUpZxXX70BVpcT9Y2NUpdCRBGE4RWFqlxtz3jJZTKpSzmvLH3bNFG870VEPcHwikLVYT67RkfpiYlQyeMYXkTUIwyvKFTtcoX9SMN2cXI5+uj1OMTwIqIeYHhFoSpX+M+u0VFfnR4/VnPEIRF1X6/C6/jx44GugwKoOgJm1+goy2DAkbpaNHu9UpdCRBGiV+E1YMAAXHbZZXjjjTfQ1NQU6JroAgghUON2RVZ46Q1o9flwhGt7EVE39Sq8du/ejZEjR2LhwoUwm8349a9/je3btwe6NuoFh8eDFp8v7GeU76iv3gAZgAPsOiSibupVeI0ePRovvPACKioq8Prrr6OyshJTpkzB8OHD8dxzz6GaD51KpurUA8r6CJhdo51GoUCGVsf7XkTUbRc0YEOhUGDmzJlYu3Yt/vznP+Po0aNYtGgRsrKycPvtt6OysjJQdVI3+aeGiqCWF9DW+tpXZZW6DCKKEBcUXjt37sTvfvc7ZGRk4LnnnsOiRYtw7NgxfPnll6ioqMD06dMDVSd1U3t4RcpzXu2yDQYcqqmGl2t7EVE3KHrzpueeew4rV65EUVERrrrqKqxevRpXXXUV5PK2LMzLy8OqVauQm5sbyFqpG6pdLqjjFNAoevVHK5lsgxGNra04YatH/+QUqcshojDXq99wK1aswK9+9SvccccdyMjI6PKY9PR0vPbaaxdUHPVctdsFYwTd72qXfWqaqAPVVQwvIjqvXoXXl19+iezsbH9Lq50QAmVlZcjOzoZKpcLs2bMDUiR1XyTNrtFRokqF1IQEHKiuwjWDh0hdDhGFuV7d8+rfvz9qulj9tq6uDnl5eRdcFPVe27yGkRdeQFvra5+VgzaI6Px6FV5CiC63O51OaCKwyyqaWJ3OiHpAuaMcYxL2V1nhO8vfLyKidj3qNly4cCEAQCaT4eGHH0ZCQoJ/n9frxbZt2zB69OiAFkg9U+N2Y3i6SeoyeiXXYISzpRknbPXol5QsdTlEFMZ6FF7/+te/ALS1vPbt2wdVh8UOVSoVRo0ahUWLFgW2Quq2Zq8XNk9TRA7YAIAcgxEAsK/KyvAionPqUXh98803AIA777wTL7zwAvR6fVCKot6pdbsBRN4zXu0SVSqkJyZin9WK6Ry0QUTn0KvRhitXrgx0HRQAVRE6u0ZH2Xoj9lZZpC6DiMJct8Nr5syZWLVqFfR6PWbOnHnOY9etW3fBhVHP1bjawysyW14AkGs04pMjRfD6fIiTc7k5Iupat8PLYDBAJpP5/5vCT7XbBblMBl0Et7xyTs20cay+DoNSUqUuh4jCVLfDq2NXIbsNw1OVywWdSg35qX9kRKIcgxEyAHutFoYXEZ1Vr/plGhsb4T41OAAASkpKsGzZMqxfvz5ghVHPVbtdEX2/CwDilUpk6vTYY+V9LyI6u16F1/Tp07F69WoAgM1mw4QJE/Dss89i+vTpWLFiRUALpO5rm10jcu93tcsxGPEvC5fTIaKz6/VKyj/5yU8AAO+++y7MZjNKSkqwevVqvPjiiwEtkLrP6nRG7NRQHeUlJaGopgZNrS1Sl0JEYapX4eV2u6HT6QAA69evx8yZMyGXyzFp0iSUlJQEtEDqvmq3O+K7DQEgz5gEr/DhAFdWJqKz6FV4DRgwAO+//z7KysrwxRdf4Oc//zkAoKqqqkcPLq9YsQIjR46EXq+HXq9HQUEBPvvsM//+pqYmzJkzBykpKdBqtZg1axasnLi1S0II1LpdMETo7Bod9dHpoZLHYY+F972IqGu9Cq+HH34YixYtQm5uLiZOnIiCggIAba2wMWPGdPs8ffv2xVNPPYVdu3Zh586duPzyyzF9+nQcOHAAALBgwQJ89NFHWLt2LTZu3IiKiorzPmMWqxqaPfB4vVFxz0shlyPbYEChlfe9iKhrMnG2KeLPw2KxoLKyEqNGjfKv67V9+3bo9Xrk5+f3uqDk5GQ888wzuO6665CWloa33noL1113HQDg0KFDGDJkCLZs2YJJkyZ163wOhwMGgwF2uz2qp7M6VleLn72xCvcXTImKIeZvH9iHgzXV+P7Oe6QuhYhCpCe/r3s9hYHZbMaYMWM6LUg5YcKEXgeX1+vFmjVr4HK5UFBQgF27dqGlpQVTp071H5Ofn4/s7Gxs2bKlt2VHraoomF2jo35JyShvcKD61JRXREQd9WpuQ5fLhaeeegobNmxAVVUVfD5fp/3Hjx/v9rn27duHgoICNDU1QavV4r333sPQoUNRWFgIlUoFo9HY6XiTyQTLOe6FeDweeDwe/88Oh6PbtUQy/7yGUXDPCwD6JyUBAP5VWYGf9x8ocTVEFG56FV533303Nm7ciNtuuw0ZGRn+aaN6Y/DgwSgsLITdbse7776L2bNnY+PGjb0+39KlS/Hoo4/2+v2RqtrlgjpOAY2iV3+kYSc5PgHJ8fHYbalkeBHRGXr1m+6zzz7DJ598gsmTJ19wASqVCgMGDAAAjB07Fjt27MALL7yAG2+8Ec3NzbDZbJ1aX1arFWaz+aznW7x4sX/RTKCt5ZWVlXXBdYa7arcrYtfxOps8YxJ2VZRLXQYRhaFe3fNKSkpCcnJwFgv0+XzweDwYO3YslEolNmzY4N9XVFSE0tJS/+jGrqjVav/Q+/ZXLKhyRf7UUKfrn5SM/VVWNHu9UpdCRGGmVy2vP/3pT3j44Yfxt7/9DQkJCb3+8MWLF+PKK69EdnY2Ghoa8NZbb+Hbb7/FF198AYPBgLvuugsLFy5EcnIy9Ho95s2bh4KCgm6PNIwlVS5nVAyT76h/UjI8Xi8O1lRjlOnsrW0iij29Cq9nn30Wx44dg8lkQm5uLpRKZaf9u3fv7tZ5qqqqcPvtt6OyshIGgwEjR47EF198gZ/97GcAgOeffx5yuRyzZs2Cx+PBtGnT8NJLL/Wm5KhndTmRZ0ySuoyAyjYYoZTLsbuyguFFRJ30KrxmzJgRkA9/7bXXzrlfo9Fg+fLlWL58eUA+L5rVuN0YbcqQuoyAUsjlyD113+vO0RdJXQ4RhZFehdeSJUsCXQddAE9rKxweT9Q849VR/6Rk7KgohxDigka1ElF06fVDyjabDf/3f/+HxYsXo66uDkBbd2F5OUeHhVp1lD3j1dGA5BRUu104GSPP6xFR9/Sq5bV3715MnToVBoMBJ06cwD333IPk5GSsW7cOpaWl/rW+KDT+PbtGdI02BIABSW2jWndWlCPLYJC4GiIKF71qeS1cuBB33HEHjhw5Ak2Hf+1fddVV+O677wJWHHWPv+UVhd2GiSoV+uj02FnJFj0R/VuvwmvHjh349a9/fcb2Pn36nHPqJgqOKpcLcTIZElUqqUsJiv5JydheflLqMogojPQqvNRqdZdzBh4+fBhpaWkXXBT1TLWrbR0veZQOaBiYnIJj9XWob2yUuhQiChO9Cq9rrrkGjz32GFpa2pZpl8lkKC0txQMPPIBZs2YFtEA6vyqXMyq7DNsNSP73fS8iIqCX4fXss8/C6XQiLS0NjY2NuPTSSzFgwADodDo88cQTga6RzqPK5YI+CgdrtEtNSERKfAK2seuQiE7p1WhDg8GAL7/8Eps2bcKePXvgdDpx0UUXdVp7i0Knyu1Canzvp+mKBAOTU7C1vEzqMogoTPQ4vHw+H1atWoV169bhxIkTkMlkyMvLg9ls5oOkEqlyOdE/KTgTJYeLQSkpeGPfHjg8nqhuZRJR9/So21AIgWuuuQZ33303ysvLMWLECAwbNgwlJSW44447cO211warTjqLVp8PtW43jFF8zwsABiWnwicEdnHIPBGhhy2vVatW4bvvvsOGDRtw2WWXddr39ddfY8aMGVi9ejVuv/32gBZJZ1frdkMgOmfX6Cg9MRFGjQbbyk/istx+UpdDRBLrUcvrH//4B/7zP//zjOACgMsvvxwPPvgg3nzzzYAVR+dndTkBIOoWojydTCbDwOQUbDvJ+15E1MPw2rt3L37xi1+cdf+VV16JPXv2XHBR1H3Vp6aGivZuQwAYnJKK/VVWOJubpS6FiCTWo/Cqq6uDyWQ6636TyYT6+voLLoq6z+pyQi6TQRcDgxgGp6TBKwR2VHDIPFGs61F4eb1eKBRnv00WFxeH1tbWCy6Kus966gHlaJ1doyNTYiKSNfHYyq5DopjXowEbQgjccccdUJ/lX/kejycgRVH3VbtcUX+/q51MJsOglBRsKiuVuhQikliPwmv27NnnPYYjDUPL6nLG1HNP+alp+Nuef8He1BT1IyyJ6Ox6FF4rV64MVh3USxanE2kJiVKXETKDU1IhAGwvP4mf9R8gdTlEJJFer6RM4aHa7YqpFkhqQiLSEhKx+SS7DoliGcMrgsXK7Bqny09NxfelJ6Qug4gkxPCKYDVuFwSi/wHl0w1JTcfx+npYnA1Sl0JEEmF4RbCqUw8ox1K3IdDW8pIBHHVIFMMYXhGsqn1qqBjrNtSp1Mg2GLGptETqUohIIgyvCGZ1uWJmdo3T5aem4YeyEgghpC6FiCTA8IpgVS4njJrYmF3jdENT01DjduNwXa3UpRCRBBheEczidMZcl2G7AckpUMnj8H3JCalLISIJMLwimMXZEHODNdqp4uIwKCUVG0uKpS6FiCTA8IpgFqcz5obJdzQsLQ07KsrR2NIidSlEFGIMrwhmdTmRpImXugzJDEszodnrxbZyLpFCFGsYXhGqqbUFDo8nplteZq0WKfHx+I6zbRDFHIZXhLI4Tz3jFcMtL5lMhqGp6dh4gve9iGINwytCWU+FV1IMt7wAYHi6CcW2epTabVKXQkQhxPCKUBYXW14AMCQ1DQq5HN+y9UUUUxheEcrqdCJeoYRG0aMl2aJOvFKJAckp+ObEcalLIaIQYnhFKIvLiaT42O4ybDci3YStJ8s4ZJ4ohjC8IpTV6YQhRmfXON2IdBM8Xi+2nCyTuhQiChGGV4SyOBtieph8R+ZELdITE/E1uw6JYgbDK0JZnLH9gHJHMpkMI9PN+Or4Uc4yTxQjGF4RyCcEqlwutrw6GGUyo8rlwoHqKqlLIaIQYHhFoNpGN7zCx5ZXBwOSU5CoVOKr48ekLoWIQoDhFYEs/geUGV7tFHI5hqWb8OXxo1KXQkQhIGl4LV26FOPHj4dOp0N6ejpmzJiBoqKiTsc0NTVhzpw5SElJgVarxaxZs2C1WiWqODxUNjgAgEPlTzPaZMbBmmqUn7o+RBS9JA2vjRs3Ys6cOdi6dSu+/PJLtLS04Oc//zlcLpf/mAULFuCjjz7C2rVrsXHjRlRUVGDmzJkSVi09i9MJhVwOrUotdSlhZViaCQq5HF+x9UUU9SSdnuHzzz/v9POqVauQnp6OXbt24ZJLLoHdbsdrr72Gt956C5dffjkAYOXKlRgyZAi2bt2KSZMmSVG25CqcDUiOj4dcJpO6lLCSoFQiPzUNnx89gtmjLpK6HCIKorC652W32wEAycnJAIBdu3ahpaUFU6dO9R+Tn5+P7OxsbNmyRZIaw0FlQwPvd53FReYM7KgoR63bLXUpRBREYRNePp8P8+fPx+TJkzF8+HAAgMVigUqlgtFo7HSsyWSCxWLp8jwejwcOh6PTK9pUNDg4TP4sRpkyAIBdh0RRLmzCa86cOdi/fz/WrFlzQedZunQpDAaD/5WVlRWgCsNHpdOJZLa8uqRXqzEwOQWfHT0idSlEFERhEV5z587Fxx9/jG+++QZ9+/b1bzebzWhubobNZut0vNVqhdls7vJcixcvht1u97/KyqJrvru2B5SdSIpneJ3NGHMGNp8shb2pSepSiChIJA0vIQTmzp2L9957D19//TXy8vI67R87diyUSiU2bNjg31ZUVITS0lIUFBR0eU61Wg29Xt/pFU1q3C60+viA8rmMzciE1+fDenYdEkUtSUcbzpkzB2+99RY++OAD6HQ6/30sg8GA+Ph4GAwG3HXXXVi4cCGSk5Oh1+sxb948FBQUxOxIw8pTDygns+V1VkZNPAalpOLjw4dw/dDhUpdDREEgaXitWLECAPDTn/600/aVK1fijjvuAAA8//zzkMvlmDVrFjweD6ZNm4aXXnopxJWGj8qGBgBgt+F5jM3IxNsH9qGu0Y3k+ASpyyGiAJM0vLozA7hGo8Hy5cuxfPnyEFQU/iqdDVDJ46BVqqQuJayNzcjEmgP78MWxo7h5+EipyyGiAAuLARvUfZUNDiTFx0PGB5TPSa/WID8lFR8cOih1KUQUBAyvCGNxOjmnYTdN7JOF7RUnOdchURRieEWY8gYHRxp20xhzBlRxcfiwiK0vomjD8Iow5Q0OjjTspnilEqNMZqw7+CNXWCaKMgyvCNLs9aLa5UIKR89126Q+WThWX8cVlomiDMMrglicDRAAw6sHhqWlw6jW4N0f90tdChEFEMMrgpSfmmSYzy11X5xcjkl9s/B+0UF4WlulLoeIAoThFUHaR82l8J5Xj0zOyobD4+F0UURRhOEVQU46HDCoNVDGxUldSkQxa3UYmJyCtw/sk7oUIgoQhlcEKW9wsNXVS1Oyc7C5rBQnbPVSl0JEAcDwiiDlDgfnNOyl8Rl9oFWp8Na+PVKXQkQBwPCKICcb7EjlYI1eUcbFYXLfbLzz4340tbZIXQ4RXSCGV4TwCQGL04mUBIZXb12Sk4cGjwcfHS6SuhQiukAMrwhR5XKi1efjMPkLkJ6YiBHpJry2eydn3CCKcAyvCMFh8oHxs34DcLiuFj+UlUhdChFdAIZXhGh/QJmza1yYwSmpyDEY8equnVKXQkQXgOEVIU46HEhUKhGvVEpdSkSTyWSYmtcfP5SV4CDnOySKWAyvCHHSYUdqQqLUZUSFcZl9kJaQiL/u2Cp1KUTUSwyvCFFityGVIw0DQiGX48oBA/H50SM4XFsjdTlE1AsMrwhRYrcjjS2vgCnom43k+AT8dTtbX0SRiOEVAVq8XlicDQyvAGpvfX1ypAg/8t4XUcRheEWAioYG+IRgt2GATc7KgUmrxdObvpe6FCLqIYZXBCh12ACALa8AU8jlmDF4CL4rPYEtZaVSl0NEPcDwigBldjvkMhkn5Q2Ci8yZ6J+UjMe++watPp/U5RBRNzG8IkCpw46U+AQo5PzjCjSZTIabho3A4doazjhPFEH42zAClNntvN8VRLnGJEzJzsGzWzah2u2Suhwi6gaGVwTgM17Bd23+UADAo99+LXElRNQdDK8IUMZnvIJOp1LjxmEj8OnRw/ji2BGpyyGi82B4hTl7UxMamj1seYXAhMw+GG0y47++/hI1brfU5RDROTC8wtwJWz0AID1RK3El0U8mk+G2kaPR6vPhj19+zjW/iMIYwyvMFZ8KL1Miuw1DQa/W4I6RY/BtSTFW7fmX1OUQ0VkwvMJcsa0eSRoNNAouhRIqI0xm/Kxffzz5/UbsrCiXuhwi6gLDK8wV2+rZZSiBmfnD0C8pCXM+/QhVLqfU5RDRaRheYe54fT3SOdIw5BRyOe69aDxafT7c89H7aGptkbokIuqA4RXGhBAosdXDpGXLSwpGjQZzxk1EUW0NFq3/HD4O4CAKGwyvMFbjdsPV0gITuw0lk2M04q7RY/HZ0cN4etN3UpdDRKcwvMJYsX+YPLsNpXRRRiZuHDYCr+zeiVWFu6Uuh4gAKKQugM6uuL4OMnAplHBwRV5/1DU24k/ffYOUhAT8clC+1CURxTSGVxgrttuQmpAIZVyc1KUQgFlDhsHR7MEf1n8Go1qDn+TkSl0SUcxit2EYK66vY5dhGJHLZJg9cgyGpKbhN598gEJLpdQlEcUshlcYO1pXx8EaYUYhl+PXY8ejj96AX324Dsfr66QuiSgmMbzClKe1FaV2GzJ1OqlLodOo4xSYO24iEpRK3Pbeu7A4G6QuiSjmMLzC1Am7DV4hkKnTS10KdSFRpcLvJ1wMj7cVv/pgHZzNzVKXRBRTJA2v7777Dr/85S+RmZkJmUyG999/v9N+IQQefvhhZGRkID4+HlOnTsWRI7Gx1tKR2hoAQKaWLa9wlRwfj3njJ6HEbsfcTz9Cq88ndUlEMUPS8HK5XBg1ahSWL1/e5f6nn34aL774Il5++WVs27YNiYmJmDZtGpqamkJcaegdqauFUaNBokoldSl0Dn31Bvxm7HhsKivB4999I3U5RDFD0qHyV155Ja688sou9wkhsGzZMvz3f/83pk+fDgBYvXo1TCYT3n//fdx0002hLDXkjtTWwsxWV0QYmpaOm4aPxOq9hRiYkopbRoySuiSiqBe297yKi4thsVgwdepU/zaDwYCJEydiy5YtZ32fx+OBw+Ho9IpEh+tq2GUYQX6ak4fLcvPwyLdfY1cll1EhCrawDS+LxQIAMJlMnbabTCb/vq4sXboUBoPB/8rKygpqncHQ7PWixGZHBkcaRpQbho5Av6Qk/O6Tj1DtckldDlFUC9vw6q3FixfDbrf7X2VlZVKX1GMnbPXwCh9bXhGmfRmVFq8X87/4BF4O4CAKmrANL7PZDACwWq2dtlutVv++rqjVauj1+k6vSHOkthYAOEw+Ahk1GvxqzFhsPVmGl3dtl7ocoqgVtuGVl5cHs9mMDRs2+Lc5HA5s27YNBQUFElYWfEW1NTBqNNBypGFEGpKahqsHDsbzWzdjd2WF1OUQRSVJw8vpdKKwsBCFhYUA2gZpFBYWorS0FDKZDPPnz8fjjz+ODz/8EPv27cPtt9+OzMxMzJgxQ8qyg+5AtRV9dQapy6AL8B8DByPPmISF6z+Fiw8wEwWcpOG1c+dOjBkzBmPGjAEALFy4EGPGjMHDDz8MAPjjH/+IefPm4d5778X48ePhdDrx+eefQ6PRSFl20B2orkKWgeEVyeLkctw56iJYnU48+cNGqcshijoyIaJ7bXOHwwGDwQC73R4R979q3W6M/78VuPei8Rif2UfqcugCfXPiON7avxd/v/Y6TM7KkbocorDWk9/XYXvPK1b9WFMFAMjWs+UVDS7NyUN+SioWf7We3YdEAcTwCjMHq6uhjlMgjet4RQW5TIbbRo5BlduF57ZukrocoqjB8AozP9ZUIUuvh1wmk7oUCpD0xERcMygff9vzL+yvsp7/DUR0XgyvMLO/yoosdhlGnSvy+qOPTo/FG9bz4WWiAGB4hZHGlhacsNk40jAKKeRy3DJiFH6srsKb+/ZIXQ5RxGN4hZED1VXwCYEcg1HqUigI+iclY0p2Dp7dsgk1brfU5RBFNIZXGCm0VEIdF4c+nBYqal2bPxQCAn/e9J3UpRBFNIZXGPmXpQI5xiTEyfnHEq10KjWmDx6Cfx48gEJLpdTlEEUs/pYMI7srK5FnTJK6DAqyS7JzkaU34JFvN8AX3XMEEAUNwytMWJwNsLqc6MfwinpymQw3DRuBvVVWrDt4QOpyiCISwytMFJ5aYLNfEsMrFgxKScX4zD54evP3nHmDqBcYXmGi0FKBlPgEGDXxUpdCITIrfxgcHg/X/SLqBYZXmNhZWYFco1HqMiiEUhIS8LN+/fHq7p046bBLXQ5RRGF4hYHGlhbstVowOCVV6lIoxK7sPwiJSiX+vOl7qUshiigMrzCws7IcrT4fwysGqRUKzMgfik+OFGFnRbnU5RBFDIZXGNhSVgaDWoMMrU7qUkgCk/pkIddoxKMbv+bQeaJuYniFgS0nSzEoJQUyziQfk+QyGW4cOgIHqqvw3qEfpS6HKCIwvCTW4PFgf5WVXYYxbkByStvQ+U0cOk/UHQwvie2oKIdXCOSnpEldCklsVv4w2D1NeGnnNqlLIQp7DC+JfV96Ainx8UjnyskxLyUhAT/vNwD/t3snSmw2qcshCmsMLwkJIfDV8WMYkW7m/S4CAPxiwEDo1Wo88f23UpdCFNYYXhI6XFeL8gYHRpnMUpdCYUIdp8CsIcPwVfExbDxRLHU5RGGL4SWhDcePQROn4GAN6mRcRh/kp6ZiycYN8LS2Sl0OUVhieEloQ/ExDElLgzIuTupSKIzIZDLcPGwkyh0OvLp7p9TlEIUlhpdEql0uFFoq2WVIXcrU6TG13wD8dcdWnLDVS10OUdhheEnk4yNFiJPLMdqUIXUpFKZ+OXAw9Go1/vubryA48wZRJwwviaw7eAAj0k1IVKmkLoXClFqhwP8bPgqby0o58wbRaRheEjhaV4sD1VWY1CdL6lIozI1IN2FCn7547LtvUOVySl0OUdhgeEngg6KDSFQqMSLdJHUpFAFuHjYSAPCfX3/J7kOiUxheIdbi9WLtgf0Ym9GHowypW7QqFW4ZPgpfFx/H2h/3S10OUVhgeIXYZ0cPo8rtwmW5eVKXQhHkooxMTM7KxqMbv0YxRx8SMbxCbWXhbgxJTUNfvUHqUijC3DRsJPRqNe777GM+vEwxj+EVQv+qrMAeqwVX5PWXuhSKQBqFAvdcNB6Ha2vw2HffSF0OkaQYXiH0/NbNyNDqOFCDei3HYMTNw0fiH/v34p0D+6Quh0gyDK8Q2VxWih/KSjBj8BDIOYM8XYCfZOfikuxc/Pc3X2HryTKpyyGSBMMrBIQQeHrTd+hnTMIYM2fUoAt38/CRGJScgt988gGKamukLoco5BheIfCP/Xuxt8qKmflDuW4XBYRCLsevx46HQa3Bbe+t5fyHFHMYXkFW7nDgyR824idZORicmiZ1ORRFEpQqzJ94MZRyOW7+5zs4Xl8ndUlEIaOQuoBo1uL1YtGXn0GjUOC6ocOlLoeikF6txsJJk7Fs2xZcv3YNVk6fiZFRsFKB1+dDjduN2kY36hob4fB40NjagqbWVnh9PviEgFwmgyouDhqFAlqVCnq1BkmaeKQlJsCg1rCXI8oxvILo8e+/xc6KciyYNBkJSqXU5VCUMmrisahgCv66YytufHcNll4xDTPyh0hdVrdUu104VF2NotoaHK+vQ7GtHqV2O6pcTni7mApLLpMhTiaDDDL4INDq83V5XlVcHDK0OvTVG5BrNCLPmIT+SckYmJKCDK2OwRYFZCLKJ0tzOBwwGAyw2+3Q6/Uh+9yXd27H05u/x20jRuOSnNyQfS7FrhavF3/fV4gtJ8tww9DheOiSy8Jq1QJXczP2Wi3YbalAoaUSe61WVLtdAP4dNmkJiUhLSEByQgKSNBro1RpolSokKJVQKxRQyDvf6RBCoMXnQ1NrC9wtLWhoboajqQn1TU2oa3Kjxt32sjqdaPZ5AbRNt5WfmoahqWkYkpaOYWnpGJicArWC/5aXWk9+XzO8AkwIgee3bsZfd2zFfwwcjOmDI+NfwBQdhBD4oawEbx/Yj5SEeDx0yWX4eb8BkrQ0bE2N2FlRjm3lJ7Gt/CQOVlfBKwQSlErkGpOQYzAix2BAlt6A1ITEoD5C4hMCtW43KpwOnHQ4UN7gQLnDgUpnAwSAOJkc/ZKSMDQtHYNTUjEwJQUDklLQR68/IzB7o6m1BRanE1anE1aXE1UuF2oa3ahvbIS9qQnOZg/cra1o8bYFrFwmg0ahQKJSBaNGg5SEBJgSteij1yNbb0BeUnJU9uZEXXgtX74czzzzDCwWC0aNGoW//OUvmDBhQrfeG8rwqm9sxKIvP8c3J47juiHDMK3/wKB+HtHZVLlcWHNgL/ZVWXFRRiZ+P6EAU7Jzghpi1W4XdpSXY0fFSWw9WYbDtTUQAFLi4zEgOQUDk1MwILmt2y5cnnX0tLbiZIMDZQ47yuz2tlBrcKDp1PRbCrkcGVodsvR6mLQ6pMTHw6CJR6JSiXiFAnK5HEK0dV82trbC1dwMu6cJ9Y2NqHa7UeVyosrlhN3j6fS56jgFjBo1tCo14k+dSymP8wdlW4vSi8bW1lMtSg/qGxvhORVuAJCh1WFwSiqGpKVhaGo6hqebkG0wRHSXaFSF19tvv43bb78dL7/8MiZOnIhly5Zh7dq1KCoqQnp6+nnfH4rwavX5sPbH/fifzT+gxefFHaMuwqgouGlOkU0Igf3VVfj48CEct9Ujx2DE9UOH4xcDBqJfUvIFnbuptQWHamqw12rBHqsFOyvKUeawAwDSExMxMDkFg5JTMSglBakJiYH4OiEjhEBdU6O/ldQ+cMTh8aCh2QNXcws8ra3+bkgAkKFt8VCNQoEEpRKJShV0KjUMGjWMmngYNW2DSZI0Ghg1GmgUPW81CSHgbGlGtcsFi8uJyoYGfwuyrqkRAKBTqTE8PR0jTGaMSDNFXKBFVXhNnDgR48ePx1//+lcAgM/nQ1ZWFubNm4cHH3zwvO8PZniVOxz48PBBvLF3DyqdDSjom4WZ+UNh1MQH9HOILoQQAofravFd6QnssVTC4/UiQ6vD+D59kJ+Shn5JSTBrdTBqNEhQqiCXAV4h0NjSAltTE+oaG1HpbECp3YYTNhsO19ag1G6DVwgo5HJk6Q3oZ0xCv6RkDExOQVJ8bPz9F0LAJwRkMpnkLckGjwcldhtK7DacsNtQZrejttENoO0e39DUdAxNS8Pg1DQMSk5B/+QU6NVqSWvuStSEV3NzMxISEvDuu+9ixowZ/u2zZ8+GzWbDBx98cN5zBCq8HJ4mlNrtOFxbg31VVmw9WYai2hoo5XKMz+yLqf36I4szxVOY83hbcaimfXRfPSoaHGjsxgz1cpkMyZp4pCUmIkOrQx+dHlkGA/rq9FyXLkw5TgVaqd2GUrsdFU4HrE4n2n/hp8QnIMdoRK7BiD56PTK0Opi0WqQnJCI1IRFGjSbkg1h68vs6rIfX1NTUwOv1wmTqPJGtyWTCoUOHunyPx+OBp0P/st3e1pXhcDh6XcefNn6NtQcPnLE9Q6vD8HQT1DI5dhYXY2evP4EotJQABmv1GJSog6ulGbamRrhbWtHsbYVAW1gp5XFIUCqgVamhV6khbx+44PWh1mZDrc2GQgm/A3WfWamCOSkVrYYkWJwNsLpcqKqvR3V9XUB/b8XJ5Hj0p5fjml4OVGv/Pd2dNlVYh1dvLF26FI8++ugZ27OysgL+WSUAtgb8rEREkeu2AJyjoaEBBsO5e7LCOrxSU1MRFxcHq9XaabvVaoXZ3PWAiMWLF2PhwoX+n30+H+rq6pCSkhLwm5YOhwNZWVkoKysL6TNkkYzXrOd4zXqO16znwuGaCSHQ0NCAzMzM8x4b1uGlUqkwduxYbNiwwX/Py+fzYcOGDZg7d26X71Gr1VCfdiPSaDQGtU69Xs//QXqI16zneM16jtes56S+ZudrcbUL6/ACgIULF2L27NkYN24cJkyYgGXLlsHlcuHOO++UujQiIpJI2IfXjTfeiOrqajz88MOwWCwYPXo0Pv/88zMGcRARUewI+/ACgLlz5561m1BKarUaS5YsOaObks6O16zneM16jtes5yLtmoX1c15ERERd4WKUREQUcRheREQUcRheREQUcRhep1m+fDlyc3Oh0WgwceJEbN++/ZzH22w2zJkzBxkZGVCr1Rg0aBA+/fTTTseUl5fj1ltvRUpKCuLj4zFixAjs3Bk9k0kF+pp5vV489NBDyMvLQ3x8PPr3748//elP3ZoyJlL05Jr99Kc/hUwmO+N19dVX+48RQuDhhx9GRkYG4uPjMXXqVBw5ciQUXyVkAnnNWlpa8MADD2DEiBFITExEZmYmbr/9dlRUVITq64REoP+edfSb3/wGMpkMy5YtC1L15yHIb82aNUKlUonXX39dHDhwQNxzzz3CaDQKq9Xa5fEej0eMGzdOXHXVVeKHH34QxcXF4ttvvxWFhYX+Y+rq6kROTo644447xLZt28Tx48fFF198IY4ePRqqrxVUwbhmTzzxhEhJSREff/yxKC4uFmvXrhVarVa88MILofpaQdXTa1ZbWysqKyv9r/3794u4uDixcuVK/zFPPfWUMBgM4v333xd79uwR11xzjcjLyxONjY0h+lbBFehrZrPZxNSpU8Xbb78tDh06JLZs2SImTJggxo4dG8JvFVzB+HvWbt26dWLUqFEiMzNTPP/888H9ImfB8OpgwoQJYs6cOf6fvV6vyMzMFEuXLu3y+BUrVoh+/fqJ5ubms57zgQceEFOmTAl4reEiGNfs6quvFr/61a86bZs5c6a45ZZbAlO0xHp6zU73/PPPC51OJ5xOpxBCCJ/PJ8xms3jmmWf8x9hsNqFWq8U//vGPwBYvkUBfs65s375dABAlJSUXXG84CNY1O3nypOjTp4/Yv3+/yMnJkSy82G14SnNzM3bt2oWpU6f6t8nlckydOhVbtmzp8j0ffvghCgoKMGfOHJhMJgwfPhxPPvkkvB1WO/3www8xbtw4XH/99UhPT8eYMWPw6quvBv37hEKwrtnFF1+MDRs24PDhwwCAPXv24IcffsCVV14Z3C8UAr25Zqd77bXXcNNNNyExsW2Rx+LiYlgslk7nNBgMmDhxYrfPGc6Ccc26YrfbIZPJgj6dXCgE65r5fD7cdtttuP/++zFs2LCA190TDK9TzrX8isVi6fI9x48fx7vvvguv14tPP/0UDz30EJ599lk8/vjjnY5ZsWIFBg4ciC+++AK//e1vcd999+Fvf/tbUL9PKATrmj344IO46aabkJ+fD6VSiTFjxmD+/Pm45ZZbgvp9QqE316yj7du3Y//+/bj77rv929rf19tzhrtgXLPTNTU14YEHHsDNN98cFXMhBuua/fnPf4ZCocB9990X0Hp7IyJm2AhXPp8P6enpeOWVVxAXF4exY8eivLwczzzzDJYsWeI/Zty4cXjyyScBAGPGjMH+/fvx8ssvY/bs2VKWL4nuXLN33nkHb775Jt566y0MGzYMhYWFmD9/PjIzM2PymnX02muvYcSIEZgwYYLUpUSM812zlpYW3HDDDRBCYMWKFSGuLjx1dc127dqFF154Abt37w74Ch29wZbXKb1ZfiUjIwODBg1CXIeVZIcMGQKLxYLm5mb/MUOHDu30viFDhqC0tDTA3yD0gnXN7r//fn/ra8SIEbjtttuwYMECLF26NHhfJkR6c83auVwurFmzBnfddVen7e3v6805I0Ewrlm79uAqKSnBl19+GRWtLiA41+z7779HVVUVsrOzoVAooFAoUFJSgj/84Q/Izc0N9Fc4L4bXKR2XX2nXvvxKQUFBl++ZPHkyjh49Cp/P5992+PBhZGRkQKVS+Y8pKirq9L7Dhw8jJycnCN8itIJ1zdxu979X7T0lLi6u03siVW+uWbu1a9fC4/Hg1ltv7bQ9Ly8PZrO50zkdDge2bdt23nNGgmBcM+DfwXXkyBF89dVXSElJCXjtUgnGNbvtttuwd+9eFBYW+l+ZmZm4//778cUXXwTle5yTJMNEwtSaNWuEWq0Wq1atEj/++KO49957hdFoFBaLRQghxG233SYefPBB//GlpaVCp9OJuXPniqKiIvHxxx+L9PR08fjjj/uP2b59u1AoFOKJJ54QR44cEW+++aZISEgQb7zxRsi/XzAE45rNnj1b9OnTxz9Uft26dSI1NVX88Y9/DPn3C4aeXrN2U6ZMETfeeGOX53zqqaeE0WgUH3zwgdi7d6+YPn161A2VD+Q1a25uFtdcc43o27evKCws7DRE3OPxBP37hEIw/p6dTsrRhgyv0/zlL38R2dnZQqVSiQkTJoitW7f691166aVi9uzZnY7fvHmzmDhxolCr1aJfv37iiSeeEK2trZ2O+eijj8Tw4cOFWq0W+fn54pVXXgnFVwmZQF8zh8Mhfv/734vs7Gyh0WhEv379xH/9139FzS8VIXp+zQ4dOiQAiPXr13d5Pp/PJx566CFhMpmEWq0WV1xxhSgqKgrmVwi5QF6z4uJiAaDL1zfffBPkbxI6gf57djopw4uzyhMRUcThPS8iIoo4DC8iIoo4DC8iIoo4DC8iIoo4DC8iIoo4DC8iIoo4DC8iIoo4DC8iIoo4DC8iCgiZTIb3339f6jIoRjC8KCYsX74cubm50Gg0mDhxIrZv3y51SWElEMFTWVkZFQuGUmRgeFHUe/vtt7Fw4UIsWbIEu3fvxqhRozBt2jRUVVVd0HlbWloCVGF0MJvNUKvVUpdBMYLhRVHvueeewz333IM777wTQ4cOxcsvv4yEhAS8/vrrPTqPTCbDihUrcM011yAxMRFPPPEEHnnkEYwePRqvv/46srOzodVq8bvf/Q5erxdPP/00zGYz0tPT8cQTT/jPI4TAI488guzsbKjVamRmZp5zZdo9e/bgsssug06ng16vx9ixY7Fz504AwKpVq2A0GvH+++9j4MCB0Gg0mDZtGsrKyjqdY8WKFejfvz9UKhUGDx6Mv//97/597WsxXXvttZDJZGddm6m5uRlz585FRkYGNBoNcnJyOq2x1rH19sgjj0Amk53xWrVqFYC25TmWLl2KvLw8xMfHY9SoUXj33Xe7+0dBxCVRKLp5PB4RFxcn3nvvvU7bb7/9dnHNNdf4f549e7a49NJLz3kuACI9PV28/vrr4tixY6KkpEQsWbJEaLVacd1114kDBw6IDz/8UKhUKjFt2jQxb948cejQIfH6668LAP4ZvdeuXSv0er349NNPRUlJidi2bds5VxoYNmyYuPXWW8XBgwfF4cOHxTvvvCMKCwuFEEKsXLlSKJVKMW7cOLF582axc+dOMWHCBHHxxRf7379u3TqhVCrF8uXLRVFRkXj22WdFXFyc+Prrr4UQQlRVVQkAYuXKlaKyslJUVVV1WcczzzwjsrKyxHfffSdOnDghvv/+e/HWW291uj7t17mhoaHTMiP/8z//IxISEsS+ffuEEEI8/vjjIj8/X3z++efi2LFjYuXKlUKtVotvv/32nH8GRO0YXhTVysvLBQCxefPmTtvvv/9+MWHCBP/PDz74oLjtttvOeS4AYv78+Z22LVmyRCQkJAiHw+HfNm3aNJGbmyu8Xq9/2+DBg8XSpUuFEEI8++yzYtCgQaK5ublb30Gn04lVq1Z1uW/lypWdglEIIQ4ePCgAiG3btgkhhLj44ovFPffc0+l9119/vbjqqqs6fbfTA/508+bNE5dffrnw+Xxd7j/bObZs2SI0Go14++23hRBCNDU1iYSEhDP+TO666y5x8803n7MGonbsNiQCsHTpUqxevfq8x40bN+6Mbbm5udDpdP6fTSYThg4d2mk1aJPJ5L/Hdv3116OxsRH9+vXDPffcg/feew+tra1n/cyFCxfi7rvvxtSpU/HUU0/h2LFjnfYrFAqMHz/e/3N+fj6MRiMOHjwIADh48CAmT57c6T2TJ0/27++uO+64A4WFhRg8eDDuu+8+rF+//rzvKS0txYwZM7Bo0SLccMMNAICjR4/C7XbjZz/7GbRarf+1evXqM74b0dkwvCiqpaamIi4uDlartdN2q9UKs9nc4/MlJiaesU2pVHb6WSaTdbnN5/MBALKyslBUVISXXnoJ8fHx+N3vfodLLrnkrANAHnnkERw4cABXX301vv76awwdOhTvvfdej2u/UBdddBGKi4vxpz/9CY2Njbjhhhtw3XXXnfV4l8uFa665BgUFBXjsscf8251OJwDgk08+6bSk/I8//sj7XtRtDC+KaiqVCmPHjsWGDRv823w+HzZs2ICCggLJ6oqPj8cvf/lLvPjii/j222+xZcsW7Nu376zHDxo0CAsWLMD69esxc+ZMrFy50r+vtbXVP4ADAIqKimCz2TBkyBAAwJAhQ7Bp06ZO59u0aROGDh3q/1mpVMLr9Z63br1ejxtvvBGvvvoq3n77bfzzn/9EXV3dGccJIXDrrbfC5/Ph73//O2QymX/f0KFDoVarUVpaigEDBnR6ZWVlnbcGIgBQSF0AUbAtXLgQs2fPxrhx4zBhwgQsW7YMLpcLd955p/+YxYsXo7y8vFtdhxdq1apV8Hq9mDhxIhISEvDGG28gPj4eOTk5Zxzb2NiI+++/H9dddx3y8vJw8uRJ7NixA7NmzfIfo1QqMW/ePLz44otQKBSYO3cuJk2ahAkTJgAA7r//ftxwww0YM2YMpk6dio8++gjr1q3DV1995T9Hbm4uNmzYgMmTJ0OtViMpKemMWp577jlkZGRgzJgxkMvlWLt2LcxmM4xG4xnHPvLII/jqq6+wfv16OJ1Of2vLYDBAp9Nh0aJFWLBgAXw+H6ZMmQK73Y5NmzZBr9dj9uzZF3qJKRZIfdONKBT+8pe/iOzsbKFSqcSECRM6DXAQovujDU8fkLBkyRIxatSoM841ffr0TtsuvfRS8fvf/14IIcR7770nJk6cKPR6vUhMTBSTJk0SX331VZef6fF4xE033SSysrKESqUSmZmZYu7cuaKxsVEI0TZgw2AwiH/+85+iX79+Qq1Wi6lTp4qSkpJO53nppZdEv379hFKpFIMGDRKrV6/utP/DDz8UAwYMEAqFQuTk5HRZyyuvvCJGjx4tEhMThV6vF1dccYXYvXt3l9fn0ksvFQDOeK1cuVIIIYTP5xPLli0TgwcPFkqlUqSlpYlp06aJjRs3dvnZRKeTCSGEhNlJRBdg1apVmD9/Pmw2m9SlEIUU73kREVHEYXgREVHEYbchERFFHLa8iIgo4jC8iIgo4jC8iIgo4jC8iIgo4jC8iIgo4jC8iIgo4jC8iIgo4jC8iIgo4jC8iIgo4vx/k5DXlUqeyvwAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1200x400 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "monte_carlo.view_histogram(kde=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Add compensator:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# distance from last surface to image plane\n",
    "tolerancing.add_compensator(\"thickness\", surface_number=2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Rerun Monte Carlo with compensator:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "monte_carlo = MonteCarlo(tolerancing)\n",
    "monte_carlo.run(num_iterations=100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAboAAAGGCAYAAAANT0RTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/TGe4hAAAACXBIWXMAAA9hAAAPYQGoP6dpAABKs0lEQVR4nO3deXxU1cE+8Gcyycwks2XPJGRlC2FHEIwoLlARrIDgrhWs1bYCFXlpFfsWVxrfWq1LEW0r4UcrRVHApSr74sKOkU0ChBBClpksZCYzSSbJzPn9ETJN2JPM5N6Zeb6fdz5vM8vNA0gezr3nnKsQQggQEREFqBCpAxAREfkSi46IiAIai46IiAIai46IiAIai46IiAIai46IiAIai46IiAIai46IiAJaqNQBfM3tdqO0tBR6vR4KhULqOEREdBFCCNTW1iIpKQkhId4bhwV80ZWWliIlJUXqGEREdIWKi4uRnJzsteMFfNHp9XoALb9xBoNB4jRERHQxNpsNKSkpnp/b3hLwRdd6utJgMLDoiIj8gLcvM3EyChERBTQWHRERBTQWHRERBTQWHRERBTQWHRERBTQWHRERBTQWHRERBTRJi27x4sUYPHiwZ41bdnY2vvzyS8/rN954IxQKRbvHr371KwkTExGRv5F0wXhycjJefvll9OnTB0II/L//9/8wefJkfP/99xgwYAAA4NFHH8ULL7zg+UxERIRUcYmIyA9JWnS33357u68XLlyIxYsXY8eOHZ6ii4iIgMlkkiIeEREFANlco3O5XFixYgUcDgeys7M9z7///vuIjY3FwIEDMX/+fNTV1V3yOE6nEzabrd2DiIiCl+R7XR44cADZ2dloaGiATqfD6tWr0b9/fwDA/fffj7S0NCQlJWH//v146qmnkJ+fj1WrVl30eDk5OXj++ee7Kz4REcmcQgghpAzQ2NiIU6dOwWq14qOPPsI//vEPbN261VN2bW3atAljx47F8ePH0atXrwsez+l0wul0er5u3Q3barVyU+cr9P6BH7C3tARalQqzR16DeK1O6khEFARsNhuMRqPXf15LXnTnGjduHHr16oV33333vNccDgd0Oh2++uorjB8//oqO56vfuEC16sdDmLf+K2RERqGizoGMyCismHYPwsPCpI5GRAHOVz+vZXONrpXb7W43ImsrLy8PAJCYmNiNiYLHQYsZz2xaj9EpqZg/egyeHHUtjlZVYv7GdVJHIyLqNEmv0c2fPx8TJkxAamoqamtrsXz5cmzZsgVr165FQUEBli9fjokTJyImJgb79+/Hk08+iTFjxmDw4MFSxg5Yf921A7HhEXhg4BAoFAqkGiNx38DBWPrD9/jViJHoFxsndUQiog6TdERnsVjw0EMPITMzE2PHjsXu3buxdu1a/OQnP4FKpcKGDRtwyy23oF+/fvif//kfTJs2DZ999pmUkQOWxWHHxsIC3JiegTCl0vP8qB4piA4Px9K8fRKmIyLqPElHdO+9995FX0tJScHWrVu7MU1wW3n4EJQhIRjVI7nd86EhIbgxLQNr8n/EvGuvRywX7BORn5HdNTrqfm4hsOLgfoxI7IGIMNV5r1+fmg4A+ODQ/m5ORkTUdSw6wn5zOUpqbRidknrB13UqFYaaEvHFsaPdnIyIqOtYdIRvi4sQHhqGXlHRF33PkHgTfqysQGktd5ohIv/CoiN8c6oImTExUIZc/D+HgfHxUCoU2FR4ohuTERF1HYsuyNU3NWFfWelllw5EhKnQJyYWG04UdFMyIiLvYNEFuT2lJWhyu5F1BWvkhsSbsP30KdgbG7shGRGRd7Dogty3xUWI1GiQqNNf9r2DE0xocrux83RxNyQjIvIOFl2Q++50MTJjYqFQKC773riICERrwrG79HQ3JCMi8g4WXRBzNjfjSGUFekZefLZlWwqFAr2io7GrhEVHRP6DRRfEjlVXodntRlpk5BV/pk90DA5YLKhravJdMCIiL2LRBbGDFjNCFAokd+B2GH2jY+ESbuSVl/kwGRGR97DogtgBixlJej3Uyivf8jRRr4c2TMXTl0TkN1h0QeyAxYwUg7FDnwlRKNA7OpoTUojIb7DoglSjy4X8ygqkGSM7/NneUTHIKy+Dy+32fjAiIi9j0QWpY1WVaHK7O1V06ZGRqG9uxokzZ7wfjIjIy1h0Qap1IkqKsWOnLgEg9Ww5HrCUezkVEZH3seiCVH5VJeK12g5NRGkVERYGk1aHAxazD5IREXkXiy5IHauuhukKtv26mFRjJPabOaIjIvlj0QWp49VVV7S/5cWkGY04XFGBZk5IISKZY9EFoVqnE2aHHYk6XaePkRYZBaerGcerq7yYjIjI+1h0QejEmWoAQJLuyndEOVeqwQgFwOt0RCR7LLogdPxs0Zm6MKILDwuDSafHQRYdEckciy4IHa+uQmx4BNShHZ9x2VYPvQFHKiu8lIqIyDdYdEHoeHVVl2Zctko2GHCkshJCCC+kIiLyDRZdEDpWXdWl05atkg0G1DY6UWav9UIqIiLfYNEFGWdzM07bbF1aWtAqWd+yq8qRysouH4uIyFdYdEGmyFoDtxBeGdFFh4cjPDQM+VW8TkdE8sWiCzJFNTUAgHht14tOcfamrRzREZGcseiCTJG1BmqlEka12ivH66E34MdKi1eORUTkCyy6IHPKWoM4rRYKhcIrx0s2GFB45gyczc1eOR4Rkbex6IJMkbUGseERXjteD70BLiFQcHYROhGR3LDogszJmpYRnbck6Vu2ETvGPS+JSKZYdEGk2e1GaW0t4iO6PhGlVURYGKI14dzcmYhkS9KiW7x4MQYPHgyDwQCDwYDs7Gx8+eWXntcbGhowc+ZMxMTEQKfTYdq0aTCbubdiZ5XW2uASbq+O6ADApNPjWBWLjojkSdKiS05Oxssvv4y9e/diz549uPnmmzF58mQcOnQIAPDkk0/is88+w8qVK7F161aUlpZi6tSpUkb2a0XWGgBAXIR3iy5Rr+OpSyKSra7t6ttFt99+e7uvFy5ciMWLF2PHjh1ITk7Ge++9h+XLl+Pmm28GAOTm5iIrKws7duzANddcI0Vkv3bKakWIQoHo8HCvHjdJZ8CWk4VwNjd3eaNoIiJvk801OpfLhRUrVsDhcCA7Oxt79+5FU1MTxo0b53lPv379kJqaiu3bt0uY1H+dstYgNiICoSHe/WNP1OvhEgInz44YiYjkRPJ/fh84cADZ2dloaGiATqfD6tWr0b9/f+Tl5UGlUiEyMrLd+xMSElBeXn7R4zmdTjidTs/XNpvNV9H9zsmaM4gN9+5pSwCefTOPV1UhMybW68cnIuoKyUd0mZmZyMvLw86dO/HrX/8a06dPx+HDhzt9vJycHBiNRs8jJSXFi2n92ymrFXFa762ha6VTqWBUa3idjohkSfKiU6lU6N27N4YPH46cnBwMGTIEb7zxBkwmExobG1Fzdm/GVmazGSaT6aLHmz9/PqxWq+dRXFzs41+B/yiptSHGi4vF20rU67nEgIhkSfKiO5fb7YbT6cTw4cMRFhaGjRs3el7Lz8/HqVOnkJ2dfdHPq9Vqz3KF1gcBNqcT9sZG3xWdVof8Km7uTETyI+k1uvnz52PChAlITU1FbW0tli9fji1btmDt2rUwGo145JFHMHfuXERHR8NgMGD27NnIzs7mjMtOKKltuVYZE+GbojPp9Pi6uAjNbrfXJ7sQEXWFpEVnsVjw0EMPoaysDEajEYMHD8batWvxk5/8BADwl7/8BSEhIZg2bRqcTifGjx+Pt99+W8rIfqv07KScGC8vLWhl0unQ7HbjtM2K9Mgon3wPIqLOkLTo3nvvvUu+rtFosGjRIixatKibEgWuklobQkNCYFBrfHL81hu5njhzhkVHRLLCc0xBomUiSjhCvHR7nnNFasKhUYbiBO9iQEQyw6ILEiU2G6J8NBEFAEIUCiTodCw6IpIdFl2QOG2z+uz6XKsErQ7HWXREJDMsuiBRUmtDtA9HdAA4oiMiWWLRBQFnczOq6ut9PqIzaXWorq+HtaHBp9+HiKgjWHRBwLOGzscjuv/OvOSojojkg0UXBLqr6BK0Z4uu5oxPvw8RUUew6IJAqc0GBYAoH5+6VIeGIjY8AgXVHNERkXyw6IJAqb0WkZrwbtmaK16rxUmO6IhIRlh0QaDMXuvz0VyreC1nXhKRvLDogkBZrR1RGt9s/XWuBK0ORdYauIXolu9HRHQ5LLogUFprQ5Sme0Z0CTodnC4Xyuy13fL9iIguh0UX4IQQKLfXduuIDgAKz/A6HRHJA4suwNU2NqK+ubnbrtHFhLdMeuGEFCKSCxZdgCs/ewqxu05dKkNCEBehRSGLjohkgkUX4MpqzxZdN43ogJYlBjx1SURywaILcOX2WigAGH10w9ULSeASAyKSERZdgCuz27ttsXirBK0OJbU2NLlc3fY9iYguhkUX4LpzxmWrBK0WLiFQbLN26/clIroQFl2Aa9n+q3uLLr51iQEnpBCRDLDoAlxZbW23zbhsFanRQK1U4mRNTbd+XyKiC2HRBbhyu71bZ1wCgEKhQLxWx7V0RCQLLLoAVut0wtHU2O0jOgBcS0dEssGiC2Bmhx0Auv0aHdAyIeUk19IRkQyw6AJYub2l6KQY0cVrdSiz18LZ3Nzt35uIqC0WXQCzSDiii9dqIQCcsnKJARFJi0UXwMwOO3QqFcKUym7/3q1LDDghhYikxqILYGa7XZLRHAAY1WpoQkNx0sqiIyJpsegCmNnh6NY9LttqWWLAzZ2JSHosugBWLsGuKG1xiQERyQGLLoCZHdKdugRaNndm0RGR1Fh0AcotBCocdYhUd//SglbxWi0sDgfqm5oky0BExKILUFX1dXAJt+QjOgAostZIloGIiEUXoMx26dbQtYrXagGw6IhIWpIWXU5ODq6++mro9XrEx8djypQpyM/Pb/eeG2+8EQqFot3jV7/6lUSJ/YeU23+10qvUCA8N5Vo6IpKUpEW3detWzJw5Ezt27MD69evR1NSEW265BQ6Ho937Hn30UZSVlXkef/rTnyRK7D/MdjtCFAoYJFpeALS9i0GNZBmIiEKl/OZfffVVu6+XLl2K+Ph47N27F2PGjPE8HxERAZPJ1N3x/JrZYYdRrUGIQiFpjrgILUd0RCQpWV2js57dFzE6Orrd8++//z5iY2MxcOBAzJ8/H3V1dVLE8ysWh0PS05at4rVcS0dE0pJ0RNeW2+3GnDlzMHr0aAwcONDz/P3334+0tDQkJSVh//79eOqpp5Cfn49Vq1Zd8DhOpxNOp9Pztc1m83l2OSq318KoVksdAwlanWeJQXhYmNRxiCgIyaboZs6ciYMHD+Kbb75p9/xjjz3m+d+DBg1CYmIixo4di4KCAvTq1eu84+Tk5OD555/3eV65M9vtSNDppI7RbuZlv9g4idMQUTCSxanLWbNm4fPPP8fmzZuRnJx8yfeOGjUKAHD8+PELvj5//nxYrVbPo7i42Ot5/YGlzoFICSeitPrvXQxqpA1CREFL0hGdEAKzZ8/G6tWrsWXLFmRkZFz2M3l5eQCAxMTEC76uVquhlsEpOyk1uVyorq+HUQbX6PQqFcJDQ1HEuxgQkUQkLbqZM2di+fLl+OSTT6DX61FeXg4AMBqNCA8PR0FBAZYvX46JEyciJiYG+/fvx5NPPokxY8Zg8ODBUkaXtcqzk3XkMKLjEgMikpqkRbd48WIALYvC28rNzcWMGTOgUqmwYcMGvP7663A4HEhJScG0adPwv//7vxKk9R+WupZ1iHIY0QHg7XqISFKSn7q8lJSUFGzdurWb0gSOirO7okh1L7pzxWu12FVSInUMIgpSspiMQt5lcTgQolBAL5NrlfEROlTUOVDHuxgQkQRYdAGoZVcUteS7orTi5s5EJCUWXQCqcDhkc30O+O/tergVGBFJgUUXgCwOBwwyOW0JADqVChFhYSw6IpIEiy4AtW7oLBcKhQIJXGJARBJh0QUgi8Mhq6IDWu5iwM2diUgKLLoA43K7UVVfJ6trdEDLhBSeuiQiKbDoAkx1fT3cQshiV5S24rVaVNbVwdHYKHUUIgoyLLoAY2ldLC67EV3LzEsuMSCi7saiCzCe7b9kNqLjEgMikgqLLsBYHA4oAFktLwBalhjoVCpOSCGibseiCzAVDgf0ajVCQ+T3R8u7GBCRFOT305C6xCKzNXRtxUdoceJMtdQxiCjIsOgCTEVdnexOW7ZK0HItHRF1PxZdgDE77PItOp0ONQ0NsDY0SB2FiIIIiy7AVDgcsltD1yqeMy+JSAIsugAihEBlnQMGmRZdwtnb9fD0JRF1JxZdAKltbITT5ZLdYvFWmtAwRGo0LDoi6lYsugBS0borikyv0QEtMy9ZdETUnVh0AcTikOeuKG3Fa3UoPMOiI6Luw6ILIBWt239p5Duia7kv3RkIIaSOQkRBgkUXQCwOBzTKUGhCw6SOclEJOi0cTU2orKuTOgoRBQkWXQCprHPIejQH/HdzZ+6QQkTdhUUXQCwO+S4taBUXoUWIQoETnJBCRN2ERRdAKuocst0VpVWYUom4iAgUckRHRN2ERRdAzHb5bujcVpxWhxOceUlE3YRFF0Aq6+pku1i8LZNWx2t0RNRtWHQBotHlQo2zQdaLxVslaHUottnQ5HJJHYWIggCLLkBUnl1DJ/fJKEDLXQxcwo1im1XqKEQUBFh0AaLCsyuKf4zoAHCHFCLqFp0quhMnTng7B3WRZ1cUPxjRRWk0UCuV3POSiLpFp4qud+/euOmmm/Cvf/0LDbyJpixU1NUhRKGA3g9GdAqFAiadHgWckEJE3aBTRbdv3z4MHjwYc+fOhclkwi9/+Uvs2rXL29moAyxn7yweolBIHeWKJGh1OF5dJXUMIgoCnSq6oUOH4o033kBpaSmWLFmCsrIyXHfddRg4cCBee+01VFRUeDsnXUZFXZ3sF4u3ZdJxiQERdY8uTUYJDQ3F1KlTsXLlSvzf//0fjh8/jnnz5iElJQUPPfQQysrKLvn5nJwcXH311dDr9YiPj8eUKVOQn5/f7j0NDQ2YOXMmYmJioNPpMG3aNJjN5q7EDkiVfrD9V1smnR5nGhpQXc/NnYnIt7pUdHv27MHjjz+OxMREvPbaa5g3bx4KCgqwfv16lJaWYvLkyZf8/NatWzFz5kzs2LED69evR1NTE2655RY4zs4gBIAnn3wSn332GVauXImtW7eitLQUU6dO7UrsgGR22P1ixmWrRF3r5s6ckEJEvhXamQ+99tpryM3NRX5+PiZOnIhly5Zh4sSJCAlp6c2MjAwsXboU6enplzzOV1991e7rpUuXIj4+Hnv37sWYMWNgtVrx3nvvYfny5bj55psBALm5ucjKysKOHTtwzTXXdCZ+QKpwOJCkN0gd44olaHVQACg4U40RST2kjkNEAaxTRbd48WL8/Oc/x4wZM5CYmHjB98THx+O9997r0HGt1pYFxNHR0QCAvXv3oqmpCePGjfO8p1+/fkhNTcX27dsvWHROpxNOp9Pztc1m61AGfySEQGV9nV+N6MKUSsRptbxOR0Q+16miW79+PVJTUz0juFZCCBQXFyM1NRUqlQrTp0+/4mO63W7MmTMHo0ePxsCBAwEA5eXlUKlUiIyMbPfehIQElJeXX/A4OTk5eP755zv2C/JzNqcTjS6XX6yhaytBq0NBNYuOiHyrU9foevXqhcrKyvOer66uRkZGRqeCzJw5EwcPHsSKFSs69flW8+fPh9Vq9TyKi4u7dDx/UOHZ/st/RnRAy8zLY1xiQEQ+1qkRnRDigs/b7XZoOrF7/qxZs/D5559j27ZtSE5O9jxvMpnQ2NiImpqadqM6s9kMk8l0wWOp1Wqo/ewHfle1bv8V6Qd3LmjLpNVjY+EJOJuboQ7t1H+KRESX1aGfLnPnzgXQsrPFggULEBER4XnN5XJh586dGDp06BUfTwiB2bNnY/Xq1diyZct5o8Hhw4cjLCwMGzduxLRp0wAA+fn5OHXqFLKzszsSPaBZ/GhD57YS9Xq4hUBhzRn0i42TOg4RBagOFd33338PoKWgDhw4AJVK5XlNpVJhyJAhmDdv3hUfb+bMmVi+fDk++eQT6PV6z3U3o9GI8PBwGI1GPPLII5g7dy6io6NhMBgwe/ZsZGdnc8ZlGxUOB9TKUGj8bFSUqNMDAAqqq1l0ROQzHfrJuHnzZgDAww8/jDfeeAMGQ9emsy9evBgAcOONN7Z7Pjc3FzNmzAAA/OUvf0FISAimTZsGp9OJ8ePH4+233+7S9w00lXUOvzttCQA6lQpGtZrX6YjIpzo1BMjNzfXKN7/Ytb62NBoNFi1ahEWLFnnlewYii8PhdxNRWiXq9Nzzkoh86oqLburUqVi6dCkMBsNldyZZtWpVl4PRlauo89+iM+n0OFp9/gxeIiJvueKiMxqNUJzdGd9oNPosEHWcxeFAchdPI0slSa/HN8VFaHK5EKZUSh2HiALQFRdd29OV3jp1Sd5R4XBgQFy81DE6JUmnR7PbjSJrDXpHx0gdh4gCUKcWjNfX16Ou7r+7zhcVFeH111/HunXrvBaMrkyjy4UaZ4NfTkYBgMSz+3NyQgoR+Uqnim7y5MlYtmwZAKCmpgYjR47Eq6++ismTJ3tmUlL3qPTTNXSt9CoV9CoVJ6QQkc90+g7j119/PQDgo48+gslkQlFREZYtW4Y333zTqwHp0lp3RfGnDZ3bUigUSNTpOaIjIp/pVNHV1dVBr29Z7Ltu3TpMnToVISEhuOaaa1BUVOTVgHRprftc+tuGzm0l6vU4eoG9U4mIvKFTRde7d2+sWbMGxcXFWLt2LW655RYAgMVi6fIicuoYi8OBEIUCej8d0QFAks6AEzVn0ORySR2FiAJQp4puwYIFmDdvHtLT0zFq1CjPvpPr1q3DsGHDvBqQLq2izgGjWo2Qs0s//FEPgwHNbjcKa3i3cSLyvk7tjHLnnXfiuuuuQ1lZGYYMGeJ5fuzYsbjjjju8Fo4ur2VXFP89bQkAPc7OvMyvqkTfmFiJ0xBRoOn0LsAmk+m8W+WMHDmyy4GoYyr8ePuvVjqVCpEaDY5W8TodEXlfp4rO4XDg5ZdfxsaNG2GxWOB2u9u9fuLECa+Eo8uzOOx+X3RAy6gunxNSiMgHOlV0v/jFL7B161b87Gc/Q2JiomdrMOp+ljoHkg3+vyVbD70BhyosUscgogDUqaL78ssv8Z///AejR4/2dh7qACEEKuvq/HppQaseegPWnTgOe2MjdG3uc0hE1FWdmnUZFRWF6Ohob2ehDqppaECz2w2jJgBOXZ5dlnKM1+mIyMs6VXQvvvgiFixY0G6/S+p+gbBYvFWiTo8QhQL5LDoi8rJOnbp89dVXUVBQgISEBKSnpyMsLKzd6/v27fNKOLo0iyNwik6lVMKk1eFIZYXUUYgowHSq6KZMmeLlGNQZng2dA+DUJdBy+vIwi46IvKxTRffss896Owd1gsXhQHhoGNTKTi+HlJUUgxFrC47BLYRf7/RCRPLSqWt0QMvtef7xj39g/vz5qK6uBtByyrKkpMRr4ejSLA6H396H7kJSDEY4mppQbLVKHYWIAkinhgL79+/HuHHjYDQacfLkSTz66KOIjo7GqlWrcOrUKc+96si3LHWBsVi8VaqxZT3g4UoL0iIjpQ1DRAGjUyO6uXPnYsaMGTh27Bg0bUYUEydOxLZt27wWji7NbLcHxESUVga1BpEaDQ5z4TgReVGnim737t345S9/ed7zPXr0QHl5eZdD0ZUJtFOXAJBsMLLoiMirOlV0arUaNpvtvOePHj2KuLi4LoeiK9N6i55AkmIwciswIvKqThXdpEmT8MILL6CpqQkAoFAocOrUKTz11FOYNm2aVwPShTkaG1HX1ARjgI3oUg1GWBwOVHEzAiLykk4V3auvvgq73Y64uDjU19fjhhtuQO/evaHX67Fw4UJvZ6QLsATQrihttU5I4aiOiLylU7MujUYj1q9fj2+//RY//PAD7HY7rrrqKowbN87b+egiKlp3RQmwEV1chBbasDDsN5djTFq61HGIKAB0uOjcbjeWLl2KVatW4eTJk1AoFMjIyIDJZIIQgrfs6SZmhx1A4I3oFAoF0oyR2G/mpCYi8o4OnboUQmDSpEn4xS9+gZKSEgwaNAgDBgxAUVERZsyYgTvuuMNXOekcFocDKqUS4aGBsStKW2nGSPzAoiMiL+nQT8mlS5di27Zt2LhxI2666aZ2r23atAlTpkzBsmXL8NBDD3k1JJ2vwmFHlEYTkCPo9MgofFlwDGa7HQk6ndRxiMjPdWhE9+9//xvPPPPMeSUHADfffDOefvppvP/++14LRxdncThgCLDTlq3Sz+6KcsDCUR0RdV2Him7//v249dZbL/r6hAkT8MMPP3Q5FF2exRFY23+1FaUJh1GtxgGLWeooRBQAOlR01dXVSEhIuOjrCQkJOHPmTJdD0eWZHQ5EBuiIrnVCyg/cZYeIvKBDRedyuRB6ickPSqUSzc3NV3y8bdu24fbbb0dSUhIUCgXWrFnT7vUZM2ZAoVC0e1xqRBlMLA5HwC0taCs9Mgo/mMsghJA6ChH5uQ5NRhFCYMaMGVBf5JSZ0+ns0Dd3OBwYMmQIfv7zn2Pq1KkXfM+tt96K3Nxcz9cX+97BpKG5CbWNzoBbWtBWr6hofHr0CAprzqBnVLTUcYjIj3Wo6KZPn37Z93RkxuWECRMwYcKES75HrVbDZDJd8TGDQYWjZXusQB7RZURGQQFgT2kJi46IuqRDRdd2ZNVdtmzZgvj4eERFReHmm2/GSy+9hJiYmG7PISeti8UD7c4FbYWHhSHZYMS+slLcPWCQ1HGIyI/JerXxrbfeiqlTpyIjIwMFBQV45plnMGHCBGzfvh1KpfKCn3E6ne1OoV7oLgv+ztJadAF86hIAekZFYU8Z71hPRF0j66K79957Pf970KBBGDx4MHr16oUtW7Zg7NixF/xMTk4Onn/++e6KKAnz2V1RIsLCpI7iU72iorG16CRqGuoRqQmXOg4R+alO3b1AKj179kRsbCyOHz9+0ffMnz8fVqvV8yguLu7GhN3DbK8N2F1R2uod1XKK+vvyMomTEJE/k/WI7lynT59GVVUVEhMTL/oetVod8DMzzQ5HQM+4bBUbEQGjWoM9pSW4Kb2n1HGIyE9JWnR2u73d6KywsBB5eXmIjo5GdHQ0nn/+eUybNg0mkwkFBQX43e9+h969e2P8+PESppae2W4P6BmXrRQKBXpHR2PH6cAblRNR95H01OWePXswbNgwDBs2DAAwd+5cDBs2DAsWLIBSqcT+/fsxadIk9O3bF4888giGDx+Or7/+OuBHbJdjdtgDesZlW/1i4rDfbIajsVHqKETkpyQd0d14442X3Pli7dq13ZjGf1gcdoxITJI6RrfIjI2FS7ixp7QEN6RnSB2HiPyQX01GIcDe2AhHUxOMQTIL0aTVIVKtwfYSnr4kos5h0fkZSxAsFm9LoVCgb0wsthefkjoKEfkpFp2fMdtbii4qSIoOADJjYnGowgJbB/dSJSICWHR+x+xwAEBQLC9olRUbB7cQ2H6aozoi6jgWnZ8xO2oRERYG9SVulxRo4rRamHQ6bD1ZKHUUIvJDLDo/Y3Y4EBUkE1HaGhiXgE0nT/D+dETUYSw6P2Ox22EMwnWEg+ITYHE4cKSqUuooRORnWHR+psxeGxS7opyrT3QM1Eoltpw8IXUUIvIzLDo/Y7bbg3In/zClElmxcdhUyKIjoo5h0fkRtxBnr9EF34gOAAYnmPB9eRkqzs48JSK6Eiw6P1JVXweXcAflZBQAGGZq2fbsq4JjEichIn/CovMj5Z7F4sFZdDqVClmxcfjiWL7UUYjIj7Do/Eh5bS0AICo8OE9dAsDwxCTsKjnN05dEdMVYdH6k3GFHaEgIdKrgW17QapgpCSEKBb48flTqKETkJ1h0fqTcXosoTThCFAqpo0hGp1JhQHwCVh4+KHUUIvITLDo/Um4PnhuuXsr1KWk4VGHBIYtZ6ihE5AdYdH6k3F7LokPLLilRmnB8cOiA1FGIyA+w6PxIWW1t0M64bEsZEoLs5BSsOfIj6puapI5DRDLHovMTQgiUO3jqstX1qWmoa27itToiuiwWnZ+obXSiobmZI7qzYiO0uDqpB97duxtNLpfUcYhIxlh0fqKsdbF4OIuu1YRefVFmr8WnR49IHYWIZIxF5ydaF4vz1OV/9TAYMDTBhL/u2oFGjuqI6CJYdH6i3F4LBQCjmkXX1uTMLBTbrFj2w/dSRyEimWLR+Ymys2voQkP4R9ZWssGIG9LS8cbO72Bx2KWOQ0QyxJ+afqLUbkNUeITUMWRpct8sKBQK/GHzBgghpI5DRDLDovMTpTZb0N6H7nK0KhWmDx6G9ScKsGw/T2ESUXssOj9RWluLGI7oLmqoKRFjM3rij19vxfbiU1LHISIZYdH5ASEEyuzcFeVy7swaiL4xsXjks9XYXXpa6jhEJBMsOj9QXV8Pp8uFaK6hu6TQkBD8esRIpBkj8dDqj/ARd00hIrDo/EKpvWUNHU9dXp5aGYrfjMzG1T2S8bsNazHzi09RbLVKHYuIJMSi8wNltTYA3BXlSqmUSkwfPAy/GDYcO04XY9w/l2Deui+RV17GWZlEQShU6gB0eaW1tQgLCYFepZI6il8Z1SMFQxMSsflkIbYUFWLVkcNIMRgwoU8mxvfsjaGmRCiC+Ca2RMGCRecHSmttiAmP4A/lTlCHhuLW3n1wS6/eyK+qxO7S01hxcD/+tnc3kvR6TM7Mwv0Dh6CHwSB1VCLyERadHyiz1yKSpy27JEShQFZsHLJi4+AWAseqq7Cr5DSW5n2Pd/fuxu19MzH3muuQYjRKHZWIvEzSa3Tbtm3D7bffjqSkJCgUCqxZs6bd60IILFiwAImJiQgPD8e4ceNw7NgxacJK6LTNhmguLfCaEIUCmTGx+NngofjTuPG4d8AgbC06ibH/XIK3d+9Es9stdUQi8iJJi87hcGDIkCFYtGjRBV//05/+hDfffBPvvPMOdu7cCa1Wi/Hjx6OhoaGbk0qrtLaWSwt8RBMaipvSe+Klm8ZhXEYvvLbjW9z38QeoqquTOhoReYmkpy4nTJiACRMmXPA1IQRef/11/O///i8mT54MAFi2bBkSEhKwZs0a3Hvvvd0ZVTJNLhcq6xwsOh9TK0MxLWsABieY8M7eXZj64fvInTwNPaOipY5GRF0k2+UFhYWFKC8vx7hx4zzPGY1GjBo1Ctu3b7/o55xOJ2w2W7uHPyu32yEARHMNXbfoEx2D+aPHwOUWuP/jD1FUUyN1JCLqItkWXXl5OQAgISGh3fMJCQme1y4kJycHRqPR80hJSfFpTl87bWtZ7MzF4t0nNkKL/8kejZAQBe5f9SHMdt7+h8ifybboOmv+/PmwWq2eR3FxsdSRuuT02cXiMTx12a0Mag2eHDUaDc3N+OXna+BsbpY6EhF1kmyLzmQyAQDMZnO7581ms+e1C1Gr1TAYDO0e/uy0zYooTTjClEqpowSd6PBwPD5iJI5UVuL3m9ZLHYeIOkm2RZeRkQGTyYSNGzd6nrPZbNi5cyeys7MlTNa9TttsHM1JKD0yCg8OHoJVRw5jzZEfpY5DRJ0g6axLu92O48ePe74uLCxEXl4eoqOjkZqaijlz5uCll15Cnz59kJGRgT/84Q9ISkrClClTpAvdzYqtVsRE8PqclLKTU3GowoIFmzdgeGISF5UT+RlJR3R79uzBsGHDMGzYMADA3LlzMWzYMCxYsAAA8Lvf/Q6zZ8/GY489hquvvhp2ux1fffUVNEF0p+1im5UTUWTg/oFDoAkLxfyNa7kxNJGfUYgA/1trs9lgNBphtVr97npdo8uFrEWv42eDh+L61HSp4wS9gxYz3ti1Ha/85FZMyxogdRyigOOrn9eyvUZHQFltLQRapruT9AbGJ2BUj2S8tG0zquu5cwqRv2DRyVgx19DJzj39B6HZ7car27+VOgoRXSEWnYyV2KwIUSi4/ZeM6NVq/LRPP6w4uB+HKyxSxyGiK8Cik7HTtTZEacIRGsI/Jjm5MT0DJp0eL2zdzIkpRH6AP0FljGvo5Ck0JAR3Zg3ArtLT2FJUKHUcIroMFp2MnbLWcDNnmRoUn4DMmFj83zfb4OL964hkjUUnY0XWGsRrOeNSjhQKBab264+j1VX4JJ87phDJGYtOpuyNjaiur0cclxbIVs+oaAwzJeKNnd+hyeWSOg4RXQSLTqZOWWsAgCM6mZvUtx+KbTZ8/OMhqaMQ0UWw6GSq6GzRxbHoZC3ZYMSIpB54c9d23sqHSKZYdDJVVFODiLAw6MJUUkehy5jUpx/MdjtHdUQyxaKTqSJrDeIitFAoFFJHoctI1OsxIqkH3t6zk9fqiGSIRSdTRTU1nIjiRyb27ovS2lqsPnJY6ihEdA4WnUydtJ7hRBQ/kmww4ipTEhbt3olmrqsjkhUWnQw5m5thtts5ovMzE/v0RbHNii+O5UsdhYjaYNHJULHNCgEuLfA3acZIDIxLwKLdO+HmHphEssGik6GimhoA4IjOD03o3QfHqquwqbBA6ihEdBaLToZOWmugUiph1GikjkId1DcmFr2jY/D27p28swGRTLDoZOjEmWqYtDqEcGmBX5rQqw/yzOXYVXJa6ihEBBadLBVUVyNBp5M6BnXSoPgEJBsMWLxnl9RRiAgsOlkqOFMFk5ZF568UCgVu7dUH206dxCGLWeo4REGPRScz1oYGVNXXw6TTSx2FumBEYg/ERkTgnb27pY5CFPRYdDJz4kw1AMDEU5d+TRkSglt69saXx496ZtESkTRYdDJzouYMAPAaXQAYnZIGnUqFv+3jqI5ISiw6mSmorkZseATUylCpo1AXqZRK3JzeEx//eAgVDofUcYiCFotOZgrOVHE0F0BuSs+AUqHAkry9UkchClosOpk5Xl3NGZcBJCJMhTFpGfjX/jzYnA1SxyEKSiw6GWlyuVBss3LGZYD5SUYvNLpc+Of+PKmjEAUlFp2MFFlr0Ox2c8ZlgDFqNBidkoYl3+9FfVOT1HGIgg6LTkbyKysBAD0MBomTkLfd0qs3rE4nPjh0QOooREGHRScj+VWViNJooFeppY5CXhYXocXIpB54d+9uNLpcUschCiosOhk5UlmBJD1Hc4FqQu++MDvsWH3ksNRRiIIKi05GjlRVoAeLLmAl6Q0YnpiEt3fvQLPbLXUcoqDBopMJR2MjTttsLLoAN7F3XxTbbPj86BGpoxAFDVkX3XPPPQeFQtHu0a9fP6lj+cTRKk5ECQapxkgMTjDhrV074OKojqhbyLroAGDAgAEoKyvzPL755hupI/lEflUlQhQKJHINXcD7ae9MFNacwX+O5UsdhSgoyH5DxdDQUJhMJqlj+Fx+VSUStDqolEqpo5CPZURFYVB8At7atQO39cmEMkT2/94k8muy/xt27NgxJCUloWfPnnjggQdw6tSpS77f6XTCZrO1e/iDlhmXHM0Fi9v79EPBmWqO6oi6gayLbtSoUVi6dCm++uorLF68GIWFhbj++utRW1t70c/k5OTAaDR6HikpKd2YuHPcQuCQxYI0Y6TUUaibtI7q3ti5nTMwiXxM1kU3YcIE3HXXXRg8eDDGjx+PL774AjU1Nfjwww8v+pn58+fDarV6HsXFxd2YuHNO1pyBvamRRRdkJvfNQmHNGXyS/6PUUYgCmqyL7lyRkZHo27cvjh8/ftH3qNVqGAyGdg+5O2AxAwCLLsikRUZimCkRr+/4jrulEPmQXxWd3W5HQUEBEhMTpY7iVQfMZsRrtdCqVFJHoW42uW8WSmtt3AOTyIdkXXTz5s3D1q1bcfLkSXz33Xe44447oFQqcd9990kdzasOWMqRaoiUOgZJoIfBgGuSU/Dmzu1wNDZKHYcoIMm66E6fPo377rsPmZmZuPvuuxETE4MdO3YgLi5O6mhe43K7cajCgvTISKmjkEQm982C1dmA3Lx9UkchCkiyXke3YsUKqSP4XGHNGdQ1NSGV1+eCVkxEBG5My8C7e3fh3oGDERsRIXUkooAi6xFdMDhgbp2IYpQ4CUnpp30yAQCv7/xO4iREgYdFJ7G95aVI1OkREcaJKMFMq1JhYp9MrDi4H8eqqqSOQxRQWHQS211yGr2jo6WOQTJwU1oGYsMj8OK2zRBCSB2HKGCw6CRU01CPY9VV6B0dI3UUkoEwpRJ39R+Ib4qLsLGwQOo4RAGDRSehvWWlAIA+USw6ajEkwYQBcfF4cdtmOJubpY5DFBBYdBLaU1qCSI2Gs+zIQ6FQ4J4Bg1Bmt+OdvbukjkMUEFh0EtpdWoLe0TFQKBRSRyEZSdTpcUvP3nh7904U1pyROg6R32PRScTZ3IwD5nL0juJEFDrfxD59EakJxx82refEFKIuYtFJZFfpaTS53ciMiZU6CsmQWhmK+wcOxneni/Hh4YNSxyHyayw6iXxddBJRmnD00Mv/7gokjYHxCbg2JRULt21Buf3i92Akoktj0Ulka9FJZMXG8focXdLdWQMRqgzBUxvW8hQmUSex6CRQbq/FseoqDIiLlzoKyZxWpcJDg4fh61NF+Of+PKnjEPklFp0Evj5VBAWA/iw6ugKD4hNwY1oGcr7ZiqNVlVLHIfI7LDoJbCs6ifTIKOh4o1W6Qnf2H4C4CC1mfvEZ6pqapI5D5FdYdN3M2dyMLSdPYGB8gtRRyI+olaF47KqrcdpmxR82b+D1OqIOYNF1s69PnYSjqQkjEpOkjkJ+JlGvx4ODhmL1kcNYtv97qeMQ+Q1Z33g1EH1x7CiS9HokcVkBdcI1ySkostbgpW1b0Ds6BqNT0qSORCR7HNF1I2dzMzacKMBwE0dz1Hl3Zg1Av9g4/Po/nyKfk1OILotF142+OVUEe1Mjhif1kDoK+TFlSAh+edXViNKE4+FPPkZprU3qSESyxqLrRh8ePoAUg5G7oVCXhYeFYfbV16DZ7cYDq1aiwuGQOhKRbLHouonZbsemwhO4PpXXVMg7osLD8eSo0bA6G/DAapYd0cWw6LrJRz8eRGhICEb1SJE6CgWQeK0Wc68Zjao6B+7+aAVKeBqT6Dwsum7gcrvx74P7MSKpByLCwqSOQwEmUafHvOzr4WhqxLQPluNwhUXqSESywqLrBl8eP4rS2lrcmJYhdRQKUPFaLZ6+dgzCw8Jw90crsK7gmNSRiGSDRedjbiHw5q4dGBiXgPTIKKnjUAAzajSYl30d+sXE4Vf/+RSvbf8WzW631LGIJMei87G1BcdwvLoKP+2bKXUUCgKa0FD8cvjVmJKZhbf37MT9qz7kdTsKeiw6H2p0ufDa9m+RFRuHXlHRUsehIBGiUOC2Ppn4n2tGo/BMNW7911J8cHA/98ekoMWi86El3+9FYc0Z3N1/oNRRKAj1jYnFgjE3Y6gpEfM3rce9H3/A2/xQUGLR+UiJzYY3d23Hzek9kWwwSh2HglREWBhmDLkKc68ZjVNWK25bvgzPbdmI6vo6qaMRdRsWnQ80uVx4Yu1/EBEWhkl9+0kdhwhZsXF4dsxNuKNff6w8fAhjlv4Df9nxLWzOBqmjEfkci84H/u/br/FDeRkeu+pqhHPdHMlEmFKJ8b36YOHN43BdShre3bMb1y35O/783TfcVYUCmkIE+BVqm80Go9EIq9UKg8H3e0wu+X4vXvp6C+7pPwjjevby+fcj6ixrQwPWnTiObadOwuV2Y1JmFqYPGcabAkvMLQRqnU7UNTWhye1CiEIBlVIJnUqN8NBQKBQKqSP6jK9+XrPovEQIgdy8fXjp6y0Y36sPpvXrH9D/QVLgqGtqxLaiImwpKkRVfR0GxifgngGD8NM+mTBqNFLHC0hCCJy22fBjpQX5VZUoOFONopoalNbaUF1fD9dFfiyrlaGI00YgxWBEemQUekfHICs2Dv3j4mBQ+/+fVVAX3aJFi/DKK6+gvLwcQ4YMwVtvvYWRI0de0We7o+gampvw7JaNWHn4EG7p2Rt3Zg1gyZHfcbndOGAxY9upkzhUYUFoSAhuSs/AT/v0w43pGdCqVFJH9Fv2xkbklZdhX1kpvi8vxQ/mctQ0tFwf1alUSNTpERehRXR4OIwaDXRhKqhDQxEaEgIhBJrcLtQ3NcPe6MSZhgZU1dfB4nCg3F6LprObAqQajRhmSsJViUkYnpiEzJhYKEP86+pU0BbdBx98gIceegjvvPMORo0ahddffx0rV65Efn4+4uPjL/t5XxadEAJbigrx3JZNKLfX4sFBQ3FtSqpXvweRFGoaGrCrpBh7ykpRWHMGKqUS1yan4sb0DFyXmoaMyCj+Y+4ihBAora3FvvJS7C0twe7SEuRXVcItBLRhYciIikZGZBQyjJFIMRphVGs6/XvpcrthdthRZK3ByZoanKw5g1M2K5rdbmjDwjDsbOldlZiEIQmJMKjVXv7VelfQFt2oUaNw9dVX469//SsAwO12IyUlBbNnz8bTTz992c/74jfO3tiIdQXHkJu3D4cqLMiKjcN9AwcjUaf3yvGJ5KSizoG88jIcsJhxrLoKzW434rVajExKxrDEJAxOSEBWbHxQblguhEBJrQ35lZU4VGHBQYsZeeYyVNa1LN9I0OrQMyoKvaNi0Ds6GiadHiE+/gdCo8uFkzVncPxMNQqqq1Bw5gwcTY1QAOgZFY2hpkQMjI9HVmw8MmNiZXV6OiiLrrGxEREREfjoo48wZcoUz/PTp09HTU0NPvnkk8sewxu/cbVOJ3aXluBwhQU7S4qxp7QETpcLWbFxuKVnbwyIi+e/bikoNDQ34Vh1NfIrK3D8TDVOWWvQ5HZDASDZYESf6BhkREV5bjBs0ukQp9UiOjwCoX52Gg1oGTHVNDSgos6BCocDpbU2lNTWoshag8Iz1ThRcwZ1TU0AAL1KhRRjJNKNkciIjELPqGhZjKDcQsDssOPEmWqcOHMGp6w1OF1r8+yDGhehRa+oaGRERSHZYDj756ZHvFaLmPAI6FSqbvv55quiC/XakXygsrISLpcLCQntZ4ElJCTgyJEjF/yM0+mE0+n0fG21WgG0/AZ21hNffY7NJws9X2vDwpBtSkSUJgLHyspwrKys08cm8kdKAJk6A3pF6FBut+F0rQ1FFjOKLGavfh8FgBAooFAooIACLf/X8kNX0fJlOwKAEICAAM7+fwHALdzwxb/ow0PDkK7TIVGvR6Q6vKUQnI0oNJtRaPbu74U36AD0N0QiU2eAxWFHqb0W5dVVMFdX4TsvHF8JBZ67aSwmZ2Z16vOtP6e9Pf6SddF1Rk5ODp5//vnznk9J8e4NTw979WhE5K8u/E/u4PUQft/lY9TW1sJo9N6OUrIuutjYWCiVSpjP+ZeR2WyGyWS64Gfmz5+PuXPner52u92orq5GTExMt55etNlsSElJQXFxcbes3/M2f87P7NJgdmkEUnYhBGpra5GUlOTV7yProlOpVBg+fDg2btzouUbndruxceNGzJo164KfUavVUJ9zXjwyMtLHSS/OYDD43X98bflzfmaXBrNLI1Cye3Mk10rWRQcAc+fOxfTp0zFixAiMHDkSr7/+OhwOBx5++GGpoxERkR+QfdHdc889qKiowIIFC1BeXo6hQ4fiq6++Om+CChER0YXIvugAYNasWRc9VSlXarUazz777HmnUf2FP+dndmkwuzSY/fJkvY6OiIioq/xvBScREVEHsOiIiCigseiIiCigseg6YNGiRUhPT4dGo8GoUaOwa9eui7536dKlLdsWtXloztk81W63Y9asWUhOTkZ4eDj69++Pd955R/LsAFBTU4OZM2ciMTERarUaffv2xRdffNGlY8ole05ODq6++mro9XrEx8djypQpyM/P94vsbb388stQKBSYM2eOD5L7JntJSQkefPBBxMTEIDw8HIMGDcKePXtkn93lcuEPf/gDMjIyEB4ejl69euHFF1/0+lZVHc1+4403nvdzRqFQ4LbbbvO8RwiBBQsWIDExEeHh4Rg3bhyOHTvm9dzezt7U1ISnnnoKgwYNglarRVJSEh566CGUlpZ2PJigK7JixQqhUqnEkiVLxKFDh8Sjjz4qIiMjhdlsvuD7c3NzhcFgEGVlZZ5HeXl5u/c8+uijolevXmLz5s2isLBQvPvuu0KpVIpPPvlE0uxOp1OMGDFCTJw4UXzzzTeisLBQbNmyReTl5XX6mHLKPn78eJGbmysOHjwo8vLyxMSJE0Vqaqqw2+2yz95q165dIj09XQwePFg88cQTXs3tq+zV1dUiLS1NzJgxQ+zcuVOcOHFCrF27Vhw/flz22RcuXChiYmLE559/LgoLC8XKlSuFTqcTb7zxhqTZq6qq2v2MOXjwoFAqlSI3N9fznpdfflkYjUaxZs0a8cMPP4hJkyaJjIwMUV9fL+vsNTU1Yty4ceKDDz4QR44cEdu3bxcjR44Uw4cP73A2Ft0VGjlypJg5c6bna5fLJZKSkkROTs4F35+bmyuMRuMljzlgwADxwgsvtHvuqquuEr///e+7nLetjmZfvHix6Nmzp2hsbPTaMTvLF9nPZbFYBACxdevWLudty1fZa2trRZ8+fcT69evFDTfc4JOi80X2p556Slx33XVez3ouX2S/7bbbxM9//vN2z02dOlU88MAD3gl9Vlf/Xv3lL38Rer3e8482t9stTCaTeOWVVzzvqampEWq1Wvz73/+WdfYL2bVrlwAgioqKOpSNpy6vQGNjI/bu3Ytx48Z5ngsJCcG4ceOwffv2i37ObrcjLS0NKSkpmDx5Mg4dOtTu9WuvvRaffvopSkpKIITA5s2bcfToUdxyyy2SZv/000+RnZ2NmTNnIiEhAQMHDsQf//hHuFyuTh9TLtkvpPUOF9HR0X6RfebMmbjtttvaHdubfJX9008/xYgRI3DXXXchPj4ew4YNw9///ne/yH7ttddi48aNOHr0KADghx9+wDfffIMJEyZImv1c7733Hu69915otVoAQGFhIcrLy9sd02g0YtSoUZL/Xb1c9guxWq1QKBQd3tbRLxaMS60ztwvKzMzEkiVLMHjwYFitVvz5z3/Gtddei0OHDiE5ORkA8NZbb+Gxxx5DcnIyQkNDERISgr///e8YM2aMpNlPnDiBTZs24YEHHsAXX3yB48eP4/HHH0dTUxOeffbZTh1TLtnP5Xa7MWfOHIwePRoDBw6UffYVK1Zg37592L17t9eydlf2EydOYPHixZg7dy6eeeYZ7N69G7/5zW+gUqkwffp0WWd/+umnYbPZ0K9fPyiVSrhcLixcuBAPPPCAV3J3Nntbu3btwsGDB/Hee+95nisvL/cc49xjtr7mDb7Ifq6GhgY89dRTuO+++zq8pyeLzkeys7ORnZ3t+fraa69FVlYW3n33Xbz44osAWopux44d+PTTT5GWloZt27Zh5syZSEpK8tm/1q+E2+1GfHw8/va3v0GpVGL48OEoKSnBK6+8csGykJOOZp85cyYOHjyIb775RoK07V0ue3FxMZ544gmsX7/+vIlNUruS33e3240RI0bgj3/8IwBg2LBhOHjwIN555x2vFZ2vsn/44Yd4//33sXz5cgwYMAB5eXmYM2cOkpKSJM3e1nvvvYdBgwZh5MiRUkfpsMtlb2pqwt133w0hBBYvXtzh47PorkBnbhd0rrCwMAwbNgzHjx8HANTX1+OZZ57B6tWrPbOMBg8ejLy8PPz5z3/2WtF1JntiYiLCwsKgVCo9z2VlZaG8vByNjY1e+f2QKrtKpfI8P2vWLHz++efYtm2bZ5Qt5+x79+6FxWLBVVdd5Xnd5XJh27Zt+Otf/wqn09nus3LKrlKpkJiYiP79+7f7XFZWFj7++OMuZ/Z19t/+9rd4+umnce+99wIABg0ahKKiIuTk5Hit6Lry98rhcGDFihV44YUX2j3f+jmz2YzExMR2xxw6dKhXcgO+yd6qteSKioqwadOmTt2hgdforkDb2wW1ar1dUNtR26W4XC4cOHDA8x9bU1MTmpqaEBLS/o9AqVTCffYW91JlHz16NI4fP94ux9GjR5GYmAiVSuWV3w+psgMt061nzZqF1atXY9OmTcjIyPBaZl9mHzt2LA4cOIC8vDzPY8SIEXjggQeQl5fnlZLzVfbW95y7jOPo0aNIS0vzSm5fZq+rq5Pl39VWK1euhNPpxIMPPtju+YyMDJhMpnbHtNls2Llzp+R/Vy+XHfhvyR07dgwbNmxATExM5wJ2aOpKEFuxYoVQq9Vi6dKl4vDhw+Kxxx4TkZGRniUDP/vZz8TTTz/tef/zzz8v1q5dKwoKCsTevXvFvffeKzQajTh06JDnPTfccIMYMGCA2Lx5szhx4oTIzc0VGo1GvP3225JmP3XqlNDr9WLWrFkiPz9ffP755yI+Pl689NJLV3xMOWf/9a9/LYxGo9iyZUu76c11dXWyz34uX8269EX2Xbt2idDQULFw4UJx7Ngx8f7774uIiAjxr3/9S/bZp0+fLnr06OFZXrBq1SoRGxsrfve730mavdV1110n7rnnngse8+WXXxaRkZHik08+Efv37xeTJ0/22fICb2ZvbGwUkyZNEsnJySIvL6/d31Wn09mhbCy6DnjrrbdEamqqUKlUYuTIkWLHjh2e12644QYxffp0z9dz5szxvDchIUFMnDhR7Nu3r93xysrKxIwZM0RSUpLQaDQiMzNTvPrqq8LtdkuaXQghvvvuOzFq1CihVqtFz549xcKFC0Vzc/MVH1PO2QFc8NF27ZFcs5/LV0Xnq+yfffaZGDhwoFCr1aJfv37ib3/7m19kt9ls4oknnhCpqalCo9GInj17it///vcd/oHri+xHjhwRAMS6desueDy32y3+8Ic/iISEBKFWq8XYsWNFfn6+13N7O3thYeFF/65u3ry5Q7l49wIiIgpovEZHREQBjUVHREQBjUVHREQBjUVHREQBjUVHREQBjUVHREQBjUVHREQBjUVHREQBjUVHRF6hUCiwZs0aqWMQnYdFR0Fh0aJFSE9Ph0ajwahRo7Br1y6pI8mKN0qqrKzMqzciJfIWFh0FvA8++ABz587Fs88+i3379mHIkCEYP348LBZLl47b1NTkpYSBwWQyQa1WSx2D6DwsOgp4r732Gh599FE8/PDD6N+/P9555x1ERERgyZIlHTqOQqHA4sWLMWnSJGi1WixcuBDPPfcchg4diiVLliA1NRU6nQ6PP/44XC4X/vSnP8FkMiE+Ph4LFy70HEcIgeeeew6pqalQq9VISkrCb37zm4t+3x9++AE33XQT9Ho9DAYDhg8fjj179gAAli5disjISKxZswZ9+vSBRqPB+PHjUVxc3O4YixcvRq9evaBSqZCZmYl//vOfntfS09MBAHfccQcUCoXn63M1NjZi1qxZSExMhEajQVpaGnJyctr9/rSOCp977jkoFIrzHkuXLgXQcguXnJwcZGRkIDw8HEOGDMFHH310pX8URB3T4e2pifyI0+kUSqVSrF69ut3zDz30kJg0aZLn6+nTp4sbbrjhkscCIOLj48WSJUtEQUGBKCoqEs8++6zQ6XTizjvvFIcOHRKffvqpUKlUYvz48WL27NniyJEjYsmSJQKAZyf3lStXCoPBIL744gtRVFQkdu7cecld/AcMGCAefPBB8eOPP4qjR4+KDz/8UOTl5QkhhMjNzRVhYWFixIgR4rvvvhN79uwRI0eOFNdee63n86tWrRJhYWFi0aJFIj8/X7z66qtCqVSKTZs2CSGEsFgsnrs3lJWVCYvFcsEcr7zyikhJSRHbtm0TJ0+eFF9//bVYvnx5u9+f1t/n2tradrdV+fOf/ywiIiLEgQMHhBBCvPTSS6Jfv37iq6++EgUFBSI3N1eo1WqxZcuWS/4ZEHUGi44CWklJiQAgvvvuu3bP//a3vxUjR470fP3000+Ln/3sZ5c8FgAxZ86cds89++yzIiIiQthsNs9z48ePF+np6cLlcnmey8zMFDk5OUIIIV599VXRt29f0djYeEW/Br1eL5YuXXrB13Jzc9uVqBBC/PjjjwKA2LlzpxBCiGuvvVY8+uij7T531113iYkTJ7b7tZ37j4FzzZ49W9x8880XvY3UxY6xfft2odFoxAcffCCEEKKhoUFERESc92fyyCOPiPvuu++SGYg6g6cuiQDk5ORg2bJll33fiBEjznsuPT0der3e83VCQgL69+/f7o7UCQkJnmuCd911F+rr69GzZ088+uijWL16NZqbmy/6PefOnYtf/OIXGDduHF5++WUUFBS0ez00NBRXX3215+t+/fohMjISP/74IwDgxx9/xOjRo9t9ZvTo0Z7Xr9SMGTOQl5eHzMxM/OY3v8G6desu+5lTp05hypQpmDdvHu6++24AwPHjx1FXV4ef/OQn0Ol0nseyZcvO+7UReQOLjgJabGwslEolzGZzu+fNZjNMJlOHj6fVas97LiwsrN3XCoXigs+53W4AQEpKCvLz8/H2228jPDwcjz/+OMaMGXPRyS3PPfccDh06hNtuuw2bNm1C//79sXr16g5n76qrrroKhYWFePHFF1FfX4+7774bd95550Xf73A4MGnSJGRnZ+OFF17wPG+32wEA//nPf5CXl+d5HD58mNfpyCdYdBTQVCoVhg8fjo0bN3qec7vd2LhxI7KzsyXLFR4ejttvvx1vvvkmtmzZgu3bt+PAgQMXfX/fvn3x5JNPYt26dZg6dSpyc3M9rzU3N3smpwBAfn4+ampqkJWVBQDIysrCt99+2+543377Lfr37+/5OiwsDC6X67K5DQYD7rnnHvz973/HBx98gI8//hjV1dXnvU8IgQcffBButxv//Oc/oVAoPK/1798farUap06dQu/evds9UlJSLpuBqKNCpQ5A5Gtz587F9OnTMWLECIwcORKvv/46HA4HHn74Yc975s+fj5KSkis6fdlVS5cuhcvlwqhRoxAREYF//etfCA8PR1pa2nnvra+vx29/+1vceeedyMjIwOnTp7F7925MmzbN856wsDDMnj0bb775JkJDQzFr1ixcc801GDlyJADgt7/9Le6++24MGzYM48aNw2effYZVq1Zhw4YNnmOkp6dj48aNGD16NNRqNaKios7L8tprryExMRHDhg1DSEgIVq5cCZPJhMjIyPPe+9xzz2HDhg1Yt24d7Ha7ZxRnNBqh1+sxb948PPnkk3C73bjuuutgtVrx7bffwmAwYPr06V39LSZqT+qLhETd4a233hKpqalCpVKJkSNHtpu8IcSVz7o8d7LFs88+K4YMGXLesSZPntzuuRtuuEE88cQTQgghVq9eLUaNGiUMBoPQarXimmuuERs2bLjg93Q6neLee+8VKSkpQqVSiaSkJDFr1ixRX18vhGiZjGI0GsXHH38sevbsKdRqtRg3bpwoKipqd5y3335b9OzZU4SFhYm+ffuKZcuWtXv9008/Fb179xahoaEiLS3tgln+9re/iaFDhwqtVisMBoMYO3as2Ldv3wV/f2644QYB4LxHbm6uEEIIt9stXn/9dZGZmSnCwsJEXFycGD9+vNi6desFvzdRVyiEEELCniWiLli6dCnmzJmDmpoaqaMQyRav0RERUUBj0RERUUDjqUsiIgpoHNEREVFAY9EREVFAY9EREVFAY9EREVFAY9EREVFAY9EREVFAY9EREVFAY9EREVFAY9EREVFA+//V/eJyTYrttAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 1200x400 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "monte_carlo.view_histogram(kde=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "By comparing the distributions, we can see that adding compensator generally improved the RMS spot size. The histogram with the compensator shows a distribution with slightly smaller values, indicating that the compensator effectively reduced the impact of the perturbations on the RMS spot size."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": ".venv",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
